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ABSTRACT 


This work constitutes a comprehensive investigation of the application 
of cryogenic technology to coils suitable for MHD applications. Detailed 
background studies are conducted in the specific areas of resistive behavior 
of metals, coil structural requirements, cryogenic refrigerator character- 
istics, properties of coolants, and cooling system design. These studies 
were utilized to develope a series of characteristic curves, Suitable for use 
in the preliminary design phase of this type of cryogenic magnet. The range 
of the investigation covers coils suitable for generating fields of up to 100 KG, 
for use with MHD channels of cross section between .01 and 10 square meters. 


Consideration is given to the optimization of operating temperature and 
coil geometry. A simplified procedure for preliminary coil analysis, uti- 
lizing the characteristic curves, is presented in the form of a flow diagram. 
Example problems indicate that ''minimum loss" coils appear highly attrac- 
tive for use with seeded- gas MHD generators of 50 Mw output, or larger. 

In general, ''optimized'' cryogenic coils should provide at least an order of 
magnitude reduction in total power requirements, compared to ‘conventional’ 
copper coils. A reasonable evaluation of present data indicates that event- 
ually, this may be improved by yet another factor of four. 
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I. INTRODUCTION 


A. Background: 

Recent developments in the field of magnetohydrodynamics point to- 
ward an increasing need for high-field electromagnets combining large 
physical size with low power requirements: In particular, the application 
of seeded-gas MHD generators to commercial power generation (1, 2) in- 
dicates a probable need for fields of the order of 20-950 kilogauss, with 
channel lengths ranging between 1 and 30 meters and widths between 0. 1 
and 5 meters (the larger sizes corresponding to the lower fields). Further, 
recent interest in the application of MHD to marine propulsion (3, 4) sug- 
gests a possible need for fields of up to 200 KG, with channel dimensions 
scaling upwards from the larger figures cited above. It has been shown 
(2,5) that the unprecedented size and high field requirements, coupled with 
the demand for operating economy, preclude the use of conventional 
electromagnets for these applications. 

The startling advance of superconductor technology within the last 
decade has led to the consideration of superconducting magnets, to the 
exclusion of all others, for applications such as the above. Indeed, in- 
vestigation of the Vanadium-Gallium system indicates the possibility that 
“hard'' superconductors may be capable of operation in fields up to 500 
KG (6). Despite this promise, however, superconducting magnets exhibit 
several less-than-desirable characteristics which have hampered their 
development, to date: 

l, The elements utilized in superconducting alloys and compounds 
are relatively rare, and correspondingly expensive. 

2. The mechanical properties of these alloys and compounds--- most 


notably, brittleness--- lead to severe problems in the manufacture and 
eo 





operation of high-field coils. 

3. The present ''standard of the industry''--- ten mil wire--- poses 
serious problems regarding the assembly of physically large coils. 

4. Due to the fact that hard superconductors become relatively poor 
conductors when quenched, the magnetic field of a suddenly quenched super- 
conducting magnet is prone to rapid collapse. This leads, at best, to the 
problem of rapid disposal of a great deal of stored energy, ana can lead, 
at worst, to local arcing and/or destruction of the coil. Subsequently, 
local impairment of the cooling system associated with such a coil--- if 
sufficient to permit a local temperature fluctuation to the vicinity of the 
quench point--- could necessitate immediate and total shutdown. This, 
in itself, could be disastrous, particularly in marine propulsion appli- 
cations. 

Do. To date, the current density attainable in large superconducting 
coils has fallen far below that predicted on the basis of short-sample 
testing, thus adding an additional complication to the design of high-field 
magnets. 

Although significant advances have been made in the search for solu- 
tions to these problems--- particularly, as regards coil protection (7), 
and determination of the factors underlying degraded performance in large 
coils (8-14)--- the successful synthesis of a large, high-field supercon- 
ducting coil has proven to be an elusive, and as yet unrealized, goal. This 
fact has become a critical roadblock to the advance of MHD technology: 
sO much so, that background studies leading to the development of a com- 
mercial MHD generator have perforce been conducted with a large, pulse- 
Operated conventional copper-coil electromagnet (15, 16). 

An alternative to the utilization of superconducting magnets for the 
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above applications would be the use of so-called ''cryogenic'’ magnets: 
that is, magnets utilizing coils wound of any of several highly-purified 
and carefully annealed metals, designed to take advantage of the marked 
resistance decrease of 'ordinary'' conductors at cryogenic temperatures, 
as contrasted with superconductivity. Such magnets share the promise of 
relatively economical generation of intense magnetic fields over large 
volumes, but not the problems cited above. Specifically: 

1. Initial material cost is less, 

2. Mechanical properties tend toward ductility. 

3. No significant restrictions need be placed on conductor size. 

4. The metals contemplated are those considered ''good'' conductors, 
even at room temperature. Thus, resistivity is not critically sensitive 
to minor temperature fluctuations. 

Oo. No serious scaling problems are known to exist. In fact, many of 
the problems associated with cryogenic magnet design become less severe 
as size is increased. 

Of course, these advantages are offset by relative disadvantages: For 
example, since cryogenic coils do have non-zero resistivity, ohmic 
heating occures within the coil dewar, and an appreciably larger refriger- 
ation system may be required for a given cryogenic coil than for a com- 
parable superconducting coil. This implies both higher operating costs 
and higher equipment costs, with the relative increase depending primarily 
on coil size. Again, larger coils would appear more competitive ona 
relative basis. However, two facts seem pertinent, here: First, even 
relatively small cryogenic coils can be designed to operate considerably 
more economically than comparable room-temperature magnets. Second, 
cryogenic magnets have been built and successfully tested in sizes larger 


than those yet attainable by superconducting technology. 
14 





Eee otate-ol -the-art: 

Interestingly enough, although at least one study indicates that large 
cryogenic coils might closely approach the anticipated efficiency of the 
yet elusive superconducting coils for large-scale MHD applications (2), 
none have been built for this purpose, and in fact, no detailed study of 
such an application has been published to date. Perhaps this can be 
attributed to the glamour of superconducting magnets. However, con- 
siderable interest has been exhibited in the development of cryogenic 
magnets in geometries suitable for controlled fusion experiments, as 
can be seen from a brief description of coils already constructed in this 
country: 

Los Alamos Scientific Laboratory (A. E. C. )(17, 18): 

Size: 2-1/2" 1.D. x 7-1/2" O.D. x 5" long. Construction: “Jelly roll”: 
120 turns of 5'' x .010'' copper foil, separated by inert spacers for 
coolant channels, Coolant: Liquid hydrogen (immersion). Field: 60 KG 
steady, 62 KG for 5 second pulse, Status: Operational. 

Project Sherwood (N. B.S.) (19, 20): 

Poze; 3 1.D. x 12’ 0,.D. x7 long. Construction: “Tape wound > 14 
stacked pancakes, each 0. 4" thick, wound of 14,000 turns of 4 mil alu- 
minum foil (insulated by 0.17 mil paper strip) and encapsulated in epoxy. 
Radial slots machined in inert spacers between pancakes for coolant flow. 
Coolant: Liquid hydrogen (forced flow). Field: Tested to 95 KG. Status: 
Structural failure at 95 KG: magnet burst when Lorentz forces exceeded 
inherent strength of epoxy-bonded Al foil. 

Lewis Research Center (N. A. S.A. ) (21, 22): 

Size: 30 cm. I.D. x 90 cm. O.D. x 60 cm. long. Construction: ‘Jelly 
roll'': 12 stacked pancakes, each with 75 turns of 5 cm. x 0.2 cm. alum- 


inum strip, separated for coolant flow. Coolant: Liquid neon (immersion), 
12 





Field: Designed for 200 KG. Status: Undergoing preliminary testing. 

Lawrence Radiation Laboratory (U. C. ) (23, 24): 

Size: 7-1/8" 1.D. x 16" O.D. x 12" long. Construction: 77 turns of 
0.625'' square conductor (sodium cast into . 010" stainless tubing), wound 
in four layers with coolant grooves and stainless support cylinders between 
adjacent layers. Coolant: Gaseous helium (forced flow), Field: Designed 
for 100 KG. Status: Failed during testing at 21 KG when Lorentz force 
extruded sodium through small crack in stainless jacket. Being rebuilt. 

The above coils share many common features: All are cylindrical 
solenoids, such as might be suitable for construction of a magnetic mirror. 
All have what has hitherto been considered a large bore, for high-field 
magnets. All but the last (designed by Taylor and Post) are of extremely 
simple design, and even the last design appears simple contrasted with 
what might be required to perform the same task with "conventional" 
electromagnets. 

Taylor and Post are evolving what might be considered the second 
generation of cryogenic magnets: Their design represents the first 
attempt to obtain the best possible performance (25-27) from a cryogenic 
magnet, as contrasted with previous attempts to obtain performance, 
period. They have attempted to optimize their design in two respects: 
first, they have attempted to identify that operating temperature which 
minimizes total power expended in removing eis losses from the coil, 
and second, they have identified the particular geometry which maximizes 


magnetic energy density per unit dissipation for a cylindrical solenoid. 


C. Application to MHD: 
The coil model which will be used throughout this analysis is the 


overlapping cylinder or "saddle coil'' model discussed by Stekly, et al. 
13 





(2) (Fig. 1). Conceptually, this model consists of two identical, parallel, 
circular cylinders, carrying a uniform, equal axial current in opposite 
directions, which have been brought together until they actually overlap. 
Since, by superposition, no current flows in the region of overlap, the 
material in this region can be removed to provide the "working area’ 

of the coil. Such a configuration confers many advantages for trans- 
verse-field MHD applications: First, referring to the axes depicted in 
Fig. 1, it can be shown that in the working area, the magnetic field pro- 


duced by the overlapping cylinders is uniform everywhere and is given by: 


Bee 0, B= 2 (1) 
where J, is the effective current density in either cylinder and a is the 
offset between the center of either cylinder and the center of the working 
area. Next, insofar as MHD channels tend to be relatively long and 
narrow, the saddle coil configuration offers a relatively efficient utili- 
zation of conductor material, in that the coil windings are so arranged 
that the individual conductors are kept close to the MHD channel which 
they serve. Finally, the configuration is practical from an engineering 
viewpoint: The indicated cylinders would, in practice, be comprised of 
a large number of individual conductors, so that the ''crescent moon’ shape 
of the coil lobes could be readily approximated by the assembly of modular 
units. In short, it seems reasonable to assume that actual MHD coils 
would closely resemble such a model, with the degree of similarity prob- 
ably controlled by a compromise between desired field homogeneity and 
ease of assembly. 

It is apparent that many of the characteristics of a particular coil 
will be determined by the application for which it is intended. For example, 
it has been shown (28, 29) that the length of a transverse-field MHD channel 
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FIGURE 1 


Saddle Coil: Side View 
A 
MHD A End 
Channel Turns 






Current flow into 
paper in this 
lobe 


Current flow out 
of paper in this 


s v 
. - 2 . e 
Zi x ; : . ; : 
\ \ 4 ; a 
l O b e . X ' VA i ¢ 
: ’ . ‘ f V4 m2 a 
‘ aK \ / : 
. e , 
4 5 P - 
\ ‘ / ‘ ; Yo f 4 
ls ‘\ S . ; a 
ry ‘ 3 \ és J 
‘ Fa 
fo . 
‘ 4 
‘ 
‘ 


Section A-A 


195 





should be 5-10 times the mean channel width, to minimize degrading end 
effects. It will be assumed here that the contemplated MHD channel is of 
rectangular cross section, with width 2w and height 2h. The largest such 


rectangle which can be inserted into a given working area has: 


-~3a t+ Vac + 8R 


4 


and 


ie tealiRe S (el eae (3) 


ch 


where R is the radius of either cylinder. The length of the coil associated 
with such a channel will be taken to be: 

ieee 10V/A (4) 
The total cross section (both lobes) of the coil is related to the above 


parameters by: 


A. = aR [TT 27- sin], J cos “a/R) (5) 

One additional concept must be included in the discussion of such 
coils: At either end of the coil, the individual conductors must pass up 
and over (or down and under) the MHD channel. It will be assumed here 
that these ‘end turns'' do not contribute to the useful magnetic field, but 
do contribute to total loss. Obviously, one would desire to minimize the 
mean length of the end turns, within the constraints imposed by the above 
parameters. Considering this, the mean length of an end turn will be 


estimated as one half the circumference of the circumscribing circle: 


that is, 
Li, =7U(R + a). 

It is then useful to define a form factor, F, as follows: 
F=1+ L/L (6) 


Thus, geometrically similar coils can be described by the common value 


of their form factor. Further, total coil volume, etc., can be estimated 
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as: 
Vtot iy ele iy 
The purpose of this work, then, is to examine the application of 

cryogenic technology to coils such as those described above. The basic 
philosophy of this investigation is to seek that information required to 
support a realistic appraisal of the desirability of cryogenic magnets for 
a particular MHD application, to identify those parameters which can be 
used most conveniently to define a coil during the preliminary design 
process, and to illuminate, in detail, the relationships between those 
parameters and the resultant coil. The range of the investigation covers 
coils suitable for generating fields of up to 100 KG, for use with MHD 
channels of cross section between ,.01and10 square meters. It might 

be noted that, aside from the apparent geometric differences, such coils 
differ from the solenoidal coils discussed earlier in that the constraints 


imposed imply approximate coincidence in size between the smallest 


coils considered here, and the largest cryogenic solenoids built to date. 
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7 Il. RESISTIVITY OF METALS 


A. Background: 

It is appropriate to begin this analysis of cryogenic magnets with a 
brief resume of the physical reasons which favor the application of cryo- 
genic technology to magnet design: Essentially, the resistivity of pure 
metals is a sufficiently strong function of temperature that refrigeration 
to the cryogenic range can reduce magnet power consumption by more 
than the additional power required to run the refrigerator. For this 
reason, cryogenic magnets can be made to operate with significantly less 
power consumption than comparable "conventional'' magnets. 

MacDonald (30) provides an excellent summary of the development 
of the underlying theory: Briefly, Mathiessen (31) first recognized that 
the resistance of metals could be separated into two components, one 
dependent upon the degree of purity, and the other upon the temperature. 
He postulated that these components were independent, and therefore 
directly additive. They have come to be generally known as ‘residual’ 
or 'impurity'' resistance and 'intrinsic'' or "ideal'' resistance, respect- 
ively. The former is attributed to scattering of conduction electrons by 
local deformation of the crystal lattice (e. g., by the presence of impurity 
atoms, or dislocations), and the latter, to scattering by thermal vibration 
of the lattice. Thus, the ''ideal'’ resistance goes to zero as the absolute 
temperature goes to zero. 

No significant advance was made in quantitative analysis of the thermal 
behavior of ideal resistance until quantum theory became fairly well 
developed: In 1929, Bloch (32), utilizing a lattice model first postulated 
by Debye (33) in 1911, evolved the first strictly theoretical treatment to 


predict accurately the low-temperature resistive behavior of "ideal" 
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nretals. Bloch's f gates: 
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where ©, the ma) temperature", is a constant, characteristic ofa 


(8) 





particular lattice. (Actually, Bloch's original solution considered only 

the low- and high-temperature limiting cases, but subsequent investigators 
have found that the above integral expression, derived by him, is essen- 
tially valid for all temperatures. For T<<@, the above expression reduces 
to 0~ (T/9)°, while for TZ, O~T/E. ) 

Within a few years, Gruneisen (34) noted that Bloch's theory could 
be modified to provide a more accurate prediction of the resistive behavior 
of real metals by use of an empirical correction factor: He suggested that 
the resistivity derived from eqn. (8) be multiplied by a factor of the form 
(1 +QRT +07), where Ofand O&, were to be determined experimentally 
for the various metals. While admitting that this technique could not be 
defended on strictly theoretical grounds, he advanced heuristic arguments 
defending the form of the indicated correction, and observed that in most 
cases, the resultant ''fit'' to experimental data was accurate to within 
about 10% of the value of the residual resistance. (Later investigators 
(35) have introduced an equivalent, though somewhat less convenient, 
concept of describing ''real'' resistive behavior by use of a variable 
Debye temperature. ) 

Concurrently with the above developments, several investigators 
became interested in the fact that the resistance of metallic conductors 
increases in the presence of a magnetic field. In 1938, Kohler (36) was 
able to show that: 


Le 
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where Eat r(0, T) is the total resistance at temperature T, in the absence 
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of any magnetic field, and Ar = r(H,T) - r(0, T) is the resistance increase 
in the presence of the field. Kohler was unable to obtain a general sol- 
ution for the functional form of Y, though he predicted that in the limit 
of small arguments, 

YiH/r.) a(H/r)?, 
while in the limit of large arguments, 

YH /r J b(H/ 2)”, 
i.e., the magnetoresistance (Ar) saturates. Obviously, magnetoresistance 
is more important at cryogenic temperatures than at ''ordinary'’ temper- 
atures, Since the reduced value of ea leads to a more rapid increase of 
magnetoresistance with field for any given metal. Further, though Kohler 
did not obtain a general solution for Y, mere recognition that the functional 
relationship exists suggests a plotting technique by which the results of 
limited testing may be generalized to predict the magnetoresistance of 
samples of different purity, or at different temperatures. Plots of 
magnetoresistance vs. H/r, have come to be known as "Kohler diagrams. 
(Equivalently, one can plot magnetoresistance vs. the ratio of mean free 
path in the lattice to the radius of the orbit of a free electron in the mag- 
netic field: However, though this latter concept may be more elegant in 
a physical sense, the former technique is generally more convenient for 
the experimenter. ) 

In 1940, Justi(37) recommended the use of a so-called reduced Kohler 
diagram, in which magnetoresistance is plotted vs. H(ro/ De) rather than 
H/r.. It was his hope that normalizing the field-free resistivity to its 
value at the Debye temperature for the particular lattice would eliminate 
differences due to the variation in ‘intrinsic’ scattering between different 
lattice structures. Unfortunately, such normalization does not cause the | 


data for all metals to fall on a single curve, though the curves for the 
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various-metals, when so plotted, do tend to fall in a rather broad band. 
Nevertheless, many experimenters have adopted the practice of plotting 
magnetoresistance data with field-free resistivity normalized in this 
manner. Actually, little is gained, and much is lost, through this tech- 
nique: In general, authorities cannot agree on the exact value of the 
Debye temperature for a particular metal. Further, since Debye temp- 
eratures for most metals fall well out of the range in which detailed 
experimental resistivity figures are available, the normalizing standard 
for a particular metal must often be generated from the theory. Aside 
from the approzimations inherent in the theory, there is the problem of 
which version to use: Bloch's original version, the Gruneisen modification 
thereto, or some more recent derivation. Thus, one is faced with data 
normalized to a rather fuzzily defined standard, and all too often, pub- 
lished data lack any indication of the normalizing standard used. Even 
when the standard can be determined, this practice obscures any direct 
comparison between different metals. About all that can be said to 
recommend the practice is that it reduces the size of the numbers involved. 
Despite continued interest in the theory of the magnetoresistive 
behavior of metals, little significant advance has been made in the last 
several years. To quote from a second article by MacDonald (38): 


'---it appears that the theory is quite unable to give any 


detailed explanation of such bizarre behavior: Wilson says 'It 
does not seem likely that any model simple enough to be tractable 
theoretically would give a magnetoresistance curve of the com- 
plexity of those actually observed. '"' 


l to 


. Limitations of the theory: 

The above relations, taken together, comprise the theoretical moti- 
vation for the use of cryogenic magnets: In theory, careful purification 
and annealing will result in a residual resistance several orders smaller 


than the room temperature resistance of a given metal. Cooling to cryo- 
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genic temperatures can make the ideal resistance arbitrarily small 
compared to the residual resistance. Then, if the magnetoresistance 
saturates at a sufficiently low value--- or, at least, does not rise too 
rapidly with increasing field--- cryogenic magnets can be made to operate 
much more economically than conventional electromagnets. 

Unfortunately, as might be inferred from the above discussion, the 
approximations inherent in the foregoing concepts preclude their direct 
application to magnet design: 

1, The theory does not quantize the relationship between degree of 
impurity, or extent of cold work, and the magnitude of the resultant 
residual resistance. Nor is anything said regarding the feasibility of 
attaining a specified degree of purity in a particular metal. 

2. Experimental data indicate that Mathiessen's rule is an best an 
approximation. Assuming its validity, one can compute the 'ideal"' 
resistance-temperature curve for a given metal by extrapolating an 
experimental curve to 0° K to determine the residual resistance, then 
subtracting that constant from the original curve. In general, curves so 
computed from observation of samples of lesser purity (34, 39) lie well 
above those computed from measurements on highly pure samples (40-43), 
at least for the liquid hydrogen-liquid helium range. Maimoni's date 
(48) indicate a similar discrepancy for aluminum samples of the same 
purity, in which the resistance ratio was varied by controlling the degree 
of anneal. One can only conclude that either the scattering due to lattice 
deformation is energy dependent, or some energy dependent coupling 
exists between intrinsic and residual resistivity. 

In essence, the original assumption of independence between these 
components of resistance implies the assumption that the density of 


‘impurity''- type scattering centers is so low as to have negligible influ- 
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ence on the overall lattice. Recent investigations may indicate the sev- 
erity of this restriction: Several metals (e.g., Ag, Au, Cu, Mg) exhibit 
an anomalous resistance minimum, below which resistance increases as 
temperature is decreased (30, 44-46). MacDonald and Pearson (47), 
working with extremely dilute copper alloys, traced this behavior to the 
presence of specific impurities: They reported that as few as 9-10 p. p.m. 
of such elements as Fe, Bi, Pb, Sn, Ga, In, Ge, and Si sufficed to pro- 
duce the phenomenon, In the absence of these elements, no minimum was 
observed. (A later paper (48) describes iron as the primary culprit, 
suggesting that the latter five elements listed serve primarily to replace 
traces of Fe originally tied up in oxides, thus releasing them into solid 
solution. ) 

3, Attempts to ‘improve the fit'’ between Bloch's law and experimental 
data--- such as that proposed by Gruneisen--- do not alter the asymptotic 
behavior predicted by that law at low temperatures. However, it has been 
found that several metals deviate from the predicted T° dependency 10m 
example, the resistivity of platinum apparently goes to zero as 7? (49, 50), 
and that of iron and nickel, as sa (51). Perhaps the most startling devia- 
tion is that reported by Mendoza and Thomas (52): Working with copper 
and gold specimens which did not exhibit the minima discussed above 
(which occurs at about 10° K for copper), they found "normal" behavior 
down to about 0. 25° K, below which the resistance of both metals rose 
asymptotically toward the ordinate as temperature was further decreased. 

4, Although Kohler predicted that magnetoresistance should saturate 
at sufficiently high fields, only two metals--- aluminum and indium (37, 
93)--- have experimentally displayed such saturation, to date. Further, 


though data from various sources correlate rather well for some metals, ! 
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when plotted on Kohler diagrams, other metals display considerable 
"scatter'. 

Lest the trees obscure the forest, it should be noted that the behavior 
of most metals conforms reasonably well to that predicted by the theory, 
especially above the liquid hydrogen range. Even in the ‘exceptional 
cases", the absolute magnitude of the deviations is small, if only because 
total resistance is small at the temperatures in question. However, from 
a practical viewpoint, a relative error of a factor of two--- or even, 10% 
--~- might be the difference between success or failure for a particular 
cryogenic magnet. Thus, it is unwise to predict quantitative results, a 
priori, from the theory alone. 

Because of these considerations, the decision was made to conduct 
a careful, detailed survey of the literature regarding the experimental 
resistive behavior of various cryogenic conductors, as a starting point 
for this study. The specific metals investigated were aluminum, copper, 
sodium, and indium: The first three have already been used in the con- 
struction of various cryogenic solenoids, and the last has been proposed 
for use in cryogenic coils, in view of its magnetoresistive behavior. The 
results of this investigation are summarized in the remaining sections of 
this chapter. TO FACILITATE DIRECT COMPARISON BETWEEN CRYO- 
GENIC CONDUCTORS, ALL RESISTIVITY FIGURES IN THIS AND SUB- 
SEQUENT CHAPTERS HAVE BEEN NORMALIZED TO THE RESISTIVITY 
OF "STANDARD" COPPER AT 300° K (WHICH IS TAKEN TO BE 1.7715 
MICRO-OHM CENTIMETERS), except that data displayed in Figs. 2-5. 


C, ‘Residual" or 'Impurity" 


resistivity: 
As a rough rule of thumb, most metals can be purified sufficiently to 


achieve a residual resistivity approximately four orders of magnitude 
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below their room temperature resistivity. Whether or not such purifi- 
cation is feasible, in a particular case, depends on the cost and complexity 
of the necessary metallurigical processing, as balanced against the con- 
templated application: 

For example, at liquid nitrogen temperatures, residual resistance is 
generally negligible compared to intrinsic resistance, for metals of ''com- 
mercial" purity. Thus, the conductors of coils intended for use with this 
coolant need not be of highest purity: Far more importance should be 
assigned to the selection of a conductor with a high Debye temperature, 
to obtain (hopefully) lower intrinsic resistance. (Approximate figures 
for the Debye temperatures of the metals covered by tlis study are: 
indium, 100° K; sodium, 202° K; copper, 333° K; aluminum, 395° K. ) 

By contrast, in the liquid neon range (i.e., about 30° K) and below, 
residual resistance becomes increasingly important: Ultimately, the 
lowest residual resistance attainable for a given metal generally determines 
the minimum operating power attainable for a coil wound of that metal. 
Obviously, then, the minimum residual resistance attainable, and the ease 
with which it can be attained, become important factors in the evaluation 
of prospective cryogenic conductors. 

In general, metals such as indium and sodium, with low Debye tem- 
peratures, can be readily purified by vacuum distilation: indium, toa 


residual resistance ratio as low as rp/Tanq - Dev Ne (GA) evecheekeasn. 


to about Iy/r = Fe on? (41,55,56). Or, normalizing these figures 


300 
to the resistance standard cited above, the ratio for indium becomes 
ro/r, = 2x 10%, and for sodium, r)/r, = 8x 10°%, It is worth noting 


that metals with such low Debye temperatures are fully annealed at room ‘' 


temperature: Thus, once residual resistivity has been decreased to a 
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particular level, it will not increase due to ''cold work" in the manufacture 
of a coil. 

Aluminum, despite its high Debye temperature, can be readily pur- 
ified by the zone refining process (57-59). This technique has proven a 
relatively economical source of aluminum with residual resistance ratios 


sa) However, both 


as low as Lo/Tao¢ =1,4x107° (ro/r. = 2,2 x 10° 
Maimoni's data (40) and the Project Sherwood report (20) indicate that the 
resistivity of aluminum is quite sensitive to cold work. If aluminum cryo- 
genic coils are to derive maximum benefit from the use of highly pure 
metal, conductors mustbe carefully annealed--- after assembly, if 
possible--- and structural support provided to ensure that working stress 
in the windings does not exceed the proportional limit. (It might be noted, 
also, that elimination of impurities tends to reduce the yield strength, 
as well as the resistivity, of most metals. ) 

Copper poses a unique problem: although methods have recently 
been perfected to zone refine this metal (60), its high conductivity com- 
plicates the use of induction heaters in the refining process, thus render- 
ing the technique relatively difficult. To date, most of the ultra-pure 
copper used in recent experiments has been prepared by a complex elect- 
rolytic process developed in 1941 (61). Minimum residual resistance 


a -4 7 -4 
300 7 Seca (ener or ro/r. = 4,8 x 10 


(17, 48,62). However, though the resistivity of pure copper is apparently 


ratios have been in the vicinity of Lo/r 


not so sensitive to cold work as that of aluminum, it is acutely sensitive 
to the presence of Fe in solid solution (as has been discussed above): so 
much so, that actual attainment of the minimum cited is a rare phenomenon. 
(As contrasted with the usual correspondence between decreased impurity 

content and lower residual resistivity, MacDonald and Pearson (47) indicate 


that the resistivity of Fe-contaminated copper can be decreased by reacting 
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the iron contaminant with oxygen, thus raising the total impurity content, 
but converting the iron to less critical iron oxide.) It would appear advis- 
able to avoid even the use of steel dies in drawing pure copper into wire, 
and even then, a copper-wound coil might be a questionable choice for 

low temperature applications. 

The above resistance figures do not, by any means, constitute an 
absolute minimum: Gold, et al. (48) report a specimen of purest natural 
copper with a resistance ratio of 3.1 x ca Even more significantly, 
Post and Taylor have recently developed a process (63) which, while 


altering neither the purity nor the degree of anneal, lowers the residual 


resistance ratio for sodium from the level cited above to ro/Tano Sex 
O° aqyaee = 4,3x ay It is at least possible that their process may 


be effective for other metals. It might also be noted that the advancing 
state-of-the-art has roughly reduced minimum residual resistance ratios 
by two orders of magnitude since the 1930's. It would seem reasonable to 
assume that at least some further improvement is probable. Again, any 
further reduction in residual resistance will generally permit further 


reduction of the minimum attainable operating cost for cryogenic coils. 


D. ‘Intrinsic'' or "Ideal" resistance: 

Figures 2-5 present a summary of the available experimental data 
on the intrinsic resistance of indium, sodium, copper, and aluminum. 

In processing this data, occasional interpretation has been necessary: 
Some of the data were originally presented in the form of curves of ''total"’ 
resistance, with no specification as to the residual resistance of the 
Specimen tested. Such curves have been graphically extrapolated to 
0° K, and the intrinsic resistance computed as discussed earlier. A few 


experimenters have presented curves normalized to their estimate of the 
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Debye temperature resistivity of the specimen tested. In such cases, the 
data have been converted in accord with the Bloch-Gruneisen theory to the 
form presented here: Where the original author specified a value for the 
Debye temperature, that value was used in the conversion. Where no 
value was indicated, the temperatures cited in the preceding section were 
used. 

Since the great bulk of resistivity data for indium, sodium, and copper 
were originally normalized to the One resistivity of the specimen tested, 
the experimental data presented in Figures 2-4 are all normalized to that 
standard. Conversely, published data for aluminum has most frequently 
been normalized to the 300° K resistivity, and is so normalized in Fig. 5. 

For comparison, the Bloch-Gruneisen curve for each metal is plotted 
along with the experimental data. For sodium, copper, and aluminum, 
the coefficients 64 and, used in computing the theoretical curve were 
those determined by Gruneisen (34). For indium, the values were deter- 
mined by the author. 

Not all available data have been included in Figs. 2-5: However, all 
recent works known to the author are included, and the data presented are 
considered fairly representative of the total information available. (In 
fact, it would seem that the trends which appear in the data presented are 
established with sufficient clarity to warrant, perhaps, greater confidence 
in the ‘overall picture" than in any single work. ) 

As is apparent in Fig. 2, experimental data for indium conform quite 
well to the Bloch-Gruneisen theory down to about 50° K. Below this temp- 
erature, the data lie above the theory: The greatest divergence occurs in 
the range AIK to 12° K, where some experimental points lie more than 
a factor of 3 above the theoretical prediction. Despite this divergence, the 


experimental curves exhibit approximate parallelism to the theoretical 
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curve: In fact, it would appear--- though available data are scanty--- that 
the experimental curves again converge toward the theoretical curve 
below 4° K. Note that, due to its low Debye temperature, the intrinsic 
resistance of indium exceeds the residual resistance figures cited earlier 
down to about 3° -4° K, Considering this, there would be little gained 
by seeking further decrease in residual resistance for this particular 
metal, Indium undergoes a field-free transition to superconductivity at 
3.41° K. However, since it is a 'soft'’ superconductor, this transition 
would presumably be quenched, in the unlikely event that indium high- 
field coils were designed for operation below that temperature. 

In many respects, the resistive behavior of sodium resembles that 
of indium, as can be seen by comparing Figs. 2 and 3. However, there is 
generally better correspondence between experimental data and the Bloch- 
Gruneisen theory for sodium. Although some divergence is again apparent 
below 50° K, the data of MacDonald and Mendelssohn coincide with the 
theory to 11° K. The greatest divergence is that reported by Woods, at 
about 10° K, and amounts to no more than a factor of 2. Again, the experi- 
mental data appear to reconverge to the Bloch-Gruneisen curve at low 
temperatures. The residual and intrinsic resistivities of pure sodium are 
of roughly equal magnitude around 0 lp K, in keeping with its higher 
Debye temperature. One would expect that at least some advantage could 
be gained from further reduction of the residual resistivity of sodium, 
since Taylor and Post (26) indicate that ''minimum power’ sodium coils 
Should be operated at about Ose 

Fig. 4 indicates remarkably close correlation between experimental 
data and theory for copper, provided that no iron contamination is present: 


in such cases, divergence is generally less than 10%. (The scarcity of data 
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points below 15° K again represents the severity of this restriction!) The 
data of Los and Gerritsen, and of Gerritsen and Linde, represent typical 
behavior in the more normal situation, when iron contamination is present. 
(These curves are labelled ''qualitative'’ behavior, due to the difficulty of 
determining any ''constant'' value for residual resistance, when total 
resistance’ rises with decreasing temperature. ) Were it not for the dif- 
ficulty of obtaining iron-free copper, these data would indicate that copper 
might be an ideal choice for low-temperature coils, since its intrinsic 
resistance decreases so rapidly. Actually, the magnetoresistive behavior 
of copper is the least desirable of the four metals considered here, as 
will be shown in a later section of this chapter. These facts, taken together, 
limit the desirability of copper cryogenic coils to the upper end of the 
100° K temperature range under investigation. 

The resistive behavior of aluminum, as depicted in Fig. 5, is unique 
in many respects: First, strong coupling between residual and intrinsic 
resistivity is apparent. The three curves plotted from Maimoni's data 


4 4 


represent specimens with ry/r = 1.476 x 10, 4.94 100, anceenc2 


300 
os nO” °, respectively, with direct correspondence between increased 
intrinsic resistivity and increased residual resistivity. This, of course, 
is in direct disagreement with Mathiessen's rule, and is the apparent 
source of the noted strain sensitivity of aluminum conductors. Aside from 
the fact that this behavior results in marked deviation from the Bloch- 
Gruneisen curve (as did Fe-contamination in copper), the nature of that 
deviation is also unique: At the upper end of the temperature span investi- 
gated, aluminum conforms even more closely than did copper to the 


behavior predicted by the theory. There is virtually no deviation between 


experimental data and the theoretical curve. However, at some tem- 
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perature between 65° K and 25° K--- depending upon the residual resis- 
tance of the specimen--- the experimental curves shear sharply away 
from the theory, after which they exhibit a ae dependence, while the 
theoretical curve rapidly approaches the usual T° dependence. This 
behavior does not appear to persist for long: Around 10° K, the experi- 
mental curves converge rapidly, and thereafter appear to parallel the 
theoretical curve, though they do not converge to it. (Fenton, et al., 
reported that their data conformed to the relationship 
O= 5.68 x 1079 41.5 x 107) 9r* 9 Chom 
throughout the range 0°-30° kK. Actually, their data points appear to fall 
increasingly under this line at temperatures below 4° K, and over it at 
temperatures above 25° K.) Because of this "sidestep", one faces the 
paradoxical situation that predictions from the theory appear to be pre- 
cisely correct, for purest aluminum, down to about 25° K, but more 
than an order of magnitude too low at temperatures much below DS Ae 

It is interesting to speculate that aluminum, above all metals, would 
profit from the further advance in metallurgical technique required to 
further reduce minimum residual resistance levels: Even with the above 
behavior, aluminum has such a high Debye temperature that its intrinsic 
resistance decreases more rapidly with temperature than was the case 
with the other metals covered here. Since the magnetoresistance of 
aluminum apparently saturates at a low value, residual resistance tends 
to dominate the performance of low-temperature aluminum cryogenic 
coils. Further reduction of residual resistance would apparently decrease 
both components of ''field-free'’ resistance, and (if Ar/r, does indeed sat- 
urate at a constant value) decrease the maximum Ar due to magnetic 
effects, as well. Such a three-fold gain could only result in a direct 


reduction of coil operating cost, again provided that care is taken to ensure 
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that ''cold work'' does not occur during assembly and operation of the coil. 

One additional conclusion can be drawn from the data plotted in Figs. 
2-5: The Bloch-Gruneisen theory is a useful asymptotic limit for the 
experimentally derived curves. At liquid nitrogen temperatures, data 
points tend to fall on the theoretical curve, or so near it, that the theory 
itself is a perfectly adequate design criterion. But to repeat the obser- 
vation made earlier, the theory should not be used as a design standard 
for cryogenic coils intended for service in the liquid neon range, or 


below. Recourse must be made to experimental data for such applications. 


EK. Zero-field resistance standards: 

If the Bloch-Gruneisen theory is not to be used as a design standard 
for cryogenic coils, the time has come to develope a standard. Given 
the above data, this is a relatively simple task, but first, it may be 
profitable to digress to a point of philosophy: 

The primary motivation for use of cryogenic coils for any application 
is the desire to reduce operating costs. Even if a particular cryogenic 
coil is not designed to seek minimum operating cost (for example, geometry 
may be altered from the ''minimum cost'' optimum to reduce coil size, or 
operating temperature altered, because of preference for a particular 
coolant, etc. ), the primary justification for accepting the greater com- 
plexity and initial cost of cryogenic coils, as compared to ‘conventional’ 
coils, is still the same. 

Bearing this in mind, the ''total field-free resistivity’ standards 
evolved here will assume that every effort has been made to achieve 
minimum residual resistivity, in each case. As was pointed out above, 


this probably would not be done for liquid-nitrogen-cooled coils: However, 
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intrinsic resistance tends to dominate the performance of such coils to 
such an extent that the error incurred by assuming too low a value for 
residual resistivity is negligible. If, on the other hand, a coil is designed 
for operation at temperatures where residual resistance is a significant 
part of total field-free resistance, the above assumption is strongly justi- 
fied: At such temperatures, total operating cost is dominated by refrig- 
eration costs. For every unit of heat generated in the coil, many units 

of work must be expended in refrigeration. It is highly probable that the 
additional cost and effort required to obtain and maintain minimum residual 
resistance would be repaid many times over the operating life of the coil. 
The importance of this fact cannot be overemphasized! 

For this reason, the residual components of the resistance standards 
presented in Tables 1-4 are the minimum figures previously cited for each 
metal. The intrinsic component corresponds to the Bloch-Gruneisen 
theory, down to the temperature at which the purest specimens tested 
diverge from that theory, and then follows the experimental trend. All 
figures have been normalized to the 1.7715 micro-ohm centimeter stand- 
ard. The design standards for all four metals are also plotted in Fig. 6, 
for ease of direct comparison. 

The detailed application of these resistance standards to the analysis 
of cryogenic saddle coils is discussed in Chapter VII, and the mean resis- 
tivity design curves resulting from that analysis are presented in Appendix 


iD. 
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TABLE 1 


FIELD FREE RESISTIVITY OF INDIUM 


0 1 2 S 4 ) 6 7 8 9 
eect 1.987 2.00 2.15 3.52 6.86 13.44 25,2 44,4 #63.4 
moe. 146,2 197:8 258 326 404 485 276 6795 T79 
892 CO eee 2d ee 2 so ed 1507 1642 iis Jats Zoe 
2200 2390 20900 2650 2790 2940 3090 £3240 33690553540 
3690 3840 3990 4140 4300 4450 4600 4760 4910 £5060 
9220 5380 5540 5700 5850 6010 6170 6330 6490 £46650 
6800 6960 7110 7280 7440 7600 7750 7810 8060 8220 
8380 8540 8700 8860 9010 9160 9320 9480 9640 £9800 
9950 10110 10260 10420 10570 10730 10880 11040 11190 11350 
P2910 117660 11820 11980 12140 123800 12460 12620 12780 12960 


Values of 10°r/r, are tabulated vs. temperature in degrees Kelvin. 


Top temperature scale: units. Left temperature scale: tens. 
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10, 30 
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904 

1182 
2020 
2940 
3930 
4930 


9910 
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4.33 


13, 23 


140. 0 


209 

1262 
2110 
3040 
4030 
9030 


6010 
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4. 39 


1 a 8, 


165.7 


619 


1343 


2200 


3130 


4130 


9130 


6110 


Jiveu subi ay 


2 
4, 40 


22.3 


195. 9 


680 

1424 
2300 
3230 
4230 
9230 


6210 


4,91 
29.1 
230 

743 

oO 
2390 
3330 
4330 
9330 


6310 


FIELD FREE RESISTIVITY OF SODIUM 


ete 
268 

810 

1590 
2490 
3430 
4430 
9420 


6410 


47.9 
308 

880 

1673 
2980 
35 30 
4530 
9910 


69510 


09, 8 
302 

953 

1758 
2670 
3630 
4630 
9610 


6600 


Values of 10°r/r, are tabulated vs. temperature in degrees Kelvin. 


Top temperature scale: 


units. 
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Left temperature scale: 
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TABLE 4 


FIELD FREE RESISTIVITY OF ALUMINUM 


Pio. 3 
300 
099 
998 
1470 


2003 


8,23 


37.8 


133. 2 


324 
event. 
1044 


Pog 


2060 


9, 39 


43.7 


147.8 


390 
666 
imeheRl 


1569 


7s NAY 


10, 80 


00, 1 


163. 0 


378 


703 


Poo 


1620 


2174 


Wan Sa) 
Dt. 3 


ie 


407 
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69. 3 


epee, 


436 


781 


i252 


1727 


2288 


17, 61 


74.1 
219 
467 
822 
Le Tie, 
io 


2349 


0. 27 


20. 6 


84. 0 
234 
498 
869 
1325 
18395 


2402 


0. 84 
24.1 
94,9 
299 
029 
908 
1374 
1890 


2459 


Values of 10%r/r. are tabulated vs. temperature in degrees Kelvin. 


Top temperature scale: 


units. 
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Left temperature scale: 
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FIGURE 6 


Field-Free Resistance Standards 
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F. Magnetoresistance: 

Basically, the comments in Section D of this chapter regarding the 
processing of experimental data apply as well to this section. In addition, 
a few of the problems associated with the collection of experimental 
magnetoresistance data bear discussion, prior to presentation of that data. 

Assuming that Kohler's rule is valid, the most convenient manner to 
collect magnetoresistance data is to obtain the purest specimen available, 
immerse it in liquid helium, and subject it to various magnetic fields 
while conducting resistance readings: The purity and low temperature 
will assure a low value for ro» and thus a high value of ce for rela- 
tively low values of H. These data can then be projected to much higher 
field intensities, for coils operating above the temperature of liquid 
helium. 

Unfortunately, aside from the obvious problems of instrumentation 
and measurement, and the difficulties inherent to any work in liquid 
helium, such collection of magnetoresistance data is beset by a number 
of unique problems: 

As zero-field resistivity is quite low, it would be convenient to work 
with a relatively long, thin wire to maintain conveniently large resistance 
values. However, size effects are known to influence magnetoresistance 
readings: The resistance of a fine wire actually decreases in the presence 
of a magnetic field, until the field intensity reaches the point that the 
Larmour radius of a free electron in the field is of the order of the radius 
of the wire. Beyond this point, magnetoresistive behavior is more in 
keeping with the theory. Still, the extent to which size effects may con- 
tinue to influence behavior, even beyond this crossover point, is difficult 


to determine. Unfortunately, the working area available to the experi- 
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menter is obviously limited by the size of both the dewar and the magnet 
providing the test field. Due to these problems, very little data have 

been taken from ‘'large’' specimens: that is, specimens approaching the 
size of conductors which might be used in cryogenic coils. In fact, the 
very small size of some specimens tested leads to problems in correlat- 
ing the available data, and questions as to whether such data can be accur- 
ately projected to coil performance. 

Orientation is a second problem: it is known that the magnetoresis- 
tance of a polycrystalline specimen rises more rapidly if the field is 
applied transverse to the specimen, than if the field is applied longitud- 
inally. The data presented here are all transverse-field data, in keeping 
with the anticipated application. Yet this factor also introduces an obvious 
element of uncertainty. 

Finally, the range of the data becomes a problem: Specimens are not 
always so pure and strain-free as might be desired, and attainment of 
sufficiently high values of H/r, can be a problem. The usual solution is 
to use a pulsed magnetic field, since few magnets exist which are capable 
of producing the necessary field intensities on a steady-state basis. Such 
dynamic readings introduce yet another element of uncertainty: The spe- 
cific heat of most metals is quite low at the temperature in question, 
Pulsed fields induce eddy currents in the specimen; eddy currents result 
in Ohmic heating; and the temperature of the specimen--- and therefore 
its resistivity--- may rise from this effect. Some experimenters have 
corrected their data for such factors, or have estimated the effect to be 
negligible. Others do not even discuss the problem. One would anticipate 
that such an effect might lead to artificially high values of Ar, particularly 
at high field intensities, but the degree of correction required to compensate 
for this may be a moot point. 


In any event, it is apparent eet eas collection of magnetoresistance 





data is a somewhat more arduous undertaking than the collection of zero- 
field resistance data, so it is perhaps not surprising that the data dis- 
played in Figs. 7-10 exhibit somewhat more scatter than was apparent in 
Pigs. 2-5! 

Figures 7 and 10 reveal marked similarity in the magnetoresistive 
behavior of indium and aluminum: For both metals, Ar it. apparently 


saturates at Hae = 107 


KG (where no r/r.), and thus never rises much 
above the value of 1.5. Further, though sources of data for indium are 
relatively limited, it is generally true for both metals that experimental 
data from different sources correlate reasonably well. (Purcell's data, 
as published in NBS Report 8481 (20), are an exception to this: His data, 
on aluminum foils, exhibit significantly higher values of Ar/ a than have 
been reported elsewhere. Further, the magnetoresistance of Purcell's 
specimens apparently increased with increasing temperature, fora 
given field strength. Due to this latter phenomenon, his data cannot be 
reduced to a single curve by application of Kohler's law, and are not 
reproduced here. The extent to which size effects may be responsible 
for such behavior is not known. ) 

Of the data reproduced in Fig. 10, the only significant deviation 
from the trend cited above is that reported by Cotti, et al.: Their data 
points indicate a possible departure from saturated behavior at high 
field intensities. Other experimenters have generally failed to observe 
such a phenomenon, and Cotti describes any such departure as ‘spurious’. 
In particular, Callaghan--- whose 3/8 diameter Specimen may be the 
largest aluminum specimen for which magnetoresistance figures are 
available--- tested to virtually the same range, with no apparent deviation 


from saturation. 


Figure 8 presents the experimental data for sodium, and is, in itself, 
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FIGURE 7 


Magnetoresistance of Indium 
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FIGURE 8 
Magnetoresistance of Sodium 
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FIGURE 9 


Magnetoresistance of Copper 
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FIGURE 10 


Magnetoresistance of Aluminum 
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ample evidence of the "bizarre behavior" alluded to earlier. Hopefully, 
the situation may not be quite so chaotic as it would seem, at first. 
White and Woods worked with very fine (less than . 35:mm) rods of 
sodium, and as a result, the magnetoresistance of their specimens was 
found, initially, to decrease with increasing field. The data points 


"crossover 


reproduced here are limited to those which fell above the 
point'' discussed earlier, but again, the extent to which size effects might 
still exert an influence is not known. In any event, the experimental 
curve exhibits a "'knee'' at about H/r, = 300 KG, as though magnetoresis- 
tance had begun to saturate, but then passes through a point of inflection 
and asymptotically approaches a linear dependency at high fields. Mac- 
Donald's data overlap the upper end of this range, and confirm the linear 
dependency at high fields: however, MacDonald's curve is displaced lat- 
erally from that of White and Woods. This may, in part, be attributable 
to the fact that MacDonald indicates his calibration to be ''approximate’’. 
The most recent data reproduced here are from Taylor and Post: 
Their curve results from measurements on a particularly large (sodium 
cast into a 20' length of 0.5" O. D. x 0.035" wall stainless tubing) and 


partucularly pure (r ~ 6500) specimen, wound into a 4-1/2" 


273/T4, 2 
diameter coil. Qualitatively, the behavior of this larger specimen is 
similar to that described above. The curve exhibits a ''knee'', though a 
much more gradual one, and apparently passes through a point of in- 
flection near the upper end of the range investigated. Whether or not the 
quantitative differences are attributable to size effects is a matter for 
speculation. Regrettably, current data do not extend to sufficiently high 
values of H/r to confirm whether or not this last curve would eventually 


approach a linear field dependency. Taylor indicates that additional test- 


ing will extend the data to higher fields. 
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Finally, Fig. 9 summarizes the experimental data for copper. Con- 
trasted with the other metals described here, the magnetoresistance 
curve for copper seems almost simple: It initially exhibits an H? depend- 
ence, then falls gradually to a linear dependence. This fact has been con- 
firmed by numerous investigators, and all reports correlate remarkably 
well, except for that of Olsen and Rinderer., Although the magnetoresis- 
tance of copper is relatively low at low fields, it rises so rapidly that 
by ET = 107 KG, copper exhibits greater magnetoresistance than all 
other metals considered. As ultra-low-temperature, ultra-high-purity 
high-field copper coils would exceed this value of H/r, by an appreciable 
margin, one faces the problem discussed above: High-purity copper 
coils are unattractive in the liquid hydrogen range because of magnetoresis- 
tance. Coils of lesser purity are unattractive in that range due to the 
overwhelming probability of iron contamination. Thus, it is advisable 
to restrict the use of copper coils to the liquid nitrogen range and above. 

The design standards derived from Figs. 7-10 are presented in 
Table 5, and plotted in Fig. 11. In order to permit analysis of high- 
field coil performance in the liquid helium range, it has been necessary 
to extrapolate the data to H/r = 10° KG. Indium and aluminum have been 
extrapolated at constant Ar; copper has been extrapolated under the 
assumption that Ar varies linearly with field intensity; and sodium has 
been extrapolated under the assumption that magnetoresistance asymptoti- 
cally approaches a linear field dependency. Extrapolated values are indi- 


cated by parentheses in the table, and broken lines in the curves. 
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TABLE 95 


MAGNETORESISTANCE DESIGN STANDARDS 


By Tr (KG) Imerys r Na: Ar/r Cu: Ar/r Al: Arye 
me 10" = = - 
2 00150 = a x 
3 00270 as iat 00110 
4 00409 a as 00197 
5 00565 - _ 00309 
5 00730 oe ee 00448 
7 00855 - 00140 00606 
8 00949 es 00182 00794 
9 01067 aa 00229 01008 
1 x 10° 01228 00100 00284 01234 
2 0418 00390 01093 0448 
3 0811 00830 0230 0840 
4 1232 01352 0378 1268 
5 1660 01915 S545 1701 
6 212 0270 0716 214 
7 258 0316 0894 256 
3 304 0382 "1080 298 
9 349 0447 1273 _ 339 
10° 393 0513 1470 563 
2 162 1101 364 676 
3 020 1613 599 872 
4 205 206 838 1. 008 
2 340 243 1.111 1.112 
6G 441 276 1. 342 1.187 
7 523 306 1.591 1. 249 
3 589 334 1.844 1. 295 
g 640 358 2.11 1. 338 
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Magnetoresistance Standards 
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III. STRUCTURAL DESIGN CONSIDERATIONS 


A. Basic concepts: 

The non-conducting structural components of any air-core electro- 
magnet must be regarded as a necessary evil: While necessary to support 
the coil against the various stresses developed, such components con- 
tribute nothing to the desired objective of generating a specified field 
throughout a specified region. Suffice it to say that if coils could be made 
to operate without such elements, they would be. Since this is generally 
impossible, this chapter will be devoted to the development of a model 
with which the minimum structural requirements of a cryogenic saddle 
coil may be estimated, as a function of the basic variables defining the 
coll. 

Note that the above does not imply that the model used provides 
answers which can be directly applied to ''real'’ coils. For one thing, 
the minimum structural fraction achievable in real coils is a function of 
such things as the ingenuity of the designer, the need to provide for ease 
of construction, etc. For another, a model which determines the min- 
imum structure ''needed" for support may quite reasonably predict that 
moO structure is necessary, ina particular instance. This, for example, 
would be true for small, "jelly-roll'' cryogenic solenoids, in which the 
configuration permits the conductors to be essentially self-supporting. 
However, even in small '‘jelly-roll" coils, it is convenient to use non- 
conducting elements to ensure centering, provide for coolant flow, etc. 
Thus, the structural volume fraction of a real" coil may not be large, 
but will not be zero. 

Nevertheless, the model proposed here provides answers which are 


useful in the preliminary design phase, which is consistent with the spirit 
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of this work. No attempt will be made to evaluate those variables which 
are primarily a function of the discretion of the designer. Corrections 
for such items can best be injected during the detailed design phase, after 
preliminary analysis has broadly defined the form of the desired coil. 

In proceeding with this analysis, it will be convenient to discriminate 
between two classes of structural elements: 

Internal structure is here defined as those structural components 
contained entirely within the coil dewar (and therefore, within the coil 
windings). The primary influence of such structure on the design is to 
reduce the packing factor attainable, and thus increase the conductor 
current required to obtain a given effective current density. (Note that 
equations (1), (3), and (5) essentially define the form of a coil in terms 
of effective current density. ) 

External structure is defined as those structural components which 
are external to the coil windings, including the dewar and those structural 
elements which penetrate the dewar. External structure has no influence 
on the effective current density, but does decrease the available working 


area, add to coil weight, and contribute to heat leak. 


B. Internal structure: 
There are two significant differences between the internal structure 
of a saddle coil and that of a solenoid, such as those described earlier: 
First, the conductors of a tape-wound or "jelly roll’ solenoid have 
structural continuity in the plane of main stress, while those of a saddle 
coil do not. The structural strength of the conductors of a tape-wound 


coil can directly support part or all of the stress load. By contrast, the 


conductors of a saddle coil run axially, while the major stresses occur 
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in the transverse plane. The latter situation is clearly less favorable. 

Second, the magnetic stresses developed in a high field solenoid tend 
to pull the coil apart. Thus, it is relatively easy to design internal 
structural elements so that all structural components are in tension. ‘The 
normalized magnetic stress distribution for various saddle coils is shown 
in Figs. 12-16. (These figures also indicate coil profile, as a function of 
the form factor defined in eqn. (6).) As can readily be seen, it is impos- 
sible to support such a stress distribution by components entirely in ten- 
sion. The point, here, is that structural components bearing compressive 
loads may be subject to elastic instabilities--- i.e., buckling--- at stress 
levels well under the yield stress of the material used. Allowance must 
be made for this in coil design. 

Two other considerations apply to the internal structure of cryogenic 
coils, in general. Mention has been made of the necessity to minimize 
residual resistance for ultra-low-temperature coils. This might neces- 
sitate designing to allow for a final anneal after the coil has been assembled. 
(Or, if the coil is to be built up from modular units--- that is, pre-formed 
“bundles” of conductors--- designing for anneal of such units, after shap- 
ing.) Such goals probably preclude the use of any structural materials 
which depend on elaborate heat treatment to develope structural strength. 
However, the designer has one consolation if ''as annealed'' metals must 
be used for internal structure. Engineering metals which do not become 
brittle at cryogenic temperatures exhibit an appreciable increase in 
strength, over room-temperature characteristics. (See McClintock and 
Hauser (89) for a discussion of factors influencing mechanical properties. ) 
Since internal structure will be maintained entirely at cryogenic temper- 
atures while the coil is operating--- and magnetic stresses do not exist 


while the coil is not operating--- advantage can be taken of the increased 
ov 
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low-temperature strength of structural metals. (For example, Swenson 
(90) found that the compressive yield strength of Type 304 stainless steel 
increased from about 75, 000 p.s.i. at room temperature to over 110, 000 


Pas.i. at 20° K. For additional information, see also (91). ) 
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C. Superposition model for Mitetiia lest euectune. 

As is apparent in Figs. 12-16, the stress distribution in a high-field 
saddle coil is rather complex. It is probable that detailed structural 
design will be a difficult, iterative analysis, with the final outcome 
strongly dependent upon the arrangement scheme initially assumed. In 
order to reduce this situation to one which is sufficiently tractable for 
preliminary analysis, a superposition model will be used. Basically, the 
illustrated stress pattern results from two separate force fields: that due 
to the interaction between the current and magnetic field of a single coil 
lobe, and that due to the interaction of the current of one lobe with the 
field of the other. Each of these force fields will be treated separately, 
and the optimum structural arrangement to support each field will be 
determined. The structural volume fraction required for the combined 
fields will be estimated as the sum of the structural volume fractions so 
determined. 

As will be seen, this is equivalent to assuming that the support 
structure of the coil is fabricated from a series of interlocking segments 
of circular shells. This probably would not be feasible, for a ''real" coil, 
due to the obvious problems of assembly. However, the circular shell 
model does indicate those regions of the coil which most need structural 
support, and provides a tentative arrangement which could be successively 
modified in an orderly fashion. 

Specifically, application of the superposition technique to this problem 
involves the following approximations: 

1, Interactions between the two force fields are ignored. (If one were 
to draw an 'X"' through the center of either cylinder, the two force fields 
would be found, generally, to oppose each other on that side of the ''X" away 


from the channel. Thus, less structure would actually be required in that 
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region than the model would predict. ) 
2. Interactions between support shells are ignored. (The shells 
required to support ''self field'' are subjected to compressive loading, 
and must be designed in accord with the appropriate elastic stability 
criteria. Actually, the interlocking shells resulting from the ‘'adjacent 
field'' analysis would provide some additional stiffening, not considered 
in the ''self field'' analysis. On the other hand, if either set of shells 
were to shift position slightly under load, part of that load would be trans- 
ferred to the other set. Such load transferral is also neglected. ) 
Basically, predictions from such a model should be reasonable, but 
conservative. At least, the model provides a systematic means to esti- 
mate the additional structure required due to the inability to avoid com- 
pressive loading. Further, it is well to be conservative at this point. 
Deviations from the idealized circular shell model--- in the interest of 
ease of assembly--- will probably necessitate a somewhat higher struc- 
tural volume fraction than would result from an ‘optimum support’ arrange- 


ment. 


D. Support of ''self field": 
The force distribution in a circular cylinder, due to interaction of 
uniform, axial current flow with the magnetic field generated by that 


current flow, is directed radially inward and has a magnitude given by: 
rJ° 
2 


To return to our original hypothesis, where the use of structure can 


f= magnetic force density = (10) 


m 
be avoided, it will be. If the cylinder is comprised of many individual 
conductors, those near the center of the cylinder are in such a region. 
Force densities are low, and directed so as to hold the coil together. It 
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is reasonable to assume that the wires in any circular shell would tend to 
"jam together’, as do the keyed stones of a circular arch. Certainly, the 
conductors could be shaped to accentuate this tendency. Stability would 
not be a problem, Since this region is essentially solidly filled with con- 
ductors. To what radius does this region extend? 

Consider a thin shell of conductors at radius r. If we assume this 
shell to be ''keyed together’, it is essentially self-supporting. The mag- 
netic force distributed within this shell is equivalent to an external pres- 
sure, directed inward, of magnitude 


LL s*r dr 


p 
a y) 


Neglecting stability, this is simply the classic ‘hoop stress'' problem. 


If we desire to limit the internal stress in the conductors to a maximum 


level O.: we have: 5 : ae ; [So r“dr 
c 
t 7, (ale 


and therefore 


Ih JAG) 


"max ks : (11) 
JV KM 


That is, the individual conductors must not be ‘keyed together" at 





fey OF the resultant hoop stress will exceed the specified stress 
limit of the conductor. If the coil must be designed to prevent self support 
in shells at greater radii, the conductors in such shells will be forced in- 
ward by the magnetic pressure, due to lack of structural continuity. The 
resultant stress must not be transferred to the shell at r = es es since 
the conductors of that shell have already been ''stressed to the limit". 
Therefore, the first structural support shell must be located at r, > ae 
Assuming the conductors located outside the first support shell all 


transfer their load inward, those located immediately adjacent to the sup- 
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port shell will be subjected to the greatest load, in the form of simple 
radial compression. Again, this must not exceed the specified limit, 
0. The position of the second support shell can then be determined by 
integrating eqn. (10) radially, subject to this condition. 


€ eee In a slab of width AO, 


2 2 
UC ea) = Le lie KS), J 








—_— The "equivalent pressure’ on the structural shell 
is thus: uae 2 3 i rj°-r,° 
p. = AF/AA = (r,°-r, )A0O/r, AS = — | ——— (12) 
© F 2 1 1 F z 


And the radius of the second structural shell must be: 
60 WES 
z Cc 
r, = r 1 


+ ee 
1 ULAR, = 


Obviously, this can be repeated until the outer radius of the cylinder 


(13) 


is reached. Eqn. (13) is used to determine the radius of each successive 
shell, since it is known that Da 0, in eqn. (12). For the last shell, eqn. 
(12) is solved for the effective pressure on the shell, by setting r, = r_, 
the outer radius of the cylinder. 

The above analysis of "effective pressure’ has ignored the thickness 
of the support shells. The reason is this: Consider a slab of height h, 
-_— | mounted on a zero-thickness shell at y = 0, and carrying 
. h a uniform current a in the -z direction, In the presenee 
of a uniform, x-directed magnetic field, BL. the slab 
will exert a pressure Pa Bld ih on the support. Now let the support 
have non-zero thickness, so that it extends from y= 0toy=t. If we 
require the same total current in the slab, the current density must be- 
come J = J/X » where A = (h-t)h. The pressure exerted on the support 
is now given by p, = BAI /A ih A) = BJh, as before. An analogous 


situation exists in the "self field'' problem above, since the magnetic field 
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at a given radius in such a composite cylindrical conductor is primarily 
determined by the "effective'’ current flow within that radius. Thus, toa 
first approximation, the location of the "next" shell is not influenced by 
the non-zero thickness of the ''present'' shell, and the appropriate current 
density to use in eqns. 10-13 is Jas the effective current density of the 
coil. Obviously, field perturbations due to the presence of non-conducting 
structure become more severe, and therefore the accuracy of this model 
decreases, as A>0. (In fact, when the stability criteria are applied to 
determine Shell thickness, it is quite possible to compute a "structural 
volume fraction" greater than unity.) However, if the analysis indicates 
that support structure alone would result in a packing factor less than, 
say, 0.25, the coil is obviously not of interest. 

Once the position of each shell has been determined, its required 
thickness must be found. Windenburg (92) summarized, in 1947, the 


equations required for this process. They are: 


2 1 
p=— Or (14) 
wae " 1.05 D vr 
ao 2. 42E (t/D) 
p= 75 a/a a ea le) 
uy Ey ABs 045(t/pi!? 
1 2E 3 
elm 1) Say JE), (16) 
LoL. Ie EN el 
Cine ' 
p= 2.340 — -1. 06( O JE)! 2 (17) 
) Dp yy 


Equations (14) and (15) apply to 'short'' cylinders, while (16) and (17) 
apply to ''long'' cylinders. Essentially, if a cylinder is very short, elastic 
instability is no problem, and (14)--- which is a slightly more conservative 
version of the hoop stress equation--- applies. As the cylinder is length- 
ened, a simple, two-lobe instability collapse is possible at stress levels 


below yield. The ‘critical pressure’ in this region is a function of the 
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elastic modulus, t/D, and L/D, as shownin (15). At still greater lengths, 
the ends are so remote that further removal has no significant influence 
on stability. At this point, (16) applies. However, as (16) was derived 
from elastic stability considerations, it does not recognize that a pro- 
portional limit exists. Thus, collapse can occur from yielding, vice 
instability, at pressures lower than predicted from (16). The appropriate 
equation to predict this mode of failure is (17). 

In using these equations to determine Shell thickness, one must 
obviously know p, t/D, L/D, the design stress limit specified for the 
structural material (which would presumably be an order of magnitude 
higher than that established for the conductor), and the elastic modulus 
and Poisson's ratio for the structural material. For ''short'' cylinders, 
one would solve both (14) and (15). The larger (t/D) so obtained would be 
correct. For 'long'' cylinders, (16) and (17) would be applied in similar 
fashion. To aid in distinguishing ''short'' from ''long'’, Windenburg defines 
the parameter: Zo 10° O L 

ee ee ee a 
K 27,500] \D 

itex— 5, (14) and (15) apply. If X 220, (16) and (17) apply, lio 
X <20, try both (14) and (15), and (16) and (17), as above, then keep the 
smaller of the two resultant values for (t/D). 

One problem arises, regarding the crossover from (16) to (17) for 
(A) high-yield-strength metals. Fora given 
metal (and a particular L and D), the simal- 
taneous solution of these two equations in - 
(t/D) volves a cubic in(t/D). Such a solution can 

assume one of three forms, illustrated in 


liege Uy Fig. 17. In case (A), there are three real, 





unequal roots. In (B), there are still three 
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real roots, but two are identical. In (C), there is only one real root, and 
two imaginary roots (i.e., no physical solution). Equations (16) and (17) 
were originally derived for mild steel. Equation (17) is an attempt to 
express case (B) in terms of the physical parameters of any metal, and 

is successful for mild steels. The physical solution in this case follows 
(16) to the point of tangency, then (17) at higher pressures. Unfortunately, 
for yield strengths in excess of about 50,000 p.s.i., substitution of the 
appropriate parameters into (16) and (17) results in case (C): That is, 

(17) no longer provides the equation of a straight line, tangent to (16) at 
the desired crossover point. 

To circumvent this situation, it is necessary to re-examine the phy- 
sical reasons underlying the necessity for ''crossover’'. Eqn. (16) becomes 
invalid at high pressures, for a particular goemetry, because real" 
metals reach their proportional limit and thereafter deviate from Hooke's 
Law. Thus, (16) predicts too low a value for (t/D), since it assumes a 
constant elastic modulus. Eqn. (17) attempts to correct for this by 
assuming that shell thickness must be increased in direct proportion to 
pressure, beyond a specified point. This specified crossover point would 
presumably lie somewhere between the proportional limit and the yield 
point (2% offset), depending on the characteristics of the metal used. 

A conservative approach will be used here. Following the solution 
of (16), the value of (t/D) so obtained, along with the known value of ae 
will be inserted in the hoop stress equation to determine whether or not 
the resultant stress in the shell exceeds the proportional limit of the 
metal in question. If so, the effective pressure corresponding to the 
proportional limit for the given shell will be determined, and the slope 


of eqn. (16) will then be found at that pressure. Assuming this to be the 
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desired crossover point, these data will be used to generate the equation 
of a straight line, tangent to (16) at the pressure corresponding to the 
proportional limit. The projection of this straight line will then be used 
to determine a corrected value for (t/D) at the effective pressure. 
Finally, after the radius and thickness of each shell have been deter- 
mined, the total volume of these ''compressive’ shells can be determined, 
and divided by the volume of the composite cylindrical conductor. It will 
be assumed herethat the compressive volume fraction of the ''crescent 
moon’ of a saddle coil lobe is approximately equal to the compressive 
volume fraction of a cylinder of the same radius, carrying the same 


effective current. 


E. Support of "adjacent field’: 

The force distribution in either coil lobe, due to interaction of uni- 
form, axial current flow with the magnetic field generated by a similar, 
but oppositely directed, current flow in the adjacent lobe, is directed 
radially outward from the center of the adjacent lobe and has a magnitude 


given by: 22 


— Mo®o 2 


PANG ; 


f 
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where Ee is the outer radius of the adjacent lobe, A natural structural 
arrangement to support such a force distribution consists of a series of 
circular arcs, concentric with the adjacent lobe. The problem here is 
Similar to that preceding. The position, thickness, and in this case, 
length of each such arc must be determined. Since this set of structural 
shells is in tension, no stability problems exist. 

To determine the location of the first tensile shell, let us refer toa 


coordinate system centered on the ''adjacent’ lobe, henceforth called the 


fol 





left lobe. Again considering the lack of 
structural continuity in the transverse plane, 
since the magnetic force considered here 


tends to "pull apart" the individual conduct- 





ors, it will be assumed that the load will be 
transferred radially outward from the center 
of the coordinate system. In the most severely stressed conductors, 
immediately adjacent to the structural shell, we have: 


Wrp2I2 |e Glas 


ar, 


oe = dF/dA - [ares ar/2r,¢ - (18) 
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where 20 is the angle subtended by the circle bounding the right lobe at 


a given radius. Referred tothe indicated coordinate system, this circle 


is: (x - 2a)* + y" = ro’, 


where 2a is the distance between lobe centers. Further, the equation 
of a circular arc, centered on the origin, of radius r, is: 


et yee x? 


Thus, the x-coordinate of the intercept of the circle bounding the 


right lobe, and above circular arc, is: 


r?+4a°-r 


eee ols (19) 
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and therefore 


1 oe ce 


QO = oe elie) = cos” 2 (20) 


4ar 
The above is a bit messy for explicit solution, but can be programmed 
readily for numerical solution. Conceptually, the procedure is the same 
as before. Shells are positioned in accord with the requirement that Dae 
0, at each shell. Shell thickness is determined from the requirement 


that: 79 








t. = gee (21) 


where d is the design stress limit for the structural material. The 
"tensile'' volume fraction then becomes: 
2 
V2 O,/O yr; a 
T.V.F. = —-—————————-,, (22) 
A. 


where A. is the cross sectional area of either lobe. 

Essentially, the above shells transfer the magnetic stresses dis- 
tributed within the windings to the outer edge of the coil: i.e., the ends 
of such shells must be "anchored" by structure bounding the outer sur- 
face of the lobe. Further, as there is a net repulsive force tending to 
separate the two lobes, provision must be made to tie the bounding 
structure of the two lobes together. However, such structure does not 
penetrate the windings, and therefore has no effect on the packing factor 
of the coil. It is therefore treated as ‘'external structure’, as defined 
above. 

It is interesting to note the role of 0, and 0. in determining the 
structural volume fraction predicted by this model. Basically, the num- 
ber of structural shells--- both tensile and compressive--- is determined 
by 0. and the thickness of each shell by O;. 

The details of the actual computational procedure followed in applying 
this model to the coils investigated are provided in Chapter VII, and the 


results of those calculations, in Appendix C. 


Fk, External structure: 
The external structure associated with a ''real'’ MHD coil is a suf- 
ficiently strong function of the intended application as to preclude quanti- 


tative discussion in a study of this nature. Nevertheless, a few points 
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bear mentioning, qualitatively: 

The phrase ‘external structure’ has been defined to include the 
structural elements bounding each coil lobe, as wellas those members . 
connecting the two lobes, the dewar required to maintain the cryogenic 
environment, and the structural members which support the assembled 
coil, both against gravity, and against any drag or thrust generated by 
the MHD channel. 

Most probably, the structure bounding and connecting the coil lobes 
would form the inner vacuum wall of the dewar, Again, advantage could 
be taken of the increased low-temperature strength of materials. 

For coils designed for operation below liquid nitrogen temperatures, 
it is likely that the dewar itself would be of ''two-stage"’ construction. 
That is, everything within the ''cold" side of the dewar would be main- 
tained at the coil operating temperature, but a ''middle layer" of the de- 
war would interpose a ''sandwich" of liquid nitrogen between the cold side 
and ambient temperature. This would reduce heat leak by conduction 
and radiation through the dewar, but far more importantly, allow greatly 
reduced heat leak by conduction along structural elements and electrical 
leads piercing the dewar, into the ''cold box'’, Most of the heat leak 
along such elements could be removed, with reduced refrigeration cost, 
from the liquid nitrogen layer. 

Most probably, the vacuum layers of the dewar would contain either 
a multi-layer ''superinsulation" (93-97), or one of the somewhat less 
effective but more readily applied opacified powders (96-98), to minimize 
heat leak, Such construction would ensure that heat leak through the dewar 
would be negligible compared to other losses. (A reasonable estimate, 


given such construction, would be .05 watts per square meter. ) 
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- From a strictly structural viewpoint, arrangements would presumably 
be made to ensure that heat leak along electrical leads was also negligible. 
McFee (99) indicates a "passive'' method of minimizing such heat leak by 
control of lead geometry. ‘Active’ methods light involve using boil-off 
vapor from the cold box to refrigerate the leads. However, for steady- 
state, closed-loop operation, it would be necessary to provide for recon- 
densation of such vapor, and its return to the cold box. 

Probably, the most efficient geometry for support structure which 
pierces the dewar would be the so-called ''bicycle-spoke'' arrangement, 
in which the coil is supported by a network of long, thin spokes. Careful 
design of such supports will ensure maximum path length from ambient 
to cryogenic temperatures, and therefore minimum heat leak, along 
each support member. Obviously, this is not the most convenient arrange- 
ment from the viewpoint of ease of construction. In small cryogenic 
coils, foR losses dominate so strongly that one might prefer to utilize 
fewer, and more easily installed supports, even at the cost of added heat 
leak. In larger coils, structural leak may be more important. Stekly 
and Woodson (100) provide a means by which the magnitude of such heat 
leak may be readily estimated: Let F be the total force to be supported 
in a given direction (weight, thrust, etc.). Then, 

ee F/O i 
where Gh is the design stress in the support material (note that since 
Supports penetrate the dewar, room temperature figures must be used), 
and A, is the required effective cross section of support structure. Since 


the rate of heat leak along such structure is given by: 


ae 

gq = LS = AT, 
L 
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where k is the mean thermal conductivity (see (101) or (102) for typical 
values), L. the mean length, and AT the temperature drop across the 


structural elements, we have: 


k F 
eee rar (23) 


OO} \L 


A first estimate of coil weight is simply F = m gV, where m is the 
mass density of the conductor, and V the total coil volume. Similarly, 
if the structure supporting the MHD channel must penetrate the dewar, 
the associated heat leak can be estimated by setting F ~(p. - p)AG, 


For a matched load, ie = Fv/2, where v is the mean flow velocity in 


n 
the channel, and an equivalent estimate is F = 2P en! In any event, 
such relations can be used to estimate the importance of support heat 
leak, and therefore, the importance of optimizing support geometry. 
Since the MHD channel operates at or considerably above room temp- 
erature, the ''working area'' of a saddle coil is not fully available for 
the intended application. Some of it must be filled by the type of dewar 
and support structure described above. Again, ‘how much’ is a strong 
function of the particular application intended, and is therefore not 
included in any of the curves generated in this study. To compensate for 
this factor, one can enter the coil geometry curves of Appendix A with 


a somewhat greater A. than is actually required for the anticipated appli- 


h 


cation. 
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IV. REFRIGERATION 


A. Efficiency: 

As was mentioned earlier, the cost of refrigeration is the dominant 
operating cost of ultra-low-temperature cryogenic coils. Such being the 
case, the time has come to discuss refrigerator efficiencies, before pro- 
ceeding further. 

Brief reference to the Carnot cycle will depict the magnitude of the 
problem. In order to maintain a liquid nitrogen environment--- that is, 
to maintain a refrigerated area at 77. ACK: given an ambient temperature 
200 K--- an ideal, Carnot refrigerator would require 2.9 watts of 
input power for every watt of heat removed from the refrigerated area. 
To maintain a liquid helium environment (i.e., 4. 2° K) under the same 
conditions would require 70. 4 watts of input power for every watt of 
refrigeration provided. As is hardly surprising, considering the working 
environment, ‘'real'’ cryogenic refrigerators require appreciably more 
power than these ‘‘ideal'' figures. 

Actually, the spread between ''real"’ and ‘ideal'' performance has 
been reduced significantly, within the last decade. Scott (62) gives a 
description of the various components used in cryogenic refrigerators, 
as well as a discussion of many of the factors limiting refrigerator per - 
formance. Suffice it to say that though the design of a refrigerator may 
be simple in theory, the design of a cryogenic refrigerator is extremely 
complex, in practice. Thus, improved performance has generally been 
attributable to improved technique--- i.e., ‘practical engineering’ '--- 
rather than any ‘‘breakthrough'' in basic concepts. Given the law of 


supply and demand, the effort expended in this field has increased in pro- 


portion to the growth of industrial application of cryogenic technology. 
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Hopefully, further improvement in the performance of cryogenic refrig- 
erators should be expected, under these circumstances. 

Fortunately, it is unnecessary to become intimately familiar with 
the design of cryogenic refrigerators in order to appreciate the two factors 
which most directly influence their relative efficiency (that is, their per- 
formance as compared to the Carnot cycle efficiency). 

First, as with many other mechanical devices, small cryogenic 
refrigerators (~1 watt) tend to be considerably less efficient than com- 
parable large refrigerators (over 100 watts), Of course, the degree of 
improvement with increased size eventually ''tapers off''. Bernert and 
Stekly (103) present curves of relative efficiency vs. size which indicate 
that relative efficiency "saturates' at a power input of approximately 1000 
KW, more or less independently of "cold side'' temperature. The signifi- 
ance of this point is that the refrigeration requirements of the coils being 
analyzed here are sufficiently great that all require ''large' refrigerators, 
by such a standard. 

Second, relative efficiency is a function of refrigeration temperature. 
The engineering problems which must be solved to maintain continuous 
operation at 4. 2° K are obviously more stringent than those confronted in 
designing for operation at 77. 4° K. Thus, not surprisingly, the relative 
efficiencies which have been attained in helium refrigerators are lower 
than those which have been attained in nitrogen refrigerators. To be 
specific, Fig. 18 has been prepared as a working estimate of the input 
power required, as a function of temperature, for 'large’, closed loop, 
continuous-operation cryogenic refrigerators. The "actual’ curve is 
based on an estimate of the present state-of-the-art, with no allowance 
for anticipated future improvements. The Carnot cycle theoretical limit, 


as well as a few specific data points, are plotted for comparison. 
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Refrigeration Power vs. Temperature 
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B. Other factors: 

The steady-state refrigeration requirements of the saddle coils 
treated here generally exceed the refrigeration requirements of known 
cryogenic applications. Again, given the law of supply and demand, refrig- 
erators of sufficient capacity to provide for these needs do not fall into 
the category of ''on-the-shelf hardware''. While it is felt that construction 
of such refrigeration systems is within the capability of existing tech- 
nology, manufacturers are BRdenerendauly discrete, regarding the speci- 
fications for systems which exist only ''on paper". Further, extrapolation 
of size, weight, and cost figures from the region of current practice 
would appear to involve a greater degree of uncertainty than was the case 
with the efficiency figures quoted above. Regrettably, little concrete 
information regarding such factors is available in the open literature. 

In any event, Bernert and Stekly (103) indicate that total refrigeration 


system mass might be estimated from the relationship: 


0,24 
I 


(where M. is system mass in Kg and Wy refrigerator input power, in KW), 


M_ = 15. 7W 
S 


for systems with an input power requirement in excess of 1O KW. Ina 
separate article (106), Bernert indicates that capital costs for 4. Dra 
systems may be estimated from: 


0,92 «< 


33 W =c_ 5 20w 0. 66 


R R 
(where) C. is system cost in thousands of dollars, and Wr is refrigeration 
capacity at 4. 2S K, in watts), for systems with a cooling capacity in excess 
of about 50 watts. Presumably, capital cost for a given cooling capacity 
should decrease as a function of temperature, probably even more rapidly 
than does system mass, considering the reduced complexity of the design 


problem. Certainly, one could scale the above cost estimate in proportion 
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to total system mass, to obtain a crude ‘first estimate’ of the capital 
cost of systems operating above 4. DOK. However, for specific information 
regarding this factor, as well as the other factors mentioned in this sec- 


tion, one would be advised to turn to the various manufacturers! repre- 


sentatives. 
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V. HEAT TRANSFER 
A. Basic concepts: 

For various reasons, it is convenient to divide the treatment of heat 
removal into two separate areas: Heat transfer, discussed in this chapter, 
deals with the flux of energy from the windings, into the coolant, across 
a fluid-solid interface. (Here, two factors are of particular importance 
to the analysis of cryogenic coils: the amount of heat transfer surface 
required, and the relationship between rate of heat generation, and temp- 
erature rise between bulk coolant, and conductor.) Coolant flow, discussed 
in the next chapter, presupposes a given heat flux into the fluid stream, 
and explores the influence of this flux on the bulk flow properties of the 
coolant. (Items of interest will be the bulk temperature profile of the 
coolant, the relationship between inlet and outlet pressures, and packing 
factor, etc. ) 

It is also convenient to categorize the overall heat removal process 
as ''single-phase", or ''two-phase". The former term implies that the 
coolant is a liquid or a vapor, and that energy is removed by sensibly 
heating that liquid or vapor. The latter term implies that the coolant is 
a Saturated liquid, and that energy is removed in the form of latent heat, 
by converting a portion of that liquid to vapor. 

The analysis of single-phase heat removal is an inherently simpler 
task than that of two-phase heat removal. For purposes of preliminary 
analysis, one can characterize a single-phase coolant in terms of ''‘mean" 
properties: Particularly, by a mean heat transfer coefficient and a mean 
specific heat. Since energy is removed in the form of sensible heat, the 


latter, coupled with the bulk flow rate, determines the mean temperature 
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rise in the bulk fluid. For the present application, there is an obvious 
interest in keeping this rise as small as practical. Suffice it to say, that 
coupled with the generally low visosities of cryogenic coolants, this leads 
to sufficiently high bulk flow rates as to ensure fully developed turbulent 
flow of the coolant. One can thus determine the mean heat transfer coef- 
ficient experimentally, or utilize an appropriate predictive correlation, 
since such correlations are reasonably accurate under these conditions. 
(Stevenson, for example, based his design of a gaseous-helium-cooled 


solenoid (107) on the McAdams correlation: 


- 


0.4 0.8 


hD/k = 0.023P 13 oe aaa (24) 
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where h is the heat transfer coefficent; D, coolant tube diameter; k, 
thermal conductivity of the coolant; Ps Prandtl number of the flow, and 
R.,, Reynolds number of the flow. See Rohsenow and Choi (108) for a 
discussion of this and other such correlations. ) Given the mean heat 
transfer coefficient, the information required at this point can readily 
be determined. While it is apparent that considerable search may be 
necessary to determine the ''optimum'" operating pressure, etc., the 
procedure is conceptually straightforward. 

By contrast, the analysis of two-phase heat removal is not straight- 
forward: in general, predictive correlations are not sufficiently reliable 
as to permit accurate computation of the heat transfer coefficient. In 
fact, it would seem prudent--- if not mandatory--- to construct a precise 
model of any proposed heat transfer scheme, and obtain experimental 
confirmation of any theoretical predictions, before entering the final 
design stage. This would be particularly advisable in the design of large 


cryogenic saddle coils, as the size and complexity of the over-all system 
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would render any significant error in the anticipated heat transfer 
coefficient a most costly mistake. 

The above notwithstanding, most of the cryogenic coils which have 
been built, to date, have utilized some form of two-phase cooling. For 
this reason, the rest of this chapter, and most of the next, will be devoted 


to a general discussion of that problem. 


B. Pool boiling characteristics: 

In a physical sense, the simplest form of two-phase heat transfer is 
that commonly called ''pool boiling’, in which a heated specimen is immer- 
sed in a large pool of liquid, and cooled by natural-convection boiling, 
with geometries sufficiently ''open" to ensure negligible interference with 
convective flow. (These conditions can be approximated with small, 
‘jelly-roll'' solenoids, which are often cooled by simple immersion. ) 

The solid line in Fig. 19 is a qualitative representation of the sort 
of heat transfer curve which one might typically observe under such 
circumstances. To briefly discuss the regions of this curve: 

Region I represents non-boiling convective cooling, and is character- 
ized by low heat flux, and low heat transfer coefficient. (Note, however, 


tf 1! e ° ° 
that the AT referred to in Fig. 19 is TT sat’ As the bulk liquid is 


subcooled, the wall temperature is only slightly greater than the sat- 
uration temperature, and '' AT" is therefore small. ) 

Region II is the zone of "nucleate boiling’ heat transfer, and is 
characterized by the formation of bubbles at discrete points distributed 
across the heat transfer surface: i.e., at surface imperfections, such 
as scratches, pits, etc. As the bubbles grow and break away from the 
surface, they are replaced by relatively cool liquid swept in from the 
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"Typical Saturated Boiling Curves for Cryogenic Liquids 
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bulk of the fluid. This bubble agitation results in a rapid increase in 
heat flux with the number of "active nucleation sites, as the primary 
flow of heat is across the surface into the relatively cool liquid. Since 
surface imperfections tend to have a random distribution in size, the 
number of active sites is a function of AT, with 'large'' sites becoming 
active before ''small”’ sites. 

Region III represents the gradual transition from ''nucleate" to ''film" 
boiling. In this region, generally called the ''metastable film boiling" 
region, increasingly greater portions of the surface are blanketed by 
solid ''sheets'' of vapor. Under such ''sheets", the cool liquid cannot 
come in contact with the surface. Since the thermal conductivity of the 
vapor is generally much lower than that of the liquid phase, heat flux 
actually decreases with increasing AT, in this region. 

Region IV is the zone of fully developed ''film boiling". Here, the 
entire heat transfer surface is covered by a stable vapor film, and as 
AT increases, heat flux again starts to rise. Finally, at very high AT, 
radiation effects become significant, and heat flux can become quite large. 

The heat flux attained at the transition between regions II and III 


'in Fig. 19) is variously termed the ''peak", "critical", or 


(point ‘a’ 
"burnout" flux. The last term is quite appropriate. If heat is generated 
at a constantly increasing rate in, say, an electrically heated device, 
upon reaching this point on the nucleate boiling curve, the operating 

point must "jump" to the corresponding flux on the film boiling curve 
(point ''b' in Fig. 19), since the characteristics of the device would not 
permit the decrease in flux required for operation in region III. However, 
as indicated by the "typical scale along the abscissa, the AT required 


to maintain the same flux in region IV is generally so large that the device 


might well melt before the operating point could reach the film boiling 
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curve. 


C. Pool boiling data: 

To turn from the general to the specific, four particular cryogenic 
coolants will be considered in this study: nitrogen, neon, parahydrogen, 
and helium. All but neon have been widely used by various experimenters, 
and if current trends continue (e. g., see (109) or (110) ), liquid neon is 
likely to be widely used in the future. As coolants, these liquids effect- 
ively span the 0°-100° K temperature range of interest here. (One-atmos- 
phere boiling points are 77.4° K, 27.2° K, 20.3° K, and 4,2° K, respec- 
tively. Of course, some control of boiling temperature is possible in 
each case by variation of operating pressures. ) 

For basic information regarding boiling heat transfer to liquid nit- 
rogen or hydrogen, one is advised to turn to either of two excellent, 
recently published works: 

Boiling Heat Transfer for Cryogenics, by J. D. Seader, W. 8S. Miller, 
and L. A. Kalvinskas (111), is a complete, 171-page survey of the liter- 
ature on heat transfer to liquid hydrogen, oxygen, and nitrogen. This 
work discusses and evaluates both predictive correlations and experimental 
data, for pool boiling, forced-convection subcooled boiling, forced-con- 
vection saturated boiling and forced-convection film boiling to these liq- 
uids, and includes a twenty-page bibliography, as well as eighteen pages 
of physical property charts. The work is particularly valuable for its 
coverage of the largely neglected topic of cryogenic forced convection. 

Nucleate and Film Pool Boiling Design Correlations for O92, No, Hg, 
and He, by E. G. Brentari and R. V. Smith (112), presents an even more 


detailed survey of pool boiling data for nitrogen and hydrogen, and, as is 
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apparent, covers helium as well. The authors discuss the various fact- 
ors which contribute most directly to ''scatter'' in observed data, and 
conclude that the general nucleate heat transfer correlation and maximum 
nucleate heat transfer correlation of Kutateladze (113) provide the best 

"a priori'' design standards, when details of surface condition, orientation, 
geometry, etc. are not known. Specifically, these correlations state: 
fe) = 8.54 x 10° %(e Jap)" °° A872 mana ca | 626 2.90 


(AT ) (25) 


(k — p 

and 

(a/A)nu pnout =. lot ooo "(gto)" x (26) 

where all properties pertain to the liquid phase, except those with sub- 
script ''v'', and g/A is heat flux in watts/em?; Che specific heat at con- 
stant pressure, joules/gm-°K; h,,» latent heat of vaporization, joules/ 
gm; fp: density, gm / cm”: k, thermal conductivity, watts/cm-°K; P, 
pressure, dynes/cm’; t, surface tension between liquid and vapor, dynes/ 
cm; [{, viscosity, gm/cm-sec; AT in se and g is the acceleration due 
to gravity in om/sec’. (The bibliography provided in their work includes 
a section listing possible sources of information regarding thermodynamic 
properties.) More conveniently, Brentari and Smith provide a comparative 
plot of these correlations and available data, and a second plot of pre- 
dicted q/A vs. AT, for various pressures (in atmospheres), for each 
liquid. 

Unfortunately, very little experimental information is available at 
this time regarding the heat transfer characteristics of liquid neon, prob- 
ably due to the hitherto limited supply. Hopefully, some data may soon 
be available as a by-product of the testing of the NASA-Lewis magnet 
(21,22), but the only experimental report known to the author at this time 
is that of Bewilogua and Mahn (114). However, one can estimate the heat 


transfer characteristics from eqns. (25) and (26), and in fact, Frederking 
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(1%5) does present a plot of q/A based on eqn. (26). (The National 


burnout’ 
Bureau of Standards appears to be currently updating its information re- 
garding the thermodynamic properties of neon. See, for example, (116) 
and (117). ) 

In any event, based on (112) and (115), one can construct a quick 
estimate of the nucleate boiling heat transfer characteristics of the four 
specified coolants, as follows: Assume the one-atmosphere "burnout" 
flux for nitrogen is 19 watts/cm”, and that is occurs at a AT of 12° K. 


at 3° K; for hydrogen, 


Similarly, for neon the figures would be 12 watts/cm 
8.5 at 3° K, and for helium, 0.77 at 0.4° K. Plot these points on a log- 
log chart, and pass a line downward from each, ata slope of 2.5, to 
obtain the one-atmosphere characteristic. (Of course, since thermo- 
dynamic properties are a function of pressure, both eqns. (25) and (26) 
are functions of pressure. In general, these characteristic lines shift 
toward lower AT's as pressure is increased. The locus of the burnout 
point, as a function of pressure, is (roughly) a concave-downward arc. 
Helium exhibits its maximum burnout flux at about one atmosphere. For 
the other liquids, the maximum is roughly: nitrogen, 37 watts/cm” at 

AT = 6.6°K (p = 10 atm. ); neon, 18.6 atl. 3°K (~8,5 atm.), and hydrogen, 
13.2 at 1. 6°K (4 atm. ). 

Again, it must be stressed that this is a rather crude "first estimate”. 
Equation (25) is a reasonable, slightly conservative ''mean', and is prob- 
ably an excellent guess for average’ conditions. However, conditions 
may deviate considerably from ‘average’, as evidenced by the fact that 
results reported by some experimenters lie more than a factor of 5 above 
or below this correlation. For this reason, it may be worthwhile to 


briefly discuss a few of the factors which lead to deviation from the norm: 
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Surface condition: As mentioned above, heat flux varies with the num- 
ber of active nucleation sites, which among other things is a function of | 
the size distribution of surface imperfections. Thus, a mirror finish 
would hinder nucleate boiling heat transfer, and data points taken from 
specimens with such a finish fall well below the norm. Conversely, an 
artificially roughened surface improves nucleate boiling heat transfer. 
(See, for example, (118).) Further, the number of active sites is a 
function of the extent to which a particular liquid ''wets" the surface. Any 
coating or contaminant which decreases ''wettability'' of the surface will 
tend to increase heat flux. Cowley, et al. (119) report that artificially 


greased surfaces can be driven to 50-100” 


higher flux levels before 
burnout, though such surfaces exhibit "below normal" flux levels at sub- 
critical AT's. Here, the dominant factor is the poor thermal conductivity 
of grease. The presence of such an “insulating layer'' (~ 0.2 mm thick) 
causes a four- to ten-fold increase in total AT (between the bulk liquid 

and the surface) at burnout. The greatest portion of this total AT occurs 
in the insulating layer, and oddly enough, results in a temperature pro- 
file which drives a greater heat flux through the insulation, into the liquid, 
than could be driven directly from the surface into the liquid. The authors 
state that such use of insulation should increase maximum heat flux in 
any situation in which heat transfer is initially limited by conditions on 
the liquid side of the interface. 

Surface history: A given surface will generally exhibit greater heat 
flux when first immersed in coolant, than after repeated measurements 
under the same conditions. This is partially attributable to the boil-off 
of surface contaminants, but is probably more dependent upon the initial 
presence of adsorbed gases. Haselden and Peters (120), having taken 


pains to ensure an initially clean surface, found that such transients per- 
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sisted for as long as 2 to 3 hours. Experimenters who do not choose to 
wait for such an interval, to ensure that consitions have stabilized, may 
find that their data predict higher heat fluxes than can be maintained for 
steady-state operation. 

Surface orientation and geometry: Although nucleate pool boiling is 
not particularly sensitive to device geometry, any features which alter 
the circulation pattern associated with ''normal" natural convection (and 
therefore cause deviation from the strict definition of ''pool boiling’) will 
have some influence on nucleate boiling heat transfer. For example, con- 
fined channels may tend to become ''clogged'' with vapor, to such an extent 
that the churning action associated with nucleate boiling is inhibited. This 
would sharply reduce attainable heat flux. Similarly, heat transfer from 
a ''bottom'" surface may be hampered in a gravity field, since bouyancy 
may impair the normal motion of bubbles away from the surface. Vari- 
ations in convective circulation can also cause minor differences between 
measurements taken near the ends of a test specimen, and those taken 
near the center. However, by far the most critical of these geometric 
factors is the marked decrease in attainable heat flux in confined channels. 
Since this is precisely the type of geometry anticipated in the cooling sys- 
tems of large coils, natural convection cooling is not an attractive solution 


for such applications. 


D. Qualitative influence of forced convection: 

The obvious solution to the problem of "confined geometry’ coolant 
ducts is to turn to forced convection. In addition to the ''microscopic'' 
processes discussed above, this adds a new form of turbulent mixing to 
the heat transfer process: i.e., the ''macroscopic'' circulation of forced 


flow. Qualitatively, the effect on the nucleate boiling curve is as depicted 
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in Fig. 19. In general, the increased heat flux attributable to forced flow 
can always be made to offset--- or more than offset--- the decrement in 
attainable flux due to confined geometry. At low qualities, forced con- 
vection produces a sharp increase in both heat transfer coefficient, and 
burnout flux. As quality increases the heat transfer coefficient increases 
even more (as shown in Fig. 20a), until transition to film boiling occurs, 
after which the heat transfer coefficient decreases sharply. The forced 
convection burnout flux decreases more or less linearly with quality, as 
shown in Fig. 20b. However, even at high qualities, one might reason- 
ably expect the forced convection burnout flux to exceed the pool boiling 
burnout flux. 

The heat flux attained under forced flow conditions is often concep» 
tually divided into its two component parts: one attributable to pool 
boiling processes, and the other attributable to forced convection. Actu- 
ally, the physics of the problem are not quite that simple. As the Reynold's 
number of the flow is increased, nucleation is gradually suppressed. As 
the quality (mass fraction of vapor) of the flow is increased, the flow of 
relatively large volumes of vapor along the tube increases the convective 
turbulence attributed to forced flow. Eventually, a point is reached where 
the contribution of microscopic" heat flux is entirely negligible, but by 
this point, the ''macroscopic'' heat flux may be much greater than that pre- 
dicted from eqn. (24). Seader, et al. (111) present a detailed discussion 
of the predictive correlations which attempt to cope with these interrelated 
complications, but in general, such correlations have not met with marked 
success, when applied to ''conventional'' coolants. Data for forced con- 
vection saturated boiling heat transfer to cryogenic coolants is so scarce 
that pne can hardly say the various correlations have been fairly tested, 
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FIGURE 20a 


Qualitative Variation of Heat Transfer Coefficient 
with Quality 


FIGURE 20b 
Qualitative Variation of Critical Flux with Quality 
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In any event, it would appear that the most successful such correlation, 

to date, is that developed by Chen (121): 
q/A = (hooF + hpp») AT, (27) 

where Apc is the forced convection heat transfer coefficient, which could, 
for example, be computed from eqn. (24); hpp is the pool boiling heat 
transfer coefficient, which could be computed by dividing eqn. (25) by AT; 
F is an empirical amplification factor which corrects for increased heat 
flux due to increased convective turbulence at high quality, and S is an 
empirical attenuation factor, which corrects for the suppression of 
nucleation at high liquid flow rates. Both (111) and (121) present plots 
of F as a function of the Martinelli two-phase parameter, X andSasa 


te 
function of RE 25 where 
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(here, x is quality, and the other parameters are as defined above), and 
Ret DG(1-x)/I (29) 

(where G is the mass flow rate divided by the tube cross sectional area). 

As a practical matter, the plots provided in (111) are, perhaps, somewhat 

easier to work with. 

It is apparent that the use of this correlation to predict heat flux ata 
Single point along a coolant duct is no minor task, unless the problem has 
been prepared in advance for computer solution. Further, since both 
eqn. (28) and eqn. (29) are functions of quality, it is necessary to evaluate 


"mean" flux is to be deter- 


eqn. (27) at many points along such a duct, ifa 
mined for purposes of preliminary analysis. Chen has tested this cor- 
relation against available data for vertical, axial, stable flow, at heat 


fluxes below the critical, for water and various hydrocarbons, and found 


an average deviation of only + 11 percent, for a quality range of 1 to 71 
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percent. This achievement is truly remarkable, considering the com- 
plexity of the problem! However, individual points may deviate consid- 


average’, and forced convection heat transfer is 


erably more than the 
more sensitive to geometry and orientation than pool boiling heat transfer. 
Considering this, and the scarcity of data for two-phase forced convection 
heat transfer to cryogenic coolants, it should be stressed again that experi- 
mental confirmation of heat transfer characteristics should be obtained 
before proceeding into the detailed design phase. 

Under the above circumstances, it is doubtful that the considerable 
labor involved in the use of Chen's correlation is justifiable for purposes 
of preliminary analysis. Instead, it is recommended here that the pro- 
cedure for determining pool boiling heat transfer characteristics, as 
described in the preceding section, be applied with no correction for 
forced flow. This procedure has the advantage of being quicker, simpler, 
and conservative, since the forced convection heat transfer characteristics 
should be as good as, or better than, this estimate. (It is noted that 
Frederking (122) specifies that heat flux under conditions of forced con- 
vection will always exceed predictions from a pool boiling analysis if 
quality is low. However, Chen's correlation is less restrictive than this, 
Since fluxes predicted from his model exceed the pool boiling figures up 


to the limiting quality of 70%. ) 
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VI. COOLANT FLOW 


A. Geometric considerations: 

While heat flux considerations may determine the amount of heat trans- 
fer surface required for a particular coil, the specification of cooling sys- 
tem geometry is far from complete. The required surface area could be 
provided in a single large coolant duct, or many small ducts, and the choice 
between these extremes must be based on considerations other than heat 
flow alone, such as the resulting packing factor, pressure drop across the 
coil, thermal gradients in the winding, ease of assembly, etc., etc. 

To deal with the last item, in particular: The natural tendency would 
be to assume some specified distribution of circular tubes, of some given 
L/D. This would be motivated in part by historical precedent, and in part 
by convenience: Most conventional heat transfer devices utilize circular 
tubes; we are accustomed to thinking in terms of circular tubes, and the 
great majority of the correlations used in heat transfer and coolant flow 
analyses are based on circular tube data. Unfortunately, circular tubes 
are not likely to be used in ‘real’ cooling systems for coils of the type 
being considered here, for precisely the same reason that they are used 
in SO many other applications: i.e., ease of manufacture. 

Whatever provisions are made for coolant flow must conform to the 
realities of coil construction. For example, it was shown in Chapter III 
that the internal support structure required to support the conductors 
against magnetic stress must run axially down the length of the cylindrical 
lobes. Thus, it is reasonable to assume that the coolant ducts must run 
axially: any other geometry would require the ducts to penetrate such 
structural members, which would be highly undesirable, if not impossible. 


Further, cooling ducts cannot be run through regions in which they would 
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be subjected to a heavy stress load. This would require that they have 
thick walls, for mechanical strength, which would extract a double pen- 
alty: reduced heat transfer, and reduced packing factor. Since structure 
has already been provided to support the conductors, it is reasonable to 
assume that coolant ducts would be located contiguous to that structure: 
stress levels are at a maximum on one Side of each support shell, but 
total stress is virtually zero on the other side. A coolant duct, run along 
the ''no-load" side of such a shell, would require a minimum of structural 
support. This fact is conducive to the use of coolant slots, vice tubes, 
or pipes. Of course, the model described in Chapter III calls for a greater 
number of support members in regions of high magnetic fields, since these 
are the regions of greatest magnetic force density. These are also the 
regions where magnetoresistance, and therefore heat generation, would 
be greatest. It is convenient that a greater number of coolant slots would 
be provided automatically in this region, under the proposed scheme. 
Thus, the basic geometry considered here will involve the use of 
axial slots for coolant flow, with one side of each slot bounded by the 
internal support structure of the coil. It will be assumed that the height 
of such slots would be much less than all other significant dimensions 
(width, length, and radius of curvature), and coolant flow will therefore 
be handled on a ''per unit width'' basis. Since the most likely structural 
material is stainless steel, and since the thermal conductivity of stain- 
less steel is two to four orders of magnitude less than that of the electrical 
conductors considered here, at the temperatures of interest, it will be 
assumed that heat enters each coolant slot only across the surface away 
from" the support shell. (So far as it involves the results of this analysis, 
one can treat the possibility of heat input through both surfaces by simply 


doubling the attainable heat flux which would otherwise be estimated from 
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the preceding chapter.) Again, the purpose of this analysis will be to deter- 
mine the relationships between coolant properties, pressure, quality, and 
heat flux, as ‘'input'' variables, and packing factor, coolant pump power, 


and axial temperature variation, as ‘output’ variables. 


B. _ Properties of coolants: 

Since coolant properties are a necessary input to this analysis, it 
would be well to consider them before proceeding further. In particular, 
since the primary effort here will be devoted to an analysis of two-phase 
cooling, the saturation properties (i.e., the thermodynamic properties of 
a saturated liquid and its vapor, at least nominally under equilibrium con- 
ditions) are of especial interest. In fact, an interesting problem exists 
regarding the processing of this data. Because the coolants considered 
(helium, parahydrogen, neon and nitrogen) will be examined at pressures 
approaching their respective critical pressures, such customary and con- 
venient concepts as "incompressible" liquids, and ''perfect'' gases, are 
not generally applicable. Further, many pertinent properties exhibit 
sufficiently non-linear behavior, within the range of interest, as to raise 
doubts regarding the accuracy of such techniques as linear interpolation 
between tabulated values. As a convenient means to avoid these compli- 
cations, a data fitting routine was used to generate a least-square poly- 
nomial fit, expressing each such property of interest as an explicit func- 
tion of pressure, in atmospheres. Specifically, the properties of specific 
volume of saturated liquid and vapor, saturation temperature, specific 
enthalpy of the saturated liquid, and heat of vaporization, were expressed 
as polynomials of the below form: 

1/0 jo) Bs 1 Ele a aap - ap” 


2 3 4 
1/0 = by/p + by + bgp + bgp” + bsp” + bep 
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TABLE 6 


SATURATION PROPERTY COEFFICIENTS 





Property Helium Parahydrogen Neon Nitrogen 
a, = 5. 45596 13. 3964 0. 767861 1. 18952 
1/0, 422 > 5. 09263 0. 796955 0. 0675544 0. 0559020 
lla, = 3.69782 -0,0606407  -0,00664262 -0, 00437871 
ay = 1, 15685 0.00500739 0,000224027 0. 000195317 
Deviation = 0. 058% 0. 041% under 0, 01% 0. 074% 
b, = 49. 3439 654, 454 101,175 195, 349 
by = 31,7682 128, 968 0. 394318 29. 5132 
1/0 (23 = ~32- 7486 -41, 6982 3, 23665 -8. 96109 
Py by = 15. 3537 6. 73209 -0, 949753 1, 49260 
be = -4, 01758 -0, 563735 0.0978116 -0. 126834 
b° = 0, 281743 0.0181262  -0,00343717 0. 00420598 
mevietien= under 0, 01% 0. 100% 0.561% 0. 059% 
c, = 2.67907 16. 4517 23. 1925 67. 7503 
cj = 2.04654 4.71658 5. 04433 11, 9973 
7 Jeg = -0. 461891 -1. 00566 -1, 25994 ~2, 66612 
cy = -0, 180713 0.144354 0. 207 938 0. 392040 
ce = 0. 158552 -0.0108406 -0.0176176 -0, 0300996 
co = -0, 0309227 0. 000323253 0,000582154 0. 000918179 
Deviation = under 0.01% 0.049% 0, 038% 0.035% 
d, = -28,9212 -293, 057 -2, 09216 9.57010 
d, = 139, 468 43. 7316 8. 43141 24. 7498 
», (dg = 7211. 281 -7, 27746 =1, 56955 -5, 51273 
L\dy = 163,231 1. 05418 0. 201394 0. 820578 
dz = -61, 7927 -0.0791551 -0,0130336 -0, 0632501 
d- = 9, 19968 0.00245513 0.000326842 0, 00193276 
Deviation = under 0.01% 0.033% 0.721% 0.134% 
e, = 58,3544 462. 481 91. 6010 214,510 
e, = -137, 426 -14, 4022 -7,50772 -18, 9988 
4 e@ = 211,270 ~1, 53707 1, 64718 4.77690 
LV ey = -163,972 0. 159357 -0, 200480 -0, 789073 
ex = 61,7927 -0, 00838393 0,0130336 0. 0632501 
35 = -9, 19968 0.00245513 -0, 000326842 -0. 00193276 
Deviation = 0, 035% 0. 015% 0. 693% 0. 148% 
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T = Cy 17 CoP * cP. + c 4p + C5. + cg 

hy = d, + dop + d.p + dP + d.p + dep 
hy, =e, + e5p + ep. + e 4p” + ep + eee” 
where specific volumes are given in em°/gm, temperature in ne and 
enthalpies in joules/gm. The coefficients a, b. , etc., are tabulated 
below. The particular functional forms indicated ensure, in all cases, that 
the accuracy of ''fit'’ exceeds the assumed accuracy of the flow models 
developed here, within the pressure ranges investigated: i.e., 1-10 atmos- 
pheres for saturated nitrogen, neon, and parahydrogen, and 1-2 atmos- 
pheres for saturated helium. (The lower limit of 1 atmosphere was set 
because it was felt that considering the hazard of leakage, no one would 
care to design a cryogenic cooling system, for a working coil, to operate 
at sub-atmospheric pressures. The upper limits were set, as a practical 
matter, considering the decrease in heat of vaporization, and increase in 
coil operating temperature, at increased operating pressures.) The ‘‘devi- 
ation’ indicated for each set of coefficients is the mean of the absolute 
value of the difference between ‘input’ property values, and property val- 
ues generated from the indicated polynomial. 

One other general point might be mentioned before turning to the par- 
ticular: If two-phase cooling is contemplated, one would presumably 
attempt to operate to a relatively high outlet quality, since pumping ‘unused’ 
liquid through the cooling loop increases pump capacity, and increases duct 
size, while contributing nothing to the removal of heat from the coil. In 
practice, outlet quality must be limited by heat transfer considerations 
(more on that, later!), but if attainable qualities even approach the range 
indicated by the upper limit of Chen's correlation, it should be recognized 


that the outlet stream, from a volumetric viewpoint, will be almost entirely 


vapor: a ''typical'' figure would be 99% vapor, for a quality of 50%. As 
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vapor densities tend to be relatively low, momentum considerations indi- 
cate that relatively high vapor outlet velocities may result from relatively 
low pressure drops: i.e., it is quite possible that the flow could ''choke’ 
at, or even before, the outlet. If this is to be avoided, outlet velocities 
must not be allowed to approach the sonic velocity in the stream. (Note 
that this would be the so-called "sonic velocity of Mach 1", which Shapiro 
(123) denotes "o%*"', not the velocity of sound under stagnation conditions. ) 
The problem here is to determine just what the limiting velocity should be. 
The only available data are for Co the sonic velocity under stagnation 
conditions. Shapiro shows that for perfect gases: 
cx/c. = V2 7+ v, 

where T = Salle = the ratio of specific heats. However, we are not deal- 
ing with perfect gases, and in any event, the sonic velocity in the mixed- 
phase coolant stream would presumably be lower than that in the vapor, 
under conditions of 100% quality. Considering these uncertainties, a 
"limiting velocity'’ has been defined for each coolant, by using the best 
available information for Co and [at the one atmosphere saturation point, 
computing c* as above, and multiplying the result by roughly 0.8. Inthe 
subsequent analyses, combinations of parameters which would result in 
flow velocities greater than this limiting velocity will not be considered. 
The values of T; Co and Viim 2re presented below, along with other spe- 
cific information of interest for each coolant. 

Helium: (Property values determined form (102,110) and (124-127). 
Saturation property polynomials developed from (127). ) 


one-atmosphere boiling point: 4. 22°K 


critical point: 5. 25°K at 2.26 atm. 
ratio of specific heats: ~1,.74 

sonic velocity: ~145 m/s 

limiting velocity: 100 m/s 


Considering its low heat of vaporization, and the high cost of 4. 2°K 
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refrigeration, liquid helium is not likely to be used as a two-phase coolant 
for cryogenic coils, except when cryogenic conductors are used to ''stabi- 
lize'' superconducting coils. (See, e.g., (14). ) However, supercritical 
helium is the logical coolant for coils designed to operate in the vicinity 
of 10°K (24,107). The National Bureau of Standards indicates that the 
optimum operating pressures for such applications lie in the band 15-20 
atm. (105). The required coolant volume fraction is not a strong function 
of pressure, in this range, since the increase in density as pressure is 
increased from 15 to 20 atm. is largely offset by a decrease in specific 
heat. 

Parahydrogen: (Property values determined from (102, 110, 124-126) 
and (128-131). Saturation property polynomials developed from (131), 
except that for liquid specific enthalpy, which was developed from (130). ) 


one-atmosphere boiling point: 20. 27°K 


critical point: 32.98°K at 12.76 atm. 
ratio of specific heats: Love 

sonic velocity: 357 m/s 

limiting velocity: 250 m/s 


Hydrogen gas occurs in two molecular forms, orthohydrogen and para- 
hydrogen, distinguished by the relative direction of nuclear spin of the two 
component atoms of the molecule. At room temperature, the equilibrium 
composition is 75 per cent ortho, 25 per cent para. However, at 20. Soke 
the equilibrium concentration is about 99% parahydrogen. Since the energy 
released during conversion is greater than the latent heat of vaporization, 
this reaction must be catalyzed to completion if liquid hydrogen is to be 
used continuously in closed-loop cooling systems. Presumably, the com- 
position would not change significantly during the brief period required for 
the liquid to pass through the coils of a magnet. For these reasons, the 


properties discussed here are those of liquid parahydrogen. 
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Basically, liquid hydrogen is a highly promising coolant for cryogenic 
coils, because of its availability and low boiling point. Although some con- 
cern has been expressed regarding the safety of this particular coolant, 
those who have worked with liquid hydrogen (notably, (19, 20) report "no 
problems", provided reasonable precautions are observed. (See NBS 
Report 7253 (132) for a list of safety precautions. ) It would seem reason- 
able to presume that the hazards involved in the use of this coolant in 
closed-loop systems would be less severe than is the case with open-loop 
systems, where large volumes of hydrogen gas must be stored or vented. 

Neon: (Property values determined from (102, 110, 133) and (116). 
Saturation property polynomials developed from (116).  ) 


one-atmosphere boiling point: 27. 09°K 


critical point: 44, 40°K at 26.19 atm. 
ratio of specific heats: ~1.7 

sonic velocity: 134 m/s 

limiting velocity: 90 m/s 


As has been discussed above, liquid neon has not been widely used as 
a cryogenic coolant, to date, because of its hitherto limited availability 
and high cost. Specific advantages attributed to this liquid as a coolant 
are a relatively high latent heat, on a per-unit-volume basis, and ease of 
temperature control through pressure regulation. In addition, liquid neon 
is the only two-phase coolant available which permits maintenance of tem- - 
peratures near|the liquid hydrogen range, without the hazards associated 
with liquid hydrogen. 

Nitrogen: (Property values determined from (102, 124-126, 128, 129) 
and (134). Saturation property polynomials developed from (134). ) 


one-atmosphere boiling point: 77. 4°K 


critical point: 126.1 K 
ratio of specific heats: 1,44 

sonic velocity: ~167 m/s 
limiting velocity: 120 m/s 


Nitrogen might be described as the ''old standby" cryogenic coolant, 
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and apparently is used in many applications primarily because its proper- 
ties are well known, it is readily available in large quantities, and it can 
be handled with relative ease. However, although various ''laboratory'' 
coils have been tested in liquid nitrogen, Tables 1-4 indicate that the tem- 
peratures attainable with this coolant are too high to permit realization of 
the full potential of cryogenic magnets. It is unlikely that large ''working"' 


coils would be designed for two-phase liquid nitrogen cooling. 


C, Single-phase flow: 

Aside from the concepts discussed in the preceding chapter, the only 
design factor essential to the preliminary analysis of single-phase cooling 
is the question of pressure drop across the coil. The customary manner 
of estimating this pressure drop would be to do so in terms of an effective 
friction factor, f: > 

p = 4£ (L/D Ov" /2), (30) 

where one would again utilize ''mean' property values, and where L is the 
estimated duct length, which would presumably be of the order of the MHD 
channel length (or, perhaps, the product of that length times the form fac- 
tor, defined in eqn. (6) ). For non-circular cross sections, D would be 
the "equivalent" or "hydraulic'' diameter, given by: 

D=4A/P (31) 
where A is the duct cross-sectional area, and P is the wetted perimeter. 
(For a rectangular slot, of the type discussed previously, 4A~4ws and P 
*>2w, so D = 2s, since slot height s » slot width w. ) 

Various estimates are available for the value of f, given fully developed 
turbulent flow conditions. Rohsenow and Choi (108) suggest: 


0,25 3 4 


FS), 0791/R, , 3x10 <R.=2 x 10 (32) 
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or (32) 


0.046/R," °°, R, 22x 10° 


f 


Stevenson (107), on the other hand, uses: 


32 


fl (33) 


i 


0.0014 + 0; 125/R,"" 
(In both eqns. (32) and eqn. (33), R. = OvD/ LL .) 

Having established an estimated bulk flow rate (from mean temperature 
rise considerations), the above Ap can be used to determine the power 
required to pump the coolant through the loop, since 

P= Fv =(AAp)v - wAp/g (34) 
where w is now the mass flow rate. Of course, this term represents 
energy added to the coolant stream, which is converted to heat by friction, 
and is therefore an added refrigeration load. However, one would antic- 
ipate that this load would normally be negligible, compared to ISR heating. 
Eqn. (34) can be used to verify this fact, as well as to ''size'' the pump 
required for a particular coil. Once again, the procedure is relatively 
straightforward, and a great deal simpler than the analysis of two-phase 


flow. 


D. Two-phase flow: 

Considering the lack of experimental data for forced-convection two- 
phase cryogenic cooling, the information required for coolant flow calcu- 
lations must be generated from an analytical model. The problem becomes 
to select that model which is most likely to correspond to the physical 
realities of the Situation, under these circumstances. Two possibilities 
will be developed and discussed here, but the ultimate choice must await 
additional testing. 

To turn, first, to background considerations: Sydoriak and Roberts 


(128) conducted one of the few studies of two-phase cryogenic coolant flow 
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in confined channels, for both natural and forced convection, and report 
close correlation between their data and "calculations based on a friction- 
less flow boiling theory'’. In essence, the pressure drop due to the always- 
low viscosity of cryogenic fluids was totally dominated by that attributable 
to acceleration of the vapor phase (as required to maintain a given flow 

rate when a quantity of relatively dense liquid is converted to relatively 
light vapor). The pressure drop associated with any such "inviscid" accel- 
eration can be computed from momentum considerations, provided certain 
basic information is available regarding flow conditions in the duct. 

Chen (121) classifies the possible flow conditions, as depicted in Fig. 
Bil: 

Zone 1, at the inlet, represents subcooled conditions, as would be 
encountered when the saturated liquid, from a heat exchanger, passes 
through a pump, with resultant increase in pressure, but little increase 
in temperature. In this region, there is either no generation of vapor, or 
formation of bubbles, at the duct wall, which decrease in size as they depart 
the thermal boundary layer. 

Zone 2Z2 represents the region of first net generation of vapor. Here, 
bubbles are formed at the heat transfer surface, break away, and move 
into the (now saturated) liquid phase. Quality is low, and aside from the 
convection down the duct, conditions are rather similar to ''nucleate pool 
boiling". 

By zone 3, sufficient vapor has been generated that the individual 
bubbles tend to merge, forming relatively large ''blobs" of vapor which 
accelerate down the duct. Nucleation is gradually suppressed. This flow 
condition is characterized as ''slug flow". 

Eventually, the center of the duct becomes filled with an unbroken 


stream of vapor, as in zone 4, and an effective separation of the phases 





21 


FIGURE 
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occurs. As Chen's work dealt with circular tubes, he called this the 


"annular flew 


regime. Although some liquid is carried with the vapor 
in mist-like droplets, the bulk of the liquid phase moves along the walls 
of the tube. 

Finally, as quality increases, the tube walls gradually become dry 
and zone 5, the ''liquid deficient zone'', is entered. All the remaining 
liquid, in this zone, is swept along in the form of droplets, dispersed 
through the vapor phase. 

Actually, the scale of Fig. 21 is distorted. Chen's correlation was 
based upon the assumption of annular flow conditions, and as mentioned 
above, is held to be valid, generally, within the range 1 to 71 per cent 
quality, for the vertical, circular-tube geometry which he investigated. 
However, the fact that he determined this to be true for one geometry is 
no assurance that it will be true for other geometries. The various re- 
gimes may be less clear-cut, with some zones missing, others enlarged, 
or mingled. The sequence depicted here is "idealized, and deviations 
from this Sequence are important in two respects: first, flow conditions 
in zone 4 are most nearly typified by separate, inviscid flow of the liquid 
and vapor phases, with no mixing between the phases. Zone 5, and per- 
haps zones 1 and 3, are better typified as inviscid flow, with the boundry 
condition of ''no slip'' between the phases. Thus, any shift in the relative 
importance of these zones would presumably require a shift in the analyt- 
ical model used to represent flow conditions. Second, when the tube walls 
become dry, the heat transfer coefficient drops abruptly, due to the low 
thermal conductivity of the vapor phase, which then contacts the heat trans- 
fer surface. If zone 5 is entered at qualities significantly below 70%, for | 
a particular geometry, the outlet quality for a coil cooled by such a system 


would perforce be limited to such lower qualities. In other words, the fact 
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that Chen has proven that heat transfer improves with quality if annular 
flow conditions are maintained, does not, a priori, guarantee that such 
conditions will be maintained. Again, conditions must be confirmed for 
the particular application contemplated. 

Bearing this in mind, two inviscid flow models will be examined 
below: 
No slip (mist flow) model: Of the two limiting conditions discussed above, 
the simpler, from an analytical viewpoint, is that assuming no slip between 
the phases. If either phase is assumed to be dispersed through, and moving 
with, the other, the resultant two-phase flow can be treated as a hypo- 
thetical, homogeneous fluid with properties which are appropriately aver- 
aged between those of the liquid and the vapor. Specifically, let us assume 
steady, inviscid, one-dimensional flow, traversing a slot of constant, but 
as yet unspecified height, of unit width, and of some length L. The vari- 
able z will be used to denote distance from the inlet of the slot (x will refer 
to quality). It will be assumed that the flow has just reached zero quality 
at the slot inlet (i.e., x=0 at z=0),and is subjected to a constant heat flux, 
q/A, between z=0 and z=L. (This last is equivalent to the assumption that 
heat transfer to ''real'' liquids, which is known to be a function of quality, 
can be expressed in terms of an effective ''mean' q/A.) With an eye toward 
eventual applications, it will be assumed that outlet pressure, and desired 
outlet quality, are specified, and that we desire to investigate the relation- 
ship between inlet pressure (or, equivalently, pressure drop across the 
slot), inlet velocity, and slot height. The axial distribution of temperature, 
velocity, etc. will also be of interest, and in particular, the solution for 
each combination of parameters must be tested to ensure that the outlet 
velocity does not exceed the limiting velocity, defined in Section B of this 
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chapter. 

An iterative approach will be used to accomplish these objectives. In 
order to investigate the influences of inlet pressure, that parameter will 
be systematically incremented, as: 

P, =P ,t+nAp n= 1,2, 3----- 

By use of the saturation property polynomials, all pertinent thermo- 
dynamic properties at the inlet and outlet can be determined, once Py and 
P, have been specified. However, the ''appropriately averaged’ density 
of the homogeneous’ fluid is required, to relate these variables to flow 
rate. To quantize this concept: one gram of coolant, in the duct, is by 
definition comprised of x grams of vapor, and (l-x) grams of liquid. The 
vapor occupies a volume of x/ OQ. cubic centimeters, and the liquid, (1-x) 
ip 1° since the total volume occupied by this gram of coolant is, by defi- 
nition, the specific volume of the hypothetical fluid, the "appropriately 


averaged’ density is defined as: 


1 

p= TD, * T-=T OD, (35) 
Since we are considering a slot of unit width, continuity implies: 

Ww = PvA = Pvys = const., (36) 
where w is the mass flow rate, v the common velocity of the two phases 
(or of the homogeneous" fluid), and s the still undetermined slot height. 
However, since s is constant, 

Pv = const. = Di; (37) 


'i'' refers to inlet properties). Under the 


(here, as elsewhere, subscript 
stipulated flow conditions, the momentum equation can be integrated di- 
rectly, between inlet and outlet, to yield : 


“Ap = p, - P; = (0,v,) AV = (Pyy,)°O/ 5-1/0) (38) 


Since the only unknown in the above is v,, we can immediately state: 
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v7 1/9, )VAIC/O | - 1/94) (39) 


The same relationship can also be solved for Vo? which can immedi- 
ately be checked against the specified limiting velocity: 


Vo Vy PL yyy (40) 


Assuming v_ < v 


O Ene aa proceed: Since inlet quality was stipulated 


to be zero, and the slot is of unit width, at a point infinitesimally removed 


from the inlet, by a distance dz, we have: 


ee (FA) dz (41) 
lv 


Here, we encounter two problems. First, although Ay is a known function 
of pressure, the pressure distribution along the slot has not yet been 
determined, so eqn. (41) cannot be integrated directly to determine x as 

a function of z. Second, though the ratio (q/A)/w is known to be constant, 
s is not known in eqn. (36), and the value of this constant is therefore 


unknown, Letting the unknown constant be K, one additional item is known: 





L 
wh h Oo 
lv lv 
O 
where ig 

fi | -1i ;_dz 
h Lf{h 
ley 9 iV 


Of course, (1/h,_) is not known, either, but the problem is now in 


convenient form for an iterative, numerical solution. Assume 

BE i fa 1 

F-] oe | Ny : My | a 

i re) 

By use of eqn. (43), eqn. (42) may be solved for a ''first guess’ for K. 
Divide the slot into N segments, of length Az= L/N. If N is sufficiently 
large (N =10 yields reasonable results), one may assume P= Peh, = sy 
etc., throughout the first segment. Eqn. (41) may then be used to determine 
the quality at the end of the first segment, after which eqn. (35) can be used 


to increment for then eqn. (37) to increment v, and eqn. (38) to increment p. 
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With p incremented, new values may be determined for the thermodynamic 
properties, and these used to advance the solution through the next segment, 
etc. When the outlet is reached, computed outlet quality can be compared 
to desired outlet quality, and the estimate for K corrected accordingly. 

A simple, linear correction was found to yield rapid convergence: 


Ie K (x lemere 


2 


Of course, when suitable convergence is obtained, the desired tem- 


Heed! <compt 


perature, pressure, and velocity profiles have also been obtained, and 
eqn. (42) may be solved directly for the hitherto unknown relationship 


between (q/A) and s. 


No mixing (annular flow) model: Essentially, the assumptions made regard- 
ing flow conditions for the no-slip model pertain as well to the no-mixing 
case, except that concerning the distribution of the phases. Here, we 

Shall assume that the liquid phase remains adjacent to the walls of the slot, 
while the vapor phase flows along the center of the slot. With neither phase 
entrained in the other, the assumption of inviscid flow implies that liquid 
and vapor can (and indeed, must!) proceed at different velocities. In- 
stead, if the pressure profile is to be one-dimensional in character, the 
bulk velocity of each phase must adjust (considering density differences) 

so that the pressure drop along, say, a streamline in the center of the duct 
--- which lies entirely in the vapor phase--- is equal, at any z, to the 
pressure drop along a streamline adjacent to the wall, lying entirely in 

the liquid phase. Thus, we now have two separate momentum equations to 
deal with, and therefore two different phase velocities, which precludes 
treating the two-phase flow as a simple "homogeneous' fluid. (At the inlet, 
assuming quality is exactly zero, the first infinitesimal quantity of vapor 


generated causes no finite change in flow conditions. Therefore, at the 
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inlet, and only at the inlet, the two phase velocities are equal. ) 

The concept of "quality'’ becomes somewhat more complex, given 
that liquid and vapor velocities are generally unequal. Actually, quality 
is a static concept, which was easily extended to the preceding dynamic 
case, because the phases moved together. To avoid confusion, the term 


will here be taken to mean: 
Sz S. 
es 2. ane ee a (44) 
+ + 
P.,s.dz Ps, dz Ps, Ps, 
where we again speak in terms of unit width, and S. is that portion of the 


slot height occupied by the vapor phase ise + s, = s, as defined above). 


l 


For convenience, the variable é will be defined as: 


A S 1-x 
l i joe 
g AL S., x p, (45) 


In terms of € , the continuity equation, written between inlet and outlet, 
now becomes: 
v,(1 + )= Vv + Vv (46) 
ri, : § Py, 14$ Py, es 
Because specific volume is a known function of pressure, the two 


momentum equations can be integrated directly between inlet and outlet: 








1 
1 2 2 
Aon = (v," - v,") (47) 
1 7l4/0)) te 3 
Ap Ss (v ay 2) (48) 
where | ie feos 
p Pe “Pj 


which can be evaluated by direct integration of the appropriate saturation 


property polynominals. Of course, as was stated above, we require 


Ap = Ap (49) 


1 Vv 


In principle, eqns. (45-49) can be reduced to three equations with the 


three unknowns Vig Weans and could therefore be solved explicitly for the 
O 
ds 





inlet and outlet velocities. In practice, this is messy, and an iterative 
approach was employed here. This problem was solved in tandem with 
the no-slip problem, preceding. The inlet velocity determined for the 
former problem was inserted in eqns. (47) and (48), which were then 
solved for ae and v5" These values were inserted in eqn. (46), which 
was solved for a corrected ''guess' for v,, etc. Convergence was found 
to be quite rapid: generally, the value of Vio converged in a Single iter- 
ation, and vie and va in two or three. 

Siecourse, once determined, Tie Ce be checked against Vim? and 
assuming vy < Yim? We again proceed. With inlet and outlet variables 
completely determined, it is again convenient to adopt an iterative, numer- 
ical solution. In this case, the energy relationship becomes: 

(q/A) dz = hy, Ow... (50) 
where Ow, is the change in the vapor mass flow rate while traversing the 
distance dz. For an incremental analysis, let 

K = GAA) (51) 

S 

Again, s is unknown, but if the two problems are solved in tandem, we can 
use the slot height determined from the no-slip model for the same set of 
thermodynamic parameters. Using an approach similar to that of the 


h 


preceding problem, at the Ke section: 


Wi Se = Ww /s+K/h 
ae Vey k-] 


Specifically, at the inlet, ie is defined to be zero. Further, we + a = Ww, 


which is still constant. To advance the solution from one station to the 


! = t= ! t= 
next, let st s_/s and s s/s (s_ +8, 1). Then, 


1 
Wie fl = au )(w /s +K/h )=v_ (l-s,') (52) 
ae aie vie Ivey Ye 

v, s,' = (1/ )(w /s - K/h ) (53) 
L. 1, as lend Ivy. 4 
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- and, to keep Ap, = Ap, between stations: 
wey de OQ We eve) (54) 
k-1 ‘k k-1 1-1 k k-1 


1 


can be used to advance the solution consecutively from one section to the 


At the inlet, Sone 0, s. = 1, Wit ), and w,/8 = P,,%:: Eqns. (52-54) | 


next, although an iterative solution is again convenient. Here, eqn. (53) 
was solved for aie under the assumption that 
Vv, -Vv 
= pO 
vy FV —————— 
k k-1 WN - (k-1) 
eqn. (52) was then solved for Vy after which eqn. (54) was used to derive 
a correction for the assumed value of ie 
At the outlet, computed and desired outlet qualities should be compared, 
and the process repeated with a new value of K, until desired convergence 
is obtained. As was the case previously (when eqn. (42) was used to com- 


pute gs), eqn. (51) can be solved for that variable, after the solution has 


converged. 


Utilization of the above models: actually, both these problems were pro- 
grammed for computer solution, but in terms of the parameter S, which 
is defined by the relationship: 
S = s/(q/A)L (55) 

The reason for this shift of variables is that the slot factor, S, is the more 
convenient parameter for the contemplated application: for a specified P, 
and Xs a single curve can be developed, giving S as a function of Mp. A 
family of such curves, for various Xo» concisely displays the general 
relationship between inlet and outlet pressure, quality, and slot height, 
for any (q/A) and L, 

Further, while the incremental models give access to a temperature 


profile, through application of the saturation temperature polynomial to the 
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computed pressure at each station, a more useful figure for preliminary 
analysis would be the mean coolant temperature, or better still, the mean 


temperature increment. The former is defined as: 


N 
1 


UOT aa ails 

wttfrezth T. (56) 
k=1 

The latter might be defined as AST. = D), = Ls but as a matter of conven- 


ience, will be defined instead as: 
AT. ; 12 - Fey I seiiiage (57) 

under the assumption that the one-atmosphere saturation temperature may 

be more readily available than the outlet temperature, if outlet pressure 

differs from one atmosphere. 

The primary advantage of dealing with the temperature increment, 
defined by eqn. (57), rather than the mean temperature, is that the mag- 
nitude of the increment clearly underlines the added ''cost'' of operating 
at elevated pressures. In computing mean coil resistivity, the first esti- 
mate for mean coil temperature would be given by: 

te! SS ot 


5,1 + ra + AT, (58) 


Mere, oe is actually the wall temperature of the cooling duct, and AT, 
is the temperature increment across the thermal boundary layer, which 
would presumably be determined from the same source as the figure used 
for mean (q/A): either a predictive correlation, as discussed in the pre- 
ceding chapter, or experimental data. This is classed as a ‘first estimate’, 
since if the coolant ducts are widely spaced, an additional correction might 
be necessary to account for thermal gradients in the conductor material. 
However, this factor is again best left for the detailed design phase. In 


general, the T defined in eqn. (58) would be the temperature used to enter 


the mean coil resistivity curves of Appendix B. 
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Correction for subcooling at inlet: Either of the above models could be 
applied to the problem at hand--- e. g., to predict coolant volume fraction 
for a cryogenic coil--- depending upon which was most representative of 
flow conditions realized in the type of slot geometry contemplated. Acaka, 
it is to be regretted that experimental data are simply not available for 
forced flow, two-phase cryogenic cooling with this type of slot geometry. 
Presumably, all efforts would be made to maintain operation in the ‘'annu- 
lar flow'' regime, in which case the no-mixing model would seem most 
pertinent. However, in the absence of experimental data, the no-slip 
model has the advantage of producing a more conservative prediction: 
Since the relatively dense liquid must accelerate with the vapor phase in 
that case, flow velocities tend to be lower for a given Ap. This, in turn, 
requires larger slots to pass the same bulk rate of flow, as would be 
required to maintain the same specified outlet quality. Until experimental 
data confirm that viscous effects are truly negligible, throughout the range 
considered here, this ''built-in conservatism’ might well be used as ‘'ad- 
vance compensation’ for possible viscous effects. 

In any event, before either model can be applied to predict flow rela- 


"real" coolant ducts, one additional correction should be intro- 


tionships in 
duced. Both models assumed zero quality at the slot inlet, but as was 
mentioned above, the zero-quality saturation point would not, in fact, be 
reached until the subcooled liquid coolant, coming from the pump(s), had 
been heated to its saturation temperature. Particularly at high pressures, 
where hb, tends to be relatively small, the energy absorbed before reaching 
the saturation point might be significant. (Note, however, that it should 
not be predominant, or this could not properly be classed as ''two-phase 


cooling’. However, the pressure ranges selected for the various coolants 
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in effect constrain the problem to avoid that possibility. ) 

What the above models actually predict is the relationship between 
pressure, slot factor, etc., and the length of that portion of the slot 
devoted to two-phase flow. Suppose zero quality is reached at the point 
Zo» rather than at the inlet. The liquid at the inlet would presumably be 


subcooled by an amount On, where 
Oh =h 


Ideal’ pumps are isentropic, but ''real'' pumps are not. Assuming that 


h (59) 


: 5 yn 
sat, P; sat, Pe pump 


the efficiency of a ''real'' pump would be of the order of 50%, we have: 
2(p.-p_) 
te 


Oh ane : P1, (60) 


Then, assuming constant (q/A), as before, 

Oh(0,v,8) = (q/A)z : 
Let the distance (L-z.) be designated L', and the slot factor computed with 
L' be S'. Then, 

t = = = t 

z/L'= HhiO.v,) TaTAYL = Ohf,¥,S 
Actually, the slot factor derived from the previous models is, correctly, 
S'. If S' were multiplied by L, vice L', too large a value would be com- 
puted for s, Presuming L would be known, while L' might not be, the 
previously computed slot factor must be corrected: 


L'S' L'S! S! 
2 See See (61) 
IS L' + Zo 1 + 6hO,v,5 


Eqns. (59-61) can be applied to correct the slot factor predicted from 
either flow model. 

Next, the mean temperature increment predicted from the preceding 
models must be re-examined, The original ae which will now be called 
Ty was based on the nonlinear, monotonically decreasing temperature 
profile between z = Zo and z= L. By contrast, as discussed previously, 
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conditions in the subcooled flow can be approximated as a linear temper- 
ature increase between z=0 and z=z.. (Note that the "hot spot’ in the 
windings of a coil cooled by such a system would occur just at Z=Z 0 where 
coolant temperature would roughly equal the coolant saturation temperature, 
at the inlet pressure.) Since a pump efficiency of 50% would imply rela- 
tively little temperature increase for the coolants considered, even when 
pumped through the maximum pressure increment investigated, the cor- 
rected AT . can be approximated as: 
eae 4s ra (at 1° Te )/2 
Nae = CG O 6S, Pj S,Po _- 
° iu 
Tm + ObPiviSiT. ,. ; Teele a ae 
1 + 6hQ,v.S' Bee 


Finally, it was originally stated that inlet pressure will be system- 


s,l 


atically incremented" by assuming P, = Py +nAp. Again, the pressures so 


specified correspond to conditions at z=z,., vice z=0. With single-phase 


o? 
flow to Z=Z os Avp could be corrected by deriving an additive term from eqn. 
(30). However, in general Zz << IL, and we choose to remain with an 
inviscid model. One correction can be applied, though, even to an inviscid 
model: The liquid density decreases between z=0 and Z=Z os Since liquid 
density is a significant function of temperature, if not pressure. This 


requires a slight flow acceleration, which results in a pressure decrement 


of the order: 


~ 2 
AP niet ~ Pi¥;) CPs, p. : Ps, p) (63) 


thus, a ''corrected'’ Ap can be determined by adding eqn. (63) tothe Ap 
specified in the incremental analyses. This has been done, though in most 
cases, the correction is insignificant. 

The details of the calculations performed, in applying the no-slip and 


no-mixing models, corrected for subcooling at the inlet, to the analysis of. 
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the two-phase cryogenic cooling problem, are discussed in the next chapter. 


The results of these calculations are presented in Appendix D. 


E. Other coolant flow considerations: 

In closing this chapter, two factors which cannot be discussed quan 
tatively deserve qualitative mention: 

The first is the problem of inlet and outlet connections to the coolant 
ducts: although the slot geometry discussed herein probably corresponds 
reasonably well to such geometries as might be employed throughout the 
greater portion of ''real'' cryogenic saddle coils, the applicability obviously 
terminates at coolant inlet and outlet connections. In fact, no simple 
solution is seen to the problem of compatability with coil geometry at 


spiraled" through several 


these points. Presumably, ducts would not be 
axial passes, to emerge at the surface of the coil, since the analysis 
indicates that increased duct length requires increased slot height, and 
therefore decreased coil packing factor. Further, the geometry of the 

‘end turns’ of the coil is not particularly compatible with conventional 

use of inlet and outlet headers. The problem is a moot one, and regretably 
beyond the scope of this study. 

The second problem has to do with the possibility of a parallel flow 
instability in two-phase cooling systems. It is probable that whatever 
connection scheme were devised would feed several ducts ‘in parallel’. 

Let us assume that steady-state operation has been established, and investi- 
gate what might happen if some local phenomenon perturbs the flow rate 
through a particular duct. 

Fig. 21a depicts two curves of flow rate vs. Ap, such as might be 
characteristic for a particular coolant and a particular geometry at con- 


stant (q/A). Basically, characteristics of class I are "safe", but those 
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of class II may lead to problems: 

If a high-capacity pump is supplying several ducts in parallel, the 
pump characteristic might imply essentially constant Ap, regardless of 
flow conditions in any one duct. Such a characteristic is shown super- 
imposed on a class II characteristic in Fig. 21b. If steady-state operation 
had been established at point a, and a sudden fluctuation caused a momen- 
tary decrease in mass flow rate ina particular duct, the operating point 
for that duct might ''jump' to point b. If the quality at point b were so 
high that flow conditions entered the ‘liquid deficient" regime, heat trans- 
fer to that duct would then drop abruptly, leading to the creation of a “hot 
spot'' extending the axial length of the coil. The problem is obvious; the 
situation can be avoided or ameliorated by careful design; but that is 


also beyond the scope of this study. 
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VII. GENERATION AND USE OF CHARACTERISTIC CURVES 


A. Curves of coil geometry: 

Generation: The curves of coil geometry, presented in Appendix A, and 
those of mean coil resistivity, in Appendix B, were generated from a single 
computer program, consisting of a calling sequence (main program), and 
several functional subroutines. A brief resume will be given here of the 
logic used in these routines, and of the numerical ranges investigated dur- 
ing the processing of data. The programs themselves are presented in 
Appendix E, All are written in Fortran, for use with the M.I. T. Compu- 
tation Center version of the Fortran Moniter System, and were run on the 
Computation Center's I. B. M. 7090 digital computer. 

MAIN, as mentioned above, is primarily a calling sequence. As per- 
tains to the coil geometry problem, MAIN reads an input data card, to 
determine the desired channel area (in square meters), desired channel 
flux density (in kilogauss), and effective coil current density (in amps/ 
Square meter). The calling sequence is ''looped", so that data from one 
card are processed completely, then another card read, etc. As processed 
here, channel area was varied through the specific values of 10, 1, 0.1, and 
0.01 oo flux density through 25, 50, 75, and 100 kilogauss, and current 
density varied as necessary to investigate form factors in the approximate 


s) 6 


range 1.3 = F=5(e.g., 10° to 10! alm? EG A. = 10, andd x10 to 


h 


toe a/m? for A. = 0.01). Customarily, B was varied through all its val- 


h 


ues, for a fixed Je and A then Je incremented, and the process repeated, 


clas 
until all coils suitable for a particular channel size had been investigated. 
SIZE performs the actual calculation of coil geometry for a specified 


set of input parameters. Eqns. (1-6) have been programmed in the ration- 
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alized MKS system of units; thus, the specified B is divided by ten to con- 
vert to a consistent set of units. SIZE uses an iterative approach to deter- 
mine the required cylinder radius for a particular coil: eqn. (1) is solved 
for cylinder offset, and a ''first guess'' is made that 

R=at Bete. | 2 
This value is inserted in eqns. (2) and (3), and the resulting computed chan- 
nel area compared to the specified channel area. Based on this comparison, 
the estimated R is varied systematically, until the deviation between spec- 


ified and computed A. is less than 0, 01%, or until the process has been 


h 
repeated 50 times. (The latter option was found to be an unnecessary pre- 
caution, within the range investigated. ) The subroutine then computes the 
width, height, and area of the optimum" channel for the resultant coil, the 
coil length, total cross section (both lobes), total surface area, and total 
volume (uncorrected for end turns), and finally, the coil form factor. This 
information is stored until MAIN calls the output routine, PRINT. 

Use of coil geometry curves: The output of the above subroutine is 
presented in Figs. A-1 to A-3, as an explicit function of the form factor, 
F, defined in eqn. (6). F has been chosen as the independent variable, for 


" Shape" 


two reasons: It provides a convenient, quantitative index to the 
of a saddle coil, and more importantly, provides an immediate indication 
of design optimization. (This latter point will be discussed in Chapter 
VIII. ) To utilize the information presented in these curves: 


l. Specify desired A B, and F., 


Chis 
2, Enter Fig. A-1 with F, and determine J required to attain that F. 


3, Enter Fig. A-2 with F. Compute w 


plotted values by 10 \/ A Gh ; 


4. Enter Fig. A-3 with F. Compute A.and Aa by multiplying plotted 
3/2 


che 2 and R by multiplying 


values by 100 A on Compute V by multiplying plotted value by 1000(A 
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5. Correct for end turns by multiplying A.and V by F. 

As stated in Chapter I, the above assumes that the coil length, L, is 
equal to 10VA_. However, the information plotted in these curves can be 
corrected for L # LOA 

l. Carry out steps 1-4, as above. 

Bee Correct A.and V by multiplying by L/10 Aun 

3. Correct F: F'=1+(F - 1)(10VA,,/L). 

4, Correct for end turns by multiplying A,and VoDy te 

Note that the optimization procedure discussed in the next chapter per- 
tains to F', vice F, when L is not equal to 10VA.. However, all parameters 


plotted as functions of F must be "looked up" as functions of F. Of the geo- 


metric parameters, only A.and V need be corrected for the ‘new’ length. 


B. Mean coil resistivity curves: 

Generation: Mean coil resistivity was computed as a function of tem- 
perature, at ten degree increments, from 5° to 95° K, for each coil proc- 
essed by MAIN. To reduce input-output time, the necessary field-free and 
magnetoresistance data for a particular metal were ‘'read in’ only once, at 
the start of each run, by having MAIN call READ before entering the com- 
puting loop. 

READ is strictly an input routine, designed to read, in sequence: 
l, Ten specified temperatures, from a single card; 
2. The field-free resistance corresponding to each specified temperature 
(see Tables 1-4), from two data cards; 
3. The magnetoresistance for the particular metal, as tabulated in Table 
o, from ten data cards. 

FIELD determines the field distribution within the upper-right-hand 
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quadrant of the coil, to permit calculation of magnetoresistance. (A modified 
version of FIELD was also used to compute the vector diagrams presented 

in Figs. 12-16.) The upper-right-hand quadrant was selected, since the 
symmetry of the field distribution indicates that the mean resistivity of any 
quadrant is equal to the mean resistivity of the entire coil. 

FIELD positions a "grid" of 5000 squares over the upper-right-hand 
quadrant of the coil. The overall size of the grid is 2R wide by R high, with 
the lower left corner of the grid positioned at the center of the MHD channel. 
Each component square has sides of dimension R/50. 

FIELD first scans the left-most vertical column of squares, working 
from bottom to top, and determines which squares have centers which fall 
within the coil quadrant. FIELD computes and saves the value of the field 


intensity at the center of these squares, using the equations: 





2 

J R 
B. = Mo" e 1 - a y (64) 

2 al) eae 

ein R? 
B = (a+x) | ———-—_——— + (a-x) (65) 
y Z (xta)? + y" 
IB| =\/B7+B7? , (66) 
x ha 


where the x-y origin is the lower left corner of the grid. After locating 

the uppermost square which falls within the coil quadrant, FIELD shifts to 
the next vertical column to the right, etc., until locating the first such 
column which lies entirely outside the coil. The field intensity at the center 
of each square which fell within the coil is saved in an array, with the loca- 
tion of the data corresponding to the location of the square in the grid. Loca- 
tions corresponding to squares which fell outside the coil are set to zero. 
FIELD also saves the location of the first column which fell entirely out- 

side the coil. 
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RHO next computes mean coil resistivity, using the field distribution 
generated by FIELD, and the resistance data supplied by READ. RHO 
sequentially examines each specified temperature. For each temperature, 
RHO determines the field-free resistivity, then scans the array of field 
intensities for non-zero entries (examining only those columns which were 
previously scanned by FIELD). RHO divides each non-zero entry by the 
field-free resistivity, and used the result as an index to find the "nearest" 
corresponding entry in the magnetoresistance table. RHO then computes 
"total resistivity at the center of the corresponding square, and assumes 
this to be representative of the total resistivity of the entire Square. Fi- 


nally, RHO computes mean coil resistivity, for the particular temperature, 


N 
i » Pp (67) 
i=l 


N being the number of locations found to have non-zero entries. 


from the relation: 


In this incremental analysis (as with all other incremental analyses 
conducted herein!), the fineness" of the grid was selected through suc- 
cessive reduction of the increment used, so that the degree of approximation 
associated with a particular increment size could be estimated. The incre- 
ment specified ensures that the accuracy of the computational technique is 
consistent with the accuracy of the physical model. 

Use of resistivity curves: The datain Tables 1-5 indicate that indium 
coils would not be competitive with coils wound from the other conductor 
materials investigated, except at temperatures below 4° K, _As it seems 
unlikely that cryogenic coils would be designed for operation below 4° Kk, 
indium was dropped from the list of potential conductor materials, and no 
mean coil resistivity curves are presented for that metal. However, curves 


of mean resistivity vs. temperature for coils wound of sodium, copper, and 
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aluminum, are presented in Figs. B-1 to B-16. These curves have been 
normalized to the 1.7715 micro-ohm-centimeter standard, to permit direct 
comparison with ''conventional' copper coils. Further, resistivity curves 
for all three metals are superimposed in each figure, to permit direct com- 
parison between these metals as cryogenic conductors. Figs. B-1 to B-4 
are for 25 kg coils, B-5 to B-8 for 50 kg coils, etc., with each figure pre- 
senting data for a particular form factor. (F and B together determine the 
geometric distribution of magnetoresistance, and therefore define mean 
resistivity, for a particular metal. However, the variation is sufficiently 
slow--- except for copper below 30° K--- that interpolation is generally 
unnecessary: to determine mean resistivity, simply enter that figure 
which most nearly corresponds to the contemplated application. 

While these curves were computed from the data presented in Tables 
1-5, approximate adjustment can be made, should it be desired to assume 
some value for any one of the three components of ''total'"’ resistivity, 
other than those used here: the plotted ''mean" resistivity can be compared 
to the tabulated ''field-free'' resistivity, to determine the percentage contri - 
bution of each component to total resistance. (Note that the 0° K entry in 
each of Tables 1-4 is the assumed residual resistance ratio for the various 
metals.) Once the percentage contribution of a component is known, the 


net effect of changing that component can readily be estimated. 


C. Structural curves: 

Generation: The curves of structural volume fraction, presented in 
Appendix C, were generated by a second computer program, consisting of 
a calling sequence (main program), and the subroutines SIZE, CSHELL, 


and TSHELL. 
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MAIN is again "looped"; uses the same input data as the previous pro- 


gram (A B, and J,), and again calls SIZE to evaluate the geometric 


ch’ 
parameters of each saddle coil so specified. MAIN then systematically 
varies the specified conductor stress limit (O.) through the values 3, 000, 
6,000, 9,000, and 12,000 p.s.i., and the structural stress limit ( oF) 
through the values 120,000, 90,000, 60,000 and 30,000 p.s.i. (All values 

of oF are ‘matched’ with a singie value of oF before the latter is incre- 
mented.) CSHELL and TSHELL are called for each combination of 0. and 
OL , to compute the compressive and tensile volume fractions, respectively. 

To explain this particular sequencing: The required number of struc- 
tural shells decreases as 0. increases, Thus, if no shells are required 
for a particular oe. none will be required for higher on Accordingly, 
if total computed volume fraction is found to be identically zero for any 
value of oy MAIN does not search to higner values, but concludes that 
the problem has been ''solved", for that particular coil, and calls for new in- 
put data, 

Similarly, shell thickness increases with decreasing one Thus, if the 
structural volume fraction 1s unreasonably high, for a given On MAIN does 
not search to lower O » but increments on immediately, to reduce the 
number of shells. ‘Unreasonably high'' is defined, in the program, asa 
volume fraction of 0.8 or greater. 

SIZE, aside from a minor reduction in the number of "output’ parameters 
computed, is the same subroutine as described above. 

CSHELL computes the compressive volume fraction by the procedure 
described in Section D of Chapter III. Since electrical variables and coil 
geometry data are expressed in metric units, while Windenburg's equations 
(and most available strength data) are in English units, considerable con- 


version of units occurs during the computational process. 
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TSHELL computes the tensile volume fraction by the process described 
in Section E of Chapter III. Specifically, TSHELL slices the coil lobe into 
shells of thickness dr = a/10, where a is the offset between the center of the 
cylinder and the center of the channel. Per unit length, the volume of each 
such shell is assumed to be 2rQ0 dr, where r is the inner radius of the shell, 
and © is defined by eqn. (20), in terms of that radius. TSHELL integrates 
eqn. (18) numerically, by ''piling up" such shells until the effective pressure 
exceeds Oo. TSHELL then "shaves" the last shell, reducing its thickness 
so that Beet 0. at the outer radius of the stack. A structural support shell 
is inserted at this radius, and its volume determined as indicated in eqn. 
(22). TSHELL continues this process until the outer edge of the coil is 
reached, then divides the accumulated volume by the volume of a unit length 


of coil lobe. 
When the above computation has been completed, TSHELL prints Ce 


Oo , the number of compressive and tensile shells, the compressive and 
tensile volume fractions, and the total structural volume fraction. Control 
is then returned to MAIN, which specifies a new combination of 0. and 
Ow or begins work on the next coil, as appropriate. 

Use of structural curves: Curves of structural volume fraction as a 
function of form factor are presented in Figs. C-1 to C-4, with each figure 
presenting the data for a particular OF Interestingly enough, for a given 
OF structural volume fraction was found to be much more sensitive to F 
and B, than to 0. Accordingly, curves are presented for each value of 
B investigated, but only for of = 3,000 and 12,000 p.s.i. The structural 
volume fraction for intermediated values of 0. can be estimated, to within 
the probable accuracy of the model, by simple linear interpolation between 


these two families of curves. 


D. Coolant curves: 


teal 





Generation: The coolant curves, presented in Appendix D, were gen- 
erated in two stages: First, curve fitting routines were utilized to develope 
the saturation property polynomials, and then, using these polynomials, 
curves of slot factor (S), mean coolant temperature increment (AT), lately 
coolant inlet velocity Fee, were generated by the techniques discussed in 
Chapter VI. 

The curve fitting routines were executive programs built around an 
M. 1. T. Computation Center SHARE program, GLSQ. Briefly, GLSQ deter- 
minesa least-Square solution to an overspecified problem (M equations in 
N unknowns, with M > N), or the exact solution to a ''properly set'’ problem 
(M equations in M unknowns), by matrix inversion techniques. 

Both GLSQ, and an example of the type of data fitting routine used in 
conjunction with it, are presented in Appendix E. The data fitting routine 
specifies the functional form to which the data are to be reduced--- e.g., 
Do) Ay qe ee Sp ax ta= and for each data point (x, y), computes 
the various powers of x. (Actually, the functional form need not be con- 
strained to a simple power series in x: it could as well involve terms in 
oo Sin x, etc.) These factors are sent to GLSQ in an array, with each 
row of the array corresponding to one data point. (Obviously, there must 
be at least as many data points as there are coefficients to be evaluated in 
the set {a,} .) GLSQ solves the matrix for the unknown coefficients {a.} : 
The example program prints these coefficients, the resultant minimum 
mean square deviation, and lists x, y 


observed’ ’generated’ Yobs gen’ 


and 100(y Ue (% deviation) for each data point. (Input required: 


obs "gen 
one card, listing the number of data points (M), and M cards, each listing 


the coordinates of a single data point.) Of the two figures which appear 


repeatedly in the DIMENSION statement, 200 is the maximum number of 
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data points which can be processed, and 7 is the number of terms in one row 
of the matrix. Both of these figures may be adjusted, as desired. 

The coolant flow program consists of a main program, the subroutines 
READS, NOSLIP, NOMIX, INLET, and PRINTS, and the functions SVOLL, 
SVOLV, TEMP, HLIQ, HTVP, AVSVL, and AVSVV: 

READS determines the saturation property coefficients (as listed in 
Table 6) for the coolant to be investigated, the three values of outlet quality 
to be examined, and the limiting velocity to be applied to the flow. 

NOSLIP reads in the outlet pressure to be assumed, the ''step size" 
to be utilized in incrementing the inlet pressure, and the number of incre- 
ments to be examined. NOSLIP then solves the no-slip flow model, exactly 
as discussed in the preceding chapter, using a 25-station incremental model. 
Iterative solutions are cycled until the one-cycle change in the control vari- 
able is less than 0.1%. Fora given P; and Po» flow conditions are solved 
for each specified Xo: After this, P; is incremented, and the process re- 
peated until all specified inlet pressures have been investigated. 

NOMIX solves the no-mixing flow model in the same manner. One 
additional point should be noted: Statements 9 and 10 of this subroutine 
generate a corrected estimate for the constant K, which appears in eqns. 
(51-53), in order to obtain convergence between computed and desired out- 
let qualities. Actually, the "loop gain" in this particular iterative solution 
is a strong function of the degree to which specified operating pressures 
approach the critical pressure of the coolant being investigated. The atten- 
uation factors of 10.0 and 5.0, which appear in statements 9 and 10, as 
listed in Appendix E, provide suitable convergence for nearly all runs 
covered here. However, if specified pressures lie quite close to the critical 


pressure--- e.g., above about 7.5 atmospheres for liquid hydrogen, or 
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Ohe atmosphere and above for liquid helium--- loop gain is so small that 
the solution converges too slowly if these factors are used. They should 
be replaced by the factors 3.0 and 1.5, respectively. (Conversely, if the 
latter values were used for, say, liquid nitrogen in the pressure range of 
1 to 10 atmospheres, the resultant correction would be excessive, and the 
solution would oscillate, or diverge. ) 

INLET corrects the flow solution for subcooling at the inlet, as indicated 
in eqns. (61-63). 

PRINTS can be used to list the output data from all three of the preceding 
analytical subroutines, or, by omission of statements 1-25, can be used to 
obtain only the ''corrected' data, generated by INLET. The latter option is 
obviously less time consuming. 

The functions SVOLL and SVOLYV utilize the saturation property coef- 
ficients, obtained from READS, to compute, respectively, the specific vol- 
ume of the saturated liquid, and specific volume of the saturated vapor, in 
accord with the functional forms listed in Section B of Chapter VI. Similarly, 
TEMP computes saturation temperature, HLIQ the specific enthalpy of the 
saturated liquid, and HTVP the heat of vaporization, for a specified pres- 
sure, AVSVL and AVSVV respectively compute the mean specific volume 
of liquid and vapor, for use in eqns. (47) and (48), for a specified P. and Po 

Use of coolant curves: Curves of 5, AT, and Vin VS: Ap for each 
of the various coolants studied are presented in Figs. D-la to D-l2c: D-1 
to D-3 pertaining to liquid nitrogen, D-4 to D-6 to neon, D-7 to D-9 to 
parahydrogen, and D-10 to D-12 to helium. 

Each figure presents a family of curves, showing variation of a par- 
ticular property with Ap, for various outlet pressures, at a specified 
quality: e.g., D-la presents data for 30%, D-1b for a 90%, and D-lc 


for aS 70%. The ''a-b-c'' sequence is similarly repeated in all other fig- : 
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ures. Similarly, the D-1 set are curves of S vs. Ap, D-2, AT. VS. 2. ps 
and D-3, Vin VS Ap. This sequence is also repeated for each coolant. 
Finally, the data generated by the no-slip and no-mixing flow models are 
displayed side-by-side, in each figure, for ease of comparison. 

In general, the plotted curves should not be extrapolated to higher values 
of Ap. Each curve is terminated, either because vee would exceed the 
specified limiting velocity, or because Pin would exceed the specified limit- 
ing pressure, were Ap increased further. However, the curves could 
certainly be extrapolated to lower values of Ap than those plotted. 

For a specified Py and Ap, the inlet velocity curves permit estimation 
of the required coolant pump power, as indicated by eqn. (34). Similarly, 
the AT. curves permit estimation of that parameter, for insertion into 
eqn. (58), and the slot factor curves permit estimation of the coolant volume 
fraction: given a preliminary estimate of q (total I7R loss), q/A, and L, 

coolant volume = V, = As = fa/(a/A)] [S(q/A)L] = qSL (68) 
and 
Cat os ae (69) 
where V is as given by eqn. (7). 

Thus, these curves supply the final link in the chain of information 

required for preliminary analysis of cryogenic saddle coils. A suggested 


procedure for conducting such an analysis is discussed in the next chapter. 
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VIII. DESIGN OF CRYOGENIC SADDLE COILS 


A. General design characteristics: 

Presumably, the first parameters which would be specified in the design 
of an MHD saddle coil would be the desired channel size, and desired field 
intensity: these requirements are set by the intended MHD application, and 
could not conveniently be altered to suit the whim of the coil designer. Since 
we have assumed L = 10 5 ae in line with Sutton's work (28, 29), this implies 


specification of B and A as the first step in the design process. 


Clan 

Assuming this has been done, the next step in the design of cryogenic 
saddle coils would be the selection. of a particular metal for the coil con- 
ductors, and the selection of a particular coolant (and, by implication, spec- 
ification of the approximate range of possible operating temperatures). ‘Next’, 
because the selection is obviously limited; the choice can be simply and di- 
rectlyrelatedto the overall design philosophy, and will have a profound 
influence on subsequent development of the design. 

Figure 23 depicts the total electric and refrigeration power requirement 
(neglecting heat leak and coolant pump power) as a function of coil operating 


2 


temperature, for a ‘'typical'' cryogenic saddle coil (B= 50 kg, A, = 0.1m’, 


fed at 
F = 2), assuming unity packing factor. The approximate range of operating 
temperatures attainable, with each of the liquid coolants examined here, is 
indicated at the bottom of the figure. Although the specified assumptions 
indicate these curves to be a ''rough'' approximation to the characteristics 
of 'real'' coils, several conclusions nevertheless appear valid: 

l. Under the same assumptions, a conventional, copper saddle coil 


would dissipate 16.7 Mw. Comparing this to the curves in Fig. 23, itis 


apparent that liquid nitrogen would not be used to cool "working" saddle coils. 
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2. If liquid nitrogen were used to cool an experimental coil, copper 
conductors would yield lowest power consumption. If coil weight were an 
important factor, aluminum conductors would seem the best compromise 
between weight and power consumption. 

3. In the liquid neon range and below, aluminum appears to be the 
"optimum" conductor material. Optimum operating temperature for aluminum 
saddle coils is about Re K--- i.e., about the same temperature recommended 
by Taylor and Post (26) for aluminum solenoids--- and the optimum is fairly 
Sharp: a i(0h temperature shift in either direction increases losses by more 
than 50%, 

4, Based on this crude model, operation at liquid neon or liquid helium 
temperatures would appear roughly equivalent, from the viewpoint of power 
consumption. However, a comparison of the slot factors presented in Figs. 
D-4 and D-10 indicates that use of the latter coolant would result in lower 
packing factors, and thus greater losses. Again, liquid helium probably 
would not be used as a coolant for cryogenic saddle coils, except when 
cryogenic conductors are used to ''stabilize'' superconductors, as discussed 
by Kantrowitz and Stekly (14). 

o. Ifa two-phase coolant is to be used, liquid hydrogen would seem the 
best choice. 

6. Supercritical helium is the only coolant which could maintain the 
‘optimum temperature for aluminum coils. 

Thus, assuming the motivation for utilizing cryogenic saddle coils is 
that of reduced operating power requirements, one would normally select 
aluminum for the coil conductors, and liquid hydrogen or supercritical helium 
for the coolant, based on these data. 

The next fundamental variable to be considered, after the specification 
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of B, A.,, the conductor metal, and the coolant, is the general ''shape"’ of 
the contemplated coil. While the specified saddle coil model seemingly 
determines this, a glance at Figs. 12-16 will indicate the wide variation of 
form still possible within this family of coils. Of course, ‘'shape'' can be 
quantitatively described by the value of the form factor, F, and Fig. 24 
displays relative power consumption as a function of F (under the same 
assumptions that applied to Fig. 23). Here, the minimum occurs at F = 
2.0. Power consumption rises gradually as F is increased above this value: 
loss density (Oe) varies inversely with F, but a comparison of Figs. A-1 
and A-3 indicates that Je decreases relatively slowly with F, above F = 2, 
while coil volume increases quite rapidly. Asa result, total power con- 
sumption and coil size both increase if F increases above 2.0. One can 
therefore immediately class form factors above that value as "not of interest". 
Conversely, total power consumption rises if F is decreased below 2, be- 
cause the sharp rise in Ji more than offsets the decrease in total volume. 
However, though Fig. 24 indicates that a coil with F = 1.34 would have twice 


the power consumption of one with F = 2, Fig. A-3 indicates that its volume 


would be smaller by more than an order of magnitude. Thus, for some appli- 


cations, form factors less than 2 might be desirable. 

As defined here, the minimum attainable F is theoretically 1.22: this 
corresponds to a sinusoidally distributed current sheet (infinite current den- 
sity), flowing axially in a cylindrical shell of zero thickness. Realistically, 
Figs. C-1 and C-4 indicate that structural requirements alone would con- 


strain ''real'’ coils to form factors above about 1. 3-1.6, depending upon the 


specified flux density, and allowable stress levels. Thus, we can define the 


‘range of interest'' of coil shapes at those corresponding to form factors 


between 1. 3 and 2.0, with the lower end of the attainable range corresponding 
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to minimum power consumption. 


B. Simplified design procedure: 


Bearing the above in mind, Fig. 25 presents a flow diagram of the 
recommended design procedure for cryogenic saddle coils. Although the 
motivation for much of the sequencing is obvious, some explanation of var- 
ous points might be useful: 

The design cycle begins with F=2, since both coolant and structural 
volume fractions are lowest at this point, within the range of forms defined 
above. The design process entails an iterative search on these factors, to 
estimate maximum attainable packing factor for a given F. Packing factor, 
of course, is determined from the relationship: 

A =1-(S. VF. +C.V.F.). (70) 
As F is decreased, the point will be reached where the sum of the coolant 
and structural volume fractions equals or exceeds unity: i.e., where no 
physical solution exists. However, if a physical solution exists at all--- 
that is, if it is possible to build any ''real" coil, for a given B, A uy con- 
ductor, and coolant combination--- the form factor of 2.0 will correspond 
toa 'realizable'' coil. Successive reduction of F can then be used to deter- 
mine the smallest ''realizable'’ coil, should that be desired. 

Once B and F have been specified, the structural volume fraction can 
be determined immediately (from Appendix C), since it is not a function of 
resistivity, temperature, coolant pressure, etc. However, to determine 
the required coolant volume fraction, total power dissipation must be known, 
and to determine power dissipation, packing factor must be known. For non- 
unity packing factors, loss density becomes Oak = OU IA i and total con- 
ductor volume becomes )FV. Thus, 
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PAAY= PVFIIA = PA =A (71) 

It is suggested that a "first guess" be made that J = 0.75 - S.V.F. For 
the smaller F=2 coils, this is reasonably accurate, and for the larger F=2 
coils, it will be found to be conservative. ie can then be estimated, the 
corresponding coolant volume fraction computed, and Xr recomputed from 
eqn. (70). Note that if the ''recomputed" ,) equals or exceeds the "assumed" 
A, the "assumed" Xr is physically realizable. However, unless the two 
are equal, the realizable" coil has not been optimized, in the sense of 
obtaining minimum power consumption for the specified form factor. Thus, 
the design should be recycled until maximum attainable packing factor has 
been determined. (The solution can be made to converge rapidly. In fact, 
after some practice manipulating Appendices C and D, a better ''first guess" 
can generally be made than that suggested. ) 

In order to determine coolant volume fraction (and, for that matter, 
mean coil operating temperature) for a two-phase coolant, one must specify 


a particular combination of x and Ap. Presumably, the simplest 


ee 266 
procedure here would be to estimate the highest outlet quality which can 
safely be maintained (since coolant volume fraction decreases as quality 
increases), select a ''reasonable'’ A p (which can be related to coolant pump 
power through eqn. (34) ), and ''search" on Po» starting at one atmosphere 
and working upward. Coolant volume fraction generally (but not always) 
decreases with increased outlet pressure, due to increased vapor density 

in the coolant stream. The question here is whether or not the increased 
packing factor presumably attainable with higher outlet pressures decreases 
power consumption more than the higher coil operating temperature (due to 
higher AT.) increases power consumption. (Or, if not, and minimum size 
is being sought, if the higher loss may be tolerable, in exchange for smaller 


or 
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size.) Thus, the flow chart indicates that once the optimum" design (mini- 
mum power consumption, minimum size, or some compromise between 
these) has been determined for a particular Xo» Pop and Ap, one should 


specify a new combination of these parameters, and recycle the design. 


C. Example problem: minimum power consumption: 

As an example of the above design procedure, let us determine the 
total power consumption for an F=2, liquid-hydrogen-cooled, aluminum 
saddle coil for use with an MHD channel requiring B=90 kg, and Aun ook Il me 


1. From Fig. A-1, J, = 5.81x 107 B/VA_, = 5.81 x 10° (50/1) 
=9,2x 10° ne 


2, L=10VA_, = 3.16 m: therefore, from Fig. A-3, V = .176L° 


MSG mm, and A_ = 2.31 L  OSnen 


3. Assume 0. = 4500 p.s.i. and 04 = 60,000 p.s.i.: from Fig. C-2, 
o. V. FF. = . 033. 

4. Assume he = latm., Ap = 0.5 atm., and ee 90%. From Fig. 
D-8b (no-slip model), AT. = 0.9 °K, 

9. Kutateladze's correlation for peak heat flux for liquid hydrogen 
indicates that the critical pool boiling heat flux, for a pressure of ~ 1.29 
atm., is about 10 eRe)! cm at AT, = 2,8 °K. (See (112). ) Presumably, 
forced convection cooling would at least equal this. In keeping with the dis- 
cussion in Chapter V, we shall assume these figures apply here. 


6. Using the above, T = 20.3 + 0.9 + 2.8 = 24°K. 


3 8 


“7. From Fig. B-7, — 


= 4,76 x eo ohm-meters. 


2.69 x10 


Oe =a Oegoex 107 °)(1. 772 x) Oe) 


8. The total electrical power loss for r = 1 would be OVFI,” = 44.9 Kw. 


9. Assume A =.75-S.V.F. =.75-.03=.72. Corrected power loss 
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becomes 44, 9/.72 = 62.3 Kw. 


10. From Fig. D-7b (no-slip model), S = 1.16 x 1074. From eqn. 


(58), coolant volume is therefore: V_ = qS(FL) =(6. 23 x Monae, iG aka = 


) 
(632) oat = 4,56 x 107° ne and C. V.F. is of the order of .005/5, which 
is negligible. 

11. At this point, it would be well to re-examine the applicability of 
the "inviscid flow'’ model. The fact that the model calls for such a minute 
coolant volume fraction in this case is related to the fact that loss density 
scales with resistivity, and the resistivity of aluminum at this temperature 
is about three orders of magnitude less than that of copper at 300° K. as- 
suming q = 46. 3 Kw and q/A = 10, the required total heat transfer surface 
would be 4630 ene. Thus, the inviscid flow model indicates that a slot 
height of about (3390/4630) = 0. 732 cm is required, from energy balance 
considerations. However, this would result in an equivalent" L/D of 
(632/1. 464) = 432, which leads to the possibility that viscous effects may 
not be negligible, in this case. 

To a first approximation, pressure drop due to viscous effects may be 
considered directly additive to that attributed to the inviscid flow model. We 
can establish an upper and lower bound on the influence of viscosity by 
asking what additional A p would be required to move the specified flow 
if a) the coolant were entirely liquid, and b) the coolant were entirely vapor. 

The viscosity of saturated liquid hydrogen, for p~1 atm., is about 126 
micropoise, and the "equivalent diameter’ of the slot is 2(0. 732) = 1. 464 cm. 
From Fig. D-9b, the inlet velocity under the specified conditions is 5. 30 m/s. 


Thus, the Reynold's number for ''all liquid" flow would be: 


4 9) 


PVD/UL = (, 0707)(530)(1. 464)/(1.26 x 10 “) = 5.86 x 10 


"3 
and the corresponding friction factor (from (108) ) would be 3.20x 10 ©, 
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From eqn. (30), Ap = 4f(L/D)O v)" [(20 ). Here, Oy =e O70 T(S30)-= 34.57 
mad atm. - 1,013 x 10° dynes/cm’. Thus, 

APs = 4(. 00320)(432)(37. 5)” /(2)(. 0707)(1, 013 x 10°) = , 054 atm. 
Similarly, the viscosity of the saturated vapor would be about 11 micropoise; 
the Reynold's number for ''all vapor'' flow would be 6.71 x 10° (since pv 
= const for the specified flow), and the corresponding friction factor would 
be 2.25x 10°. Thus, 

AP a, = 4(. 00225)(432)(37. 5)” /(2)(1. 34 x 107°)(1, 013 x 10°) = 3,24 atm. 
Obviously, this flow is not "inviscid"! However, we can adjust our model, 
by demanding that the total pressure drop be of the order originally speci- 
fied. Very approximately, viscous pressure drop scales as (L/D)?: Let 
us try increasing D by a factor of 4. To maintain roughly the same outlet 
quality, we shall maintain constant mass flow rate by decreasing (Ov) by 
the same factor. The 'new'' bounds become: 


Ap ins 4(. 00320)(108)(9. 38)7/(2)(. 0707)(1. 013 x to) = , 000845 atm. 


ml 


Ap_ = 4(. 00225)(108)(9. 38)" /(2)(1. 34 x 107°)(1. 013 x 10°) = . 0506 atm. 


ma 
This more acceptable! (Ironically, coolant volume fraction now becomes 

. 0135/5. 56 = . 002, which is still negligible!) With the corrected figures, 
maximum attainable packing factor is of the order of 0. 96, and corresponding 
electric power consumption would be 46. 8 Kw. 

12. Assuming surface heat leak is of the order of .0o watts/m“, total 
heat leak from this source is .05(23. 1) 1 watt, which is negligible. 

13. Assume that the thermal conductivity of the support structure is 
about 0.1 watts/cm-°K, and that the length of support elements is of the order 
of the coil radius (which, from Fig. A-2, is .186 L = 0.588 m). Let the 
specified stress limit in the supports be 30,000 p.s.i. The net cross sec- 
tion required to support the coil mass would be of the order of: 
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2 2 


Mg _ 2700(5. 56)(2.205) 1 1 in? = 71 om? 


7 O 30, 000 
Since the coil would presumably be supported in all three directions, let 
us increase this figure to ~ 29 sian. Then, the net structural heat leak 
would be of the order of . 1(25/58. 8)(276) = 11.68 watts, which is negligible, 

14. The flow energy supplied by the pump represents hydraulic dis- 
sipation, and presumably is removed as heat by the refrigerator. Assuming 
a pump efficiency of 50%, this is of the order of 2 (F.¥) - 2(V A p/L)v, .. 
Using the adjusted flow figures, dissipation from this source is of the order 
of 2(. 0135). 5 x 1.013 x 10°)(1, 325)/6, 32 = 287 watts. 

15. Although the coil windings are at about 24° K, the refrigerator 
must remove heat from the coolant at 20.3° K. From Fig. 18, refrigerator 
power is about 70 watts/watt. Total operating power required is thus 
71(4, 71 x 10%) = 3.34 Mw. 

16, If the hypothesized MHD channel were utilized for a seeded-gas 
MHD generator, one might assume a density of power generation in the 
channel of ~ 150 Mw/m°. Since channel volume is 0. 316 m°, SuUchearcene 
erator would produce about 47.4 Mw of electrical power. Thus, the total 
operating power consumption of this coil would amount to about 7% of the 
generator output. Similar calculations indicate that at this temperature, 
all 50 kg, F=2 saddle coils share the basic characteristics discussed above: 
i,e., coolant volume fraction, heat leak, and hydraulic dissipation do not 
significantly influence the design, within the range of sizes covered by this 
investigation. Further, since structural volume fraction has been found 
to be a function of B, F, 0. and Og only, packing factor is virtually inde- 
pendent of coil size. Thus, since coil volume and channel volume scale 
together for constant F, and Je scales as Ua, for constant B, the per- 


centage of generator output consumed by such coils would scale as 1/A p> 
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for B=50 kg, F=2 saddle coils, in the range . 01 = 10 square meters. 

Several factors bear discussion at this point. The first is that one might, 
or might not, ever realize the theoretically attainable packing factor of 0. 96. 
Probably, insulation requirements, and simple engineering convenience, 
would reduce this to about 0.7. However, it was previously indicated that 
for coils of the same packing factor, the power consumption of a cryogenic 
aluminum coil would be about an order of magnitude less than that of a com- 
parable "conventional" coil. With loss densities running three orders 
higher, it does not appear likely that ‘conventional coils could approach 
the packing factors attainable with near-optimum-temperature cryogenic 
coils. Thus, whether or not the theoretical limit of 0.96 were attained, one 
might expect that a relative assessment of "'attainable'' cryogenic and con- 
ventional coils would favor the former, even more than might be assumed 
from Fig. 23. 

Next, even if the packing factor realized in an ''as built" coil were appre- 
ciably lower than that predicted by the theoretical model, the underlying 
physical characteristics of the design should remain roughly as described 
here. Most specifically, one would expect the percentage power consumption 
of ‘minimum loss'' MHD generator coils to scale inversely with channel 
cross section, under the geometric assumptions made here. Under these 
circumstances, it appears that such coils would be suitable for MHD power 
generation, with generators of 50 Mw output, or larger. Since MHD gener- 
ators must probably produce well over 50 Mw to be commercially compet- 
itive, due to other considerations, such coils would appear highly attractive 
for power generation applications. 

Finally, since coolant volume proved to be such an insignificant item in 
this example problem, it is apparent that one need not strive for maximum 


attainable q/A in cooling system design. Our ''conservative’ estimate was 
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that we could obtain a heat flux of 10 ete ena ata AT. Ol Zac he 
Suppose this AT) could be reduced to, say, l a a either through achieve- 
ment of an improved heat transfer coefficient (through forced convection), 
or through the luxury of providing additional heat transfer surface, to re- 
duce q/A. Presumably, we would have two choices: either to operate the 


refrigerator 1.8° "hotter", or to operate the coil 1. 8° " 


colder". Assuming 
that attainable refrigerator efficiencies do, in fact, conform to some smooth 
curve, and are not a ''step function", Fig. 18 indicates that the former 
option would reduce refrigeration power to about 61.6 watts/watt, and there- 
fore reduce total power consumption by about 12%. The latter option would 
reduce mean coil resistivity to about 3.6 x oes ohm-centimeters, and 
therefore reduce total power consumption by about 20%. This again reflects 
the marked temperature sensitivity of the optimum operating point, and 
indicates that one would do well to experiment with different heat transfer 


solutions. It also indicates the importance of obtaining realiable heat trans- 


fer data for the particular surface and geometry contemplated. 


D. Example problem: 'lightweight"' coil: 
For comparison with the above, it might be interesting to determine the 
characteristics of an F=1. 3, liquid-hydrogen-cooled, aluminum saddle coil 


for the same application: 


r/ 2 


a/m“, and V becomes 


veo 16 m°, From Fig. B-5, mean coil resistivity is 5.11 x Noe ohm-meters, 


1. From Appendix A, oe now becomes 7. 3 x 10 


and from Fig. C-2, the structural volume fraction becomes . 145. Other 
significant 'input'' parameters are unchanged. 


2. Electrical power consumption for Xr =1 becomes 111.7 Kw. 
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“3. Slot factor is unchanged, but assumed slot length now becomes 
FL=1, 3(3.16) = 4.11 m. Coolant volume becomes 6230 om, and Ger. 
becomes . 02, for a net packing factor of - 835, and resultant electrical 
loss of 134 Kw. 

4. Checking slot height, L/D=445. We again must enlarge the slots 
by a factor of 4, to reduce L/D to about 110. Adjusted coolant volume 
fraction becomes . 08, for a net packing factor of - 829, and an electrical 
loss of 135 Kw. 

5. Assuming the same refrigerator efficiency, total power consumption 
would be 9.58 Mw. 

To compare this coil with that preceding: the volume of the ''minimum 
loss'' (F=2) coil was 17.6 times as great as that of the MHD channel. By 
contrast, the volume of the "lightweight'' (F=1. 3) coil is equal to that of 
the MHD channel, The former coil would weigh approximately 16.5 tons, 
the latter, about 0.93 tons. This reduction in size and mass (and presumably, 
initial capital cost) has been purchased at the expense of an increase in 


power consumption (and therefore, operating cost) by a factor of 2. 9. 
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IX. CONCLUSIONS 


1. ''Minimum loss"' cryogenic saddle coils appear economically fea- 
sible for use with MHD generating channels with an output of 50 Mw or 
greater. The total power requirement of such coils comprises less than 
10% of generator output at 50 Mw, and the percentage of output power con- 


sumed by the coil scales approximately as Lj Ape for the specified geometric 


model, 


2. The optimum conductor for such saddle coils is aluminum. The 
optimum operating temperature (to minimize total power consumption) is 
about 16 °K. At this temperature, total power consumption of cryogenic 
aluminum saddle coils would be about one order of magnitude less than that 
of comparable ''conventional'' copper saddle coils, with the same packing 
factor. However, reduced loss density (the resistivity of aluminum at 
this temperature is about three orders less than that of ''standard''’ copper 
at 300 °K), and the increased strength of structural materials at cryogenic 
temperatures, should permit realization of much higher packing factors 
in the cryogenic coils. Since total loss varies inversely with packing fac- 
tor, the comparison between attainable cryogenic and ''conventional" coils 
favors the former, even more strongly. 

3. Temperature optimization is fairly critical: a shift of 10° K in 
either direction will cause, roughly, a 50% increase in total operating 
power requirements. Therefore, supercritical helium or liquid hydrogen 
would be the preferred coolants. (A recent NBS study (105) indicates that 
optimum operating pressures for the former coolant fall in the range of 15 
to 20 atm.) At constant packing factor, liquid-helium and liquid-neon-cooled 
coils would be about equally remote from the optimum, in terms of total 


-_ 
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power consumption. However, the low heat of vaporization and low heat 
flux attainable with the former might lead to lower packing factors. What- 
ever the coolant used, geometric considerations and coil size indicate that 
a forced convection cooling system would be required. The advantage most 
frequently cited for the use of liquid neon for immersion-cooled coils--- i.e., 
better packing factors, because of higher latent heat of vaporization, ona 
Seiametric basis--- is inapplicable to forced-convection cooling systems. 
The low sonic velocity reported for neon sharply restricts attainable flow 
rates, due to choking considerations. Thus, for forced-convection cooling 
systems, given the same heat flux, outlet quality, etc. , liquid hydrogen 
permits higher packing factors than liquid neon. 

4. The "shape" of the saddle coil model has been quantized in terms of 
the form factor, F = 1+ L/L, which compares the mean length of the end 
turns (assumed to contribute nothing to the "useful" field) to the 'working" 
length of the coil. Forms of interest fall in the range 1, 3 = F220 
"minimum loss’ coil is that corresponding to F=2, at which point coil vol- 
ume is nearly 18 times the volume of the associatedMHD channel. If F is 
increased above 2, both coil size and total loss increase. Obviously, such 
shapes are not of interest. If F is decreased below 2, total loss increases, 
but coil size decreases: at F=1. 3, coil volume equals the volume of the 
associated MHD channel, and losses are higher by about a factor of 3 (for 


90 kg, liquid-hydrogen-cooled, aluminum saddle coils). Minimum attainable 
F is largely determined by structural considerations and falls roughly in the 


range 1.3 =F , = 1.6, depending upon the application. Mean coil resistivity 


gai. 
is a very weak function of F, except for copper coils below 60 7K, 


o. otructural Volume fraction has been found to be a function of the 
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desired channel flux density, the coil form factor, the specified structural 
stress limit, and the specified conductor stress limit, only. It is a fairly 
strong function of the first three, and a weak function of the fourth. Basi- 
cally, the conductor stress limit determines the number of structural ele- 
ments required, and the structural stress limit, their thickness. Since the 


‘cold work", it is essential that 


resistivity of aluminum is sensitive to 
stress levels in aluminum cryogenic conductors be kept below the propor- 
tional limit for that material. Since yield strength is generally decreased 
by the same metallurgical processing used to minimize the resistivity of 
cryogenic conductors, considerable thought must be given to the design of 
structural supports for aluminum coils. If possible, aluminum coils should 
be designed to permit final anneal after assembly, or at least, anneal of 
individual conductor "bundles" after final shaping. 

6. At or near optimum temperature, heat removal from cryogenic 
aluminum coils is not a significant problem. For example, for a liquid- 
hydrogen-cooled aluminum coil, the heat transfer surface and coolant 
volume required from an energy-balance viewpoint are insignificant. How- 
ever, coolant ducts must be distributed through the windings, in sufficient 
number of ensure reasonable temperature homogeneity. Analysis of coil 
geometry and stress distribution indicates that the most likely coolant duct 
geometry would involve the use of slots (vice tubes, pipes, etc. ), located 
adjacent to structural support members. Assuming the support structure 
would be stainless steel, the thermal conductivity of that structure would 
be two to four orders of magnitude less than that of the coil conductors, at 
the temperatures investigated here. Therefore, it seems reasonable to 
assume that heat flux into the coolant would occur only across the slot 


surface ‘away from'' the support structure. 
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7. At or near optimum temperature, heat generation may be so low 
that the coolant slot height must be ‘sized’ by consideration of viscous 
effects, rather than by use of the inviscid-flow, energy-balance models pre- 
sented herein. However, when such is the case, coolant volume fraction 
tends to be so small that the choice between these alternatives does not 
strongly influence the ultimate outcome of the design analysis. 

8. At or near optimum temperature, the performance of aluminum 
cryogenic coils is dominated by refrigeration costs. In this study, it was 
estimated that currently, 70 watts of refrigerator power would be required 
to remove each watt of heat generated at 20. 3 °K. Of the various heat 
sources, electrical losses dominate. (Fora 50 kg, F=2, liquid-hydrogen- 
cooled aluminum coil, suitable for use with a 47 Mw MHD generator, elec- 
trical losses were computed to be 46.8 Kw, surface heat leak, 1 watt, 
structural heat leak, ~6 watts, and hydraulic dissipation, ~ 300 watts. 
Heat leak through the electrical input connections was not estimated, but 
presumably, proper design of these leads can ensure that this leak is neg- 
ligible, compared to lag heating.) An ideal, Carnot refrigerator would 
require only 13.8 watts of power for every watt of heat 'pumped"’ from 
20. 3° to 300° K. There is obviously ''room for improvement", and past 
trends indicate that further improvement in refrigerator efficiency may be 
expected. Any improvement in refrigerator efficiency would obviously be 
reflected in reduced operating costs. 

9. Due to the strong coupling between the residual resistivity of 
aluminum, and the other components of ''total'’ resistivity, improved met- 
allurgical techniques might yield a significant reduction in the resistivity 
of that metal, at 16° K. (In particular, the process which Post and Taylor 
have utilized to reduce the residual resistivity of sodium (63) may hold 


some promise.) If the Bloch-Gruneisen curve does, in fact, serve asa ’ 
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lower-bound for the resistivity of aluminum, a net reduction of about a 
factor of two may be possible. 

10. Considering 8. and 9. above, one might reasonably hope for an 
eventual reduction in the minimum power requirements of cryogenic saddle 
coils, by as much as a factor of 4 below that currently attainable. 

11. This study has been limited to 'air-core' coils, in keeping with 
the model defined in Chapter I. However, it is possible that the basic 
saddle coil configuration, for a smaller F=2 coils investigated here, could 
feasibly be modified for use with an iron core. Ina sense, this might buy 
the first 20 kg ''for free''’. By such means, 50 kg saddle coils might be 
made attractive for use with smaller MHD generators than the 50 Mw mini- 


mum established here. 
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X. RECOMMENDATIONS 


Every effort has been made to ensure mutual independence between the 
various sections of this work, in order that the value of the whole need not ° 
depend upon complete acceptance of each component part. Specifically, 
there are several areas where existing data are scanty, or controversial, 
or where the ''state-of-the-art'' is in a process of apparent change: 

1, Magnetoresistance is known to be sensitive to size effects, and one 
might be advised to obtain experimental data for the conductor geometry 
contemplated. If this is done, the possible inaccuracies inherent in ‘pulsed’ 
measurements should be weighed carefully, before proceeding. If possible, 


such data should be obtained before entering the detailed design phase. 


2. Cryogenic refrigerators have been improved significantly within 


the last several years. Since the early 1960's, the projected power con- 
sumption of large, liquid-helium refrigerators has dropped from 1000 watts 
per watt to 500 watts per watt, and further improvement seems likely. Here, 
the best source of information would be the various manufacturers' repre- 
sentatives. 

3. Heat transfer is probably the greatest point of uncertainty in this 
work, While the heat transfer problem does not dominate optimum” saddle 
coil design, both resistivity and refrigerator efficiency are sufficiently 


strong functions of temperature that AT,, and AT. should be established 


bl 
with some accuracy, if possible. Presumably, one could afford the luxury 
of designing for low heat flux and low pressure drop, to minimize these tem- 


perature increments. With two-phase coolants, it would be valuable to know 


the maximum heat flux which could be tolerated, without risking burnout. 
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Again, experiments to obtain these data, for the coolant and duct geometry 
contemplated, should be conducted as early in the design process as possikle. 
In particular, experimental data for forced-convection, two-phase cryogenic 
boiling, with the type of slot goemetry discussed here, are nonexistant at 

the present time. 

In short, if better information can be obtained, in any of these areas, it 
can readily be substituted for the pertinent section of this work, and the 
analysis modified accordingly. 

Finally, one point might bear repeating, in closing. The only justi- 
fication for the added cost and complexity of cryogenic coils, as compared to 
"conventional" coils, is that of reduced operating power requirements. Even 
if an optimum-temperature, F=2 coil is considered undesirable for a par- 
ticular application, the basic motivation for using some other cryogenic 
coil, rather than a ''conventional” coil, is still the same. Logically, then, 
the design analysis should sieze every opportunity to ensure that the cryogenic 
coil specified delivers its full potential. For example, the use of metal 
foils is convenient in the construction of ''tape-wound' solenoids. However, 
if preliminary experiments indicate that the magnetoresistance of such foils 
lies appreciably above the bulk of amassed experimental data, then that con- 
cept should be discarded, and some other conductor geometry investigated, 
before proceeding to the construction stage. Considering the extreme pen- 
alty paid for every watt of heat dissipated within the coil, the expenditure 
of additional time, effort, and capital before building the coil should be 


easily justifiable on the basis of reduced operating cost! 
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APPENDIX A 


Coil Geometry Curves 


List of Symbols 


a Offset between center of channel and gs of either coil cylinder (m) 
A ch MHD channel cross-sectional area (m? ) 


A el Total coil cross-sectional ee (both lobes) faa ) 


as 


Total coil surface area (m* ) 
Channel flux density (KG) 


S 


B 

J Effective current density (amps/m°) 
IL, Channel (or coil) length (m) 

R Coil cylinder radius (m) 

V Total coil volume (m°) 

W 


a MHD channel width (m) 
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FIGURE A-1 


Effective Current Density for given Form Factor 
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ho ayn ang es for given Form Factor 
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FIGURE A -3 


As Aap and V for given Form Factor 
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Resistivity Curves 


List of Symbols 


Channel flux density (KG) 

Coil form factor (dimensionless) 

Mean coil resistivity (ohm-meters) 

Standard resistance: 1.7715 x fone ohm-meters 


Mean coil temperature (°K) 
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FIGURE B-1 
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for B= 29 kg 
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FIGURE B-2 


R/R, for B = 25 kg, F = 1.6 
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FIGURE B-3 


R/R, for B = 25 kg, F = 2,0 
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FIGURE B-4 


R/R, for B= 25 kg, F = 2.4 
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FIGURE B - 6 
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FIGURE B - 7 


for B = 50 kg, F = 2.0 
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FIGURE K-38 


R/R, for B = 50 kg, F = 2.7 
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FIGURE B-9 
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for B= 75 kg, F = 1.3 
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FIGURE B-10 





for B= 75 kg, F = 1.6 
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FIGURE p-1i1 


R/R, for B= 75 kg, F = 2.0 
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FIGURE B-12 


for B= 75 kg, F = 2,3 
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FIGURE B-13 


R/R, for B= 100 kg, F = 1.4 
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FIGURE B-14 


100 kg, F = 1.7 


R/R_ for B 
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FIGURE B-15 


R/R, for B = 100 kg, F = 2.1 
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FIGURE 4B-16 


fons = 100%ke, F — 2.7 
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APPENDIX C 


Structural Curves 


List of Symbols 


B Channel flux density (KG) 
F Coil form factor (dimensionless) 

S.V.F. Structural volume fraction (dimensionless) 
C. Specified stress limit for coil conductors (kpsi) 


Oq Design stress limit for coil structure (kpsi) 
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APPENDIX D 


Coolant Curves 


List of Symbols 


Py Coolant outlet pressure (atm. ) 
Ap Coolant pressure drop (inlet to outlet) (atm. ) 
S Slot factor (cm /watt) 


AT_ Mean temperature increment to correct for deviation of temperature of 
bulk liquid from l-atm. saturation temperature (°K) 


in Coolant inlet velocity (m/s) 
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MATN PROGRAM (COIL GEOMETRY/RESISTIVITY ANALYSIS) 


DIMENSTON TEMP(10) sRTEMP(10) »RMAG(995) »BTAB( 100950) sAVRHO(10) 
COMMON ACHANs WCHANs HCHANs ALTHs ACOILs RADIUS» OFFSETs TEMP s 
IRTEMPs RMAGs BIABs AVRHOs RSQs LBITABs VCOILs SCOIL»s ENDFAC 
motel READ 

READ 2s Als Bs CUR 

FORMAT(2F1365s E1344) 

meet SIZE (Als Bs CUR) 

Beet FIFLD (CUP) 

meet RHO 

CALL PPINT (AlsBsCUR) 

eo 1 1 

FAI 
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SURROUTINE RFAD 


DIMENSTON TEMP(10) »RTEMP( 10) sRMAG(995)2BTAB(100950) sAVRHO(10) 
COMMON ACHAN»s WCHANs HCHANs ALTHs ACOILs RADIUS» OFFSET» TEmp, 
IRTEMPs RMAGs BIABs AVRHOs RSQs LBTABs VCOIL» SCOIL»s ENDFAC 
RFAD ls (TEMP (1)s7T=1910) 

FORMAT(19174) 

READ 2s (RTEMP(1)sT=1910) 

FORMAT(5E11¢4) 

RFAD 3s ((RMAG(Is9J)sT=199)sJ=195) 

FORMAT (5E6114¢4/4E11¢4) 

PFETURN 

FNN 





BUBROUTINE SIZE (Als Bs CUR) 


DIMENSTON TE’P (10) x9 TEMP( 10) sPMAG(995) sBTAB( 100950) sAVRHO(10) 
COMMON ACHANs WCHANs HCHANs ALTHs ACOIL»s RADIUSs OFFSET» TEMP, 
IRTFMP >» RMAGs BIABs AVRHOs RSQ» LBTABs VCOITL» SCOIL» ENDFAC 
=A ) 

OFFSET=8/(1.25664E-5*CUR) 

Mmee4 [=1550 

RANTUS=OFFSFT+SQPTF(A)/2¢ 

MSQ=RANTUS#RPANIUS 

S=SORTFI(OFFSET#OFFSET+8.#RSQ)—-32*#OFF SET 
ACHAN=S*SOQRTF(RSQ-(OFFSET+H+S/4e)* (OFFSET+S/4.) ) 
DELTA=A1—ACHAN 

IF (A1/190000,-ABSF(DELTA)) 29595 

IF (NFLTA) 329594 

A=zA-O..9*DELTA 

A=A+DELTA 

WCHAN=S/2.6 

HCHAN=ACHAN/WCHAN 

ALTH=10 e*SQRTF (ACHAN) 

GAMMA=ACOSF(OFFSET/RADIUS ) 

ACOIL=2 e*RSQ* (3e14159—-( 2e X¥GAMMA—SINF(2e*GAMMA ) ) ) 
MeQtTL=ACOIL#ALTH 

weep =12.56637*RADIUS*#ALTH 
ENDFAC=3~-14159*(RADIUS+OFFSET)/ALTH + le 

RETURN 

FAD 
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meeRQUTINF FIFLOD (CUR) 


DIMENSTON TEMP(10) »RTEMP(10) »RMAG(995) sBTAB( 100950) sAVRHO( 10) 
COMMON ACHANs WCHANs HCHANs ALTHs ACOITL» RADIUS» OFFSET» TEMP, 
IRTEMPs RMAGs BIABs AVRHO»s RSQs LBTABs VCOIL» SCOIL» ENDFAC 
De} 1=15100 

Be 1 J=1550 

mapBi(lsJ)=0, 
X1=PANTUS-OFFSFT 
X2=RANTUS+0F FSET 
DFLTA=RADIUS/50. 

mer =PFLIA/s2-. 
BMAG=60e28318E=-/*CUR 

meee.) f=159100 

AT=] 

men 1 *DFLTA-SET 

IF(X2-X) 119292 

XMO=X-OFFSET 

XMONSQEXMN#XMNO 

XPO=X+0OFFSFT 

XPOSQ=XPOXXPO 
W2Z=SORTF(RSQ=XMOSQ) 

TF(X—=X1) 3797 
Y1=SQRTF(PSO=]XPOSQ) 

NN 6 JH1559 

AJ=J 

men JHeNFLTA-SFT 

IF(Y=Y1) 69494 

MPC Y2-Y) 1099535 

RATTO#R2 SN/(XPOSOQ+Y*Y) 
BX=BMAG¥(1-e-RATIO)*Y 
BY=BMAG# (XPOXPATTO=XMO) 
RTAR( TsJ)=SQRTF(RX*#¥RX4+BYXBY ) 
PONTINUF 

eo) TO 10 

moe 9 J=1559 

AJ=J 

Sean IJeDELTA=-SFT 

IF(Y2-Y) 109898 

ef O=RSQ/(XPOSQ+Y*Y) 
BX=BMAG*(1e¢-RATIO) *Y 
RBY=RMAG*(XPO*#RAT TO—XMO) 
BTAB( IsJ)=SQRTF( BX*¥BX+BY*#BY ) 
m7) AB=] 

RPETURN 

FAD 
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SUBROUTINE RHO 


DIMENSION TEMP(10) »RTEMP( 10) »RMAG(995) 9BTAB( 100950) sAVRHO(10) 
COMMON ACHANs WCHANs HCHANs ALTHs ACOILs RADIUS» OFFSET» TEMPs 
1RTFMPs RMAGs BTIABs AVRHOs RSQ» LBTABs VCOILs SCOIL»s ENDFAC 
BaA=2S59/7500. 

NN 8&8 T=1910 

NASUM=0, 

PDA moM=0) ¢ 

Bo 7 J=lsLBTAR 

NO 7 K=1959 

TF(RTAR(U9K)) 79791 

ial 

M=1 

IN=19.*BIAR(JU9K)/RTEMP(T) 

TFC TN=$19) 29494 

RNASUM=RDASUM+RTIFMP(1)*DA 

mo 1!) 6 

IN=IN/10 

M=M+ ] 

IF(1O00=IN) 39395 

b=(1TN+5)/10 

RDASUM=RDASUM+RTEMP(1T)*#(RMAG(L9M)+1¢)#*DA 

DASUM=DASUM+DA 

CONT INUF 

AVRHO(1T)=RDASUM/DASUM 

RETURN 

FAD 
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SURROUTINE PRINT (AlsBsCUR) 


NDIMENSTON TE&MP(10) »9RTEMP(10) »9RMAG(995)9BTAB(100950)sAVRHO(10) 
COMMON ACHANs WCHANs HCHANs ALTHs ACOILs RADIUSs OFFSETs TEMP, 
YRTEMP»s RMAGs BIABs AVRHOs RSQs LBTABs VCOIL» SCOIL»s ENDFAC 

Mme (LIST) 1ls193 

LIST=4 

PRINT 2 

FOPMAT(1H1) 

PRINT 4s AlsReCUR 

FORMAT (1HC/34HO DESIRED CHANNEL AREA (M2) =9F9e4s6Xs26HDESIRE 
Meme UX DENSTTY (KG)=9F5e1296Xs32HDESIRED CURRENT DENSITY (A/M2) =s 
o-9.3) 

PRINT 5S» ACHAN »WCHAN sHCHAN 

FOPMAT( 35H COMPUTED CHANNEL AREA (M2) =9F8e496X923HOPT CHANN 
BEL WIDTH (M) =sF8e496Xs28HOPTIMUM CHANNEL HEIGHT (M) =9F13e4) 
PRINT 69 ALTHsRADIUS sOFFSET 

FOPMAT (33H COMPUTED COIL LENGTH (M) =9F100¢496Xs21HCYLINDER R 
MmeprUS (M) =sFiO0e4s6Xs21HCYLINDER OFFSET (M) =9F20e4) 

meen? 7s ACOIL»sVCOIL»sSCOIL 

FORMAT (33H TOTAL CROSS SECTION (M2) =sFl0.4s6X%s SRO oie 
IME (M3) =9F126¢4s6X925HIONTAL SURFACE AREA (M2) =9F160e4) 

PRINT 8» FNNDFAC 

FOPMAT(1IHOs48Xs19HCORRECTION FACTOR =93F 124) 

PRINT 9 

FORMAT( HO s30X972HMEAN COIL RESISTIVITY AS A FUNCTION OF TEMPERATU 
1RE FOR ALUMINUM WINDINGS) 

PRINT 1Os(TEMP(1T)sT=1910) 

FORMAT(12HO T (DEG K)9101I12) 

PRINT lls (AVRHO(7)sT=1910) 

FORMAT(12H RHO RATIO910E12¢4) 

mo l=LIST=1 

RFTURN 

ENN 
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MATN PROGRAM (STRUCTURAL ANALYSIS) 


COMMON AlsBsCURsSIGMAC »sSIGMAD sALTHsACOIL sRADIUS sRSQsOFFSET s9ENDFAC 9 
INCOMP sVCOMP sCFRACT sNTENSsVTENS os TFRACT sVFRACT 

PRINT 1 

FORMAT(1H1) 

READ 3s AlsBsCUR 

FORMAT (2F134e59E1304) 

met SIZE 

PRINT 49 AlsBseCUR sRADIUS sENDFAC 

FORMAT(1HOU/7H ACHAN=9F60293H M2910X s92HB=9F60e1l93H KGslO0Xs2HJ=9E9e39 
15H A/M2910Xs2HR=9F8e492H Ms1lOXs7THENDFAC=9F7e4/1H ) 
DO 7 T=194 

Al=]I 

SIGMAC=AI*¥3000¢ 

‘DO 6 J=Hls4 

AJ=5—-J 

SIGMAD=AJ*300006 

ment CSHELL 

meet TSHELL 

IF (VFRACT=—e8) 598797 

IF (VFRACT) 89896 

BON TINUE 

CONTINUE 

mo TO 2 

FND 
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meeeQOUTINE SIZE 


COMMON AlsBsCUR»SITGMAC »sSIGMAD sALTHsACOIL »sRADIUS »RSQsOFFSET sENDFAC 9 
INCOMP »VCOMP sCFRACT »NTENSsVTENSosTFRACT »VFRACT 
a= A | 

MPFSET=B/(1.25664E-5*CUR) 

DO 4 1=1%950 

PANIUS=OFFSFT+SQRTF(A)/2.6 

RSN=RADTUS#RADIUS 
S=SQRTF(OFFSFT*#OFFSET+8.2*RSQ)-3-*0F FSET 
ACHAN=S*SQRTF (RSO-(OFFSET4+S/42.)*(OFFSET4+S/4,)) 
DELTA=A1=-ACHAN 

MomA1/10000,-ABSFI(DELTA)) 23595 

IF(DFELTA) 39594 

A=A--e9*PFLTA 

A=A+DELTA 

me e=10,.*SQRTF CACHAN) 

GAMMA=ACOSF(OFFSFT/RADIUS) 

ACOITL=2 e*RSQ*(3214159-( 26 *GAMMA-SINF(2¢*GAMMA ) )) 
FNDFAC=3e.14159¥*¥ (RADIUS+OFFSET)/ALTH + le 

RETURN 

FANT 
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BURPROUTINE CSHELL 


COMMON AlsBeCURsSTGMACsSIGMADsALTHsACOIL sRADIUS »sRSQsOFFSET sENDFAC + 
INCOMP »VCOMP sCFRACT sNTENSsVTENSoTFRACT sVFRACT 

VCOMP=0, 

NCAMD =f) 

RL IM=SORTF(2e*SIGMAC#1201325E+05/(14.696%1425664E-06))/CUR 
IF (RLIM-RADIUS) 1916916 

PIN=RLIM 

ROUTHERINKX(C (CL etG0#SITGMAC#1 201325E405/ (14 6696%*1 e25664F-O6*CURXCURER I 
mNeRIN)) **(12/306)) 

IF (ROUT-RADTUS) 339394 

PRESS=SIGMAC 

eo 19 & 

PREFSS=14,696%1 e25664E-O6*#CUR*FCUR#(RSQO#RADIUS“-RIN#®RIN*RIN) 7(1201325 
1—£+05*6e*RIN) 

Met SITGMAN/37500.)*(ALTH/(2.#RIN) ) 

TF (X=296) 6910910 

m= 25 25*PREISS/STGMAD 

Beret =7,532F+O07*(TNIF#2,5)/ CALTH/(2e*RINJ-—-O045*SQRTF(TN1)) 
IF (PCRIT—PRFSS) 899939 

TNL=TN1I¥((PRESS/PCRIT)**04) 

me) [0 7 

Meee -X) 19513913 

miP= (PRESS/50,186F+06) **(1./3-6) 

DLIM=.96386*SIGMAD#TN2 

IF (PLIM-PRFSS) 11912912 

TNL IM=(PLIM/50e186EF+06)## (10/36) 
BmoPF=150.,558F+06*TNLIMXTNLIM 

TINTH=TNLIM-PLIM/SLOPE 

PINT==$StOPEX TINT 

MNA=(PRESS=—PINI)/SLOPE 

TF (X-20,.) 12914514 

IF (TN1=-TN2) 13914914 

MEOMP=(2.*RIN)*#INI 

mo 17) 15 

TCOMP=(2,*RIN)*IN2 

VCOMP=VCOMP+4+3 .14159%(2e*RIN)*IT COMP 

NCOMP=NCOMP+4+] 

PIN=ROUT 

IF (RIN-RANIUS) 2916916 

CFRACT=VCOMP/(3414159*RSQ) 

RETURN 

FN 
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BUBROUTTNE TSHELL 


COMMON AlsBsCURsSIGMAC sSIGMADsALTHsACOIL s RADIUS sRSQsOFFSET sENDFAC 5 
INCOMP sVCOMP sCFRACT »sNTENSos VTENS 9 TFRACT s»VFRACT 

VTFNS=0. 

NTENS=9 

m=aRANTUS 

me=., LXOFFSF IT 

miei M=RADIUS+2.*OFFSET 

Mmmonstl=1.25664E-06*RSQ*CUR*CUR 

FSUM=0. 

K=(R#R-RSO+42*OFFSET*OFFSET)/(4.*OF FSET ) 

mG LE=ACOSF(X/R) 

NE=FCONST*#ANGLE#DR 

FSUM=FSUM+DF 

P=R+NR 

mee CR=RLIM) 339696 

PRESS=14.696*FSUM/(1.01325E+05#2.*RX ANGLE ) 

IF (PRFSS=SIGMAC) 253594 

R=R-(FSUM=24.*R*ANGLE*SIGMAC#] -01325E+05/14.696)/(FCONST*ANGLE ) 
TTFNS=R*¥SITGMAC/SIGMAD 

VTENS=VTEFNS+2.*R¥ANGLF*¥TTFNS 

NTFENS=NTFNS+1 

So 19 1 
mee ACTH2Z.*VIFNS/ACOIL 

WEPACT=CFRACT+IFRACT 

PRPINT 7s SITGMACsSIGMAD»VFRACT »NCOMPsCFRACT sNTENSs TFRACT 

FORMAT( 8H SIGMAC=9F8els4H PSI98Xs7HSIGMAD=9F9e194H PSI 28Xs7HVFRACT 
L=9F60398Xs 3HNC=91398Xs BHVC= 9F 60 398X989 3ZHNI =91398X9 BHVI = oF Ge 3) 

PFTURN 

FAD 
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SAMPLE MAIN PROGRAM (DATA FITTING) 


DIMENSION X(200)sY(200) sYGEN( 200) sDELTA( 200) sPRCENT( 200) sA(20097)5 
me (7)sIL(7) 

RFAN ls L 

FORMAT(T3) 

mea 2s (X(JT)s YC(I)s TtlesL) 

FORMAT(2F10.44) 

NO 3 JTeHlolt 

memes 1 )=1.0 

A(Ts2)=X(T) 

A(Ts3)=ACT92)4X( J) 

Mrs G4)=A(Ts2)*X(T) 

Mets 5 )=AlCT94)#X(7) 

Mr 96) =A T95)¥*X(T) 

pays / )=Y(T) 

M=a& 

Bs [L 

m=. 1]F-5 

CALL GLSQ(AsZsTL»NoaMsALPHAsE19E1 ) 

DO 4 J=HloeL 

mee—=ACT)*X()) 

YGFN(1)2=Z2(01)4202)*X07T)42(3)#X242(4)*X2#X(7T)42(5) ¥X2#X24Z2(6)4X2*#X2% 
mX(T) 

DELTA(TI)=Y(I)-YGEN(T) 

MeGeNT(7)=100,.*DFLTACI)/YGEN(T) 

SUM=L 

RMS=SORTF(ALPHA*ALPHA/SUM) 

PRE INT 55 PMS 

FAPMAT(1H196Xs120HASSUMING Y = ACO) + Al1)X + Al2)X2 + AL3)X3 + Al 
mens + A(5)X5s THE BELOW COEFFICIENTS MINIMIZE RMS DEVIATION 10 A 
BYVALUE /1H s6Xs3HOF 9E14e681lHe) 

PRINT 65 (Z(7T)sT=H194) 

FORMAT(1HOs6Xs7HA(O) = 9E14e6910Xs7HA(1) = 29E14e6910X97HA(2) = 9El 
14.6310Xs7HA(3) = 961406) 

PRPINT 7s 2(5)s 2(6) 

FORMAT( 1H 937X97HA(G) = 9E14e6910X97HA(5) = 9£E1406) 

PRINT R 

FORPMAT(1HO/1HO942X »47HGENERATED AND OBSERVED DATA ARE COMPARED BEL 
1OW, ) 

PPINT 9 

FAOMAT( JHO922KXs1THX918Xs]TOHY OBSEPVEDs14X911HY GENER ATEDs14X s9HDEVI 
TATTONS10Xs 7HPFRCFNT /1HO) 

PPINT 105 (X(T)s YlT)9 YGEN(CT)s DELTA(C TI)» PRCENT(I)s I=l19L) 
FOPMAT(1H 914X9E14e¢63910XsF14e6910X SELIG e6s1lOX9E14e6910X 9F 702) 

GO TO 109 

END 
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4) 


35 


3] 


42 
39 
pI 
59 


52 


SURROUTINE GLSQ(AsXsTLeNeMsALPHA9E1 9E2) 


DIMENSTON A(20087) 9X(7)sTL(7) 
MM=M+1] 

a = |! 

AN 60 J=],MM 

me J)=0 

T=1 

DO 34 K=1,MM 

m= 1+! 

DO 4 J=IIoN 

IF(ABSF(A(JU9K) )-E1)49496 
TI=SORTFI (ACU 9K) ) #¥24+0 ACT 9K) )**2) 
S=A(J9K)/T1 

me( [1sK)/T I 


MO 5 L=KsMM 


T2=C#*A( ToL) +S¥*A (JU oL) 
A(JsL)=-S¥A(CT 9L)+C*#A(JoL) 
A(TeLl)=T2 

me) 1+!) 

SONTINUF 

BEC ABOF (A(T 9K) )—-E2)39398 
IL(K)=T 

T=[+1 

mONT TNUF 

X (MM) =-1.0 

my] =M 

mo 35 T=19M 

X(T)=020 

DO 30 J=loM 

IFC TL(11)) 39230931 
m=), 0 

me=1 1+) 

mey_(II) 

NO 32 K=LL»MM 
S=S+A(T 5K )#¥X(K) 
meTT)=—-S/ACTsITI) 
IT=IT-1 

IFC TL OMM) 50251950 
ALPHA=0,0 

me 10 52 

T=TL (MM) 

ALPHA=A(T »MM) 

RETURN 

END 





MAIN PROGRAM (COOLANT FLOW ANALYSIS) 


DIMENSION A(4)9B(6)9C(6)9D(6) 9E (6) sQUAL (3) sPIN(10) s0PNS(1093)sDPNM 
1(1093)s9TIN(10)9TAVNS(1093 ) »TAVNM( 1093) 9DTNS(1093) 9DTNM(1093) 9DENSI 
2N(10) 9VINNS( 1093) sVINNM( 1093) 9CVINNS(1093} »9CVINNM( 1093) s9VOUTNS(10¢9 
33)9VLONM(1093) 9VVONM(1093) sHVAVNS( 1093) sHVAVNM(1093) 9SLOTNS(1093) 9 
4SLOTNM(1093)s9CSLTNS(1093) »CSLTNM(1093) 

COMMON AsBeCeoD2FsQUuAL »VLIMsPOUT sDELTAP »sNOPINC®9PINSDPNS »sDPNMsTINo TA 
1VNS» TAVNMsDTNS»DINMs TOUT sDENSIN »DNLOUT »DNVOUT sVINNSsVINNMsCVINNS 9 C 
2VINNMs VOUTNS » VLONM 2 VVONMs HVAVNS s HVAVNM 9 SLOTNS 9 SLOTNM sCSLINS9CSLTNM 

marl RFANDS 

CALL NOSLIP 

CALL NOMIX 

meet INLFT 

CALL PRINTS 
moO 10 1 

FENN 
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SURROUTINF READS 


DIMFNSION A(4)9B(6)2C(6)9D(6)9E(6) sQUAL(3) sPIN(10) sDPNS(1093)sDPNM 
1(1093)sTIN(10)sTAVNS(1093)sTAVNM(1093) sDTNS(1083)sDTNM(1093)sDENSI 
9N(10) sVINNS(1093) sVINNM( 1093) sCVINNS(1093) sCVINNM(1093) sVOUTNS(105 
33)sVLONM(1093) sVVONM( 1093) sHVAVNS( 1093) sHVAVNM( 1093) 9SLOTNS(1093) 5 
4SLOTNM(1093)sCSLTNS(1093) sCSLTNM(1093) 

COMMON AsBsCsDsEsQUALsVLIMsPOUT sDELTAP sNOPINCs PINs DPPNS sDPNMoTINo TA 
LVNSsTAVNMesDINS sDTNMo TOUT »9DENSIN ea DNLOUT sDNVOUT sVINNSsVINNMsCVINNSC 
2VINNM s VOUTNS »s VLONM s VVONM 9 HVAVNS s HVAVNM sSLOTNS 9 SLOTN™Y sCSLINSsCSLINM 

PFAD ls (A(T)sI=194) 

RFAD le (R(T)sI=196) 

RFAD ls (C(I)sT=156) 

RFAD ls (D(1T)sTI=1:26) 

RFAD ls (E(T)sI=196) 

MmORMAT(5F14e6) 

RFAD 2s (QUAL(1I)sI=1]1 93) 

FORMAT(3F10.44) 

PFAN 3s VLIM 

FOPMAT(F1%.4) 

RETURN 

FNN 


~ 235 





SUBROUTINE NOSLIP 


DIMENSION A(4)s8(6)92C(6) 2D(6) 9E (6) sQUAL (3) 3PIN( 10) 9D PNS(1093) sDPNM 
1(10s3)9TIN(10) »9TAVNS (1053) sTAVNM(1093) 9DTNS(1093) 93DTNM(1093) sDENSI 
9N(10) sVINNS(1093) »VINNM( 1093) sCVINNS( 1093) s9CVINNM( 1093) 2VOUTNS(109 
33) sVLONM( 1093) sVVONM(10 93) sHVAVNS( 1093) 9HVAVNM(1093) sSLOTNS(1093) 9 
4SLOTNM(1093)9CSLTNS(1093) sCSLTNM(1093) 

COMMON AsReCeDsEsQUAL sVLIMsPOUT sDELTAP »sNOPINCsPIN»sDPNSsDPNMsTINge TA 
TVNSs TAVNMs DINSsDINMs TOUT sDENSIN»DNLOUT sDNVOUT »VINNSsVINNMSCVINNS 9 C 
2VINNMsVOUTNS » VLONMs VVONMs HVAVNS 9HVAVNM 9 SLOTNS 9 SLOTN™M 9 CSLTINSsCSLITNM 

RFAD 190s POUT »sDELTAP sNOPINC 

momMAT(2F10e49110) 

TOUT=TEMP (POUT ) 

mmeOoUT=1,./SVOLL(POUT) 

mmyvOUTH1,./SVOLV(POUT ) 
mo 7 T=lsNOPTNC 

Al=I 

BeNiC=AT*#DELTAO 

PIN( T)=POUT+P INC 

TIN(CI)=TEMP(PIN(T)) 

BYOLIN=SVOLL(PIN(T)) 

PENOIN(IT)=le/SVOLIN 

DO 6 M=l153 

SPVOUT=QUAL (M)/DNVOUT+4+(1¢-QUAL(M))/DNLOUT 

VINNS(T 9M) =SVOLIN*SQRTF(1201325E+06*PINC/ (SPVOUT=SVOLIN)) 

VOUTNS(T »sM)=HVINNS CI 9M)4+1201325E+06*PINC#SVOLIN/VINNS( TI 9M) 

IF (VLIM-VOUTNS(T9M)) 192922 

SEOTNS(1T9M)=0, 

oe) TO 6 

BON ST=,C08X¥QUAL(M)/(1LE/HIVP(PIN(I1))4+1-/HI VP (POUT) ) 

G=VINNS(1T9M)/SVOLIN 

m=O 

moUM=06 

HNVSUM=0, 

MTP OO=PIN(T ) 

DO 4 JHl1l925 

TSUM=TSUM+TFMP(POESS) 

HINV=ls/HIVP( PRESS) 

HNV SUM=HNVSUM4+HINV 

X=X+CONST#HINV 

SPVOLHEX¥SVOLV(PRESS)4+(1-¢=X )*SVOLL(PRESS) 

PRFSS=PIN(T)-G*(SPVOL#G-VINNS(JT 9M) )/1.201325E+06 

mONSTECONST#QUAL(M)/X 

QCHECK=ARSF((X—-QUAL (M) )/QUAL (M) ) 

TF( 2-OO01=QCHECK) 39595 

MAYVNS (1 »M)=TSUM/25. 

HVAVNS( IT »M)=HNVSUM/25, 

SLOTNS( I »M)=HVAVNS(T 9M) /(QUAL(M)#G) 

CONTINUF 
CONTINUF 

RETURN 

END 
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SURROUTINE NOMIX 


DIMENSION A(4)9R(6)9C(6) 9D(6) 9E (6) sQUAL (3) 9PIN( 10) sD PNS(1093)sDPNM 
1(1093)sTIN(19)sTAVNS(1053)9TAVNM(1093) sDINS(1093) sDTNM(1093) sDENSI 
9N(10) »VINNS(1093)sVINNM( 1093) sCVINNS( 1093) sCVINNM(1093) sVOUINS(109 
323) sVLONY (1093) 9 VVONM(10 93) sHVAVNS( 1093) sHVAVNM( 1093) 3SLOINS(1093) 9 
4SLOTNM(1093)sCSLTNS(1093)sCSLINM(1093) 

COMMON AsBsCsDsEsQUAL sVLIMsPOUT sDELTAP »sNOPINCs PINs DPNS »sDPNMsTINe TA 
TVNSsTAVNMsDINS »DINMs TOUT sDENSIN»s DNLOUT sDNVOUT sVINNSsVINNMSCVINNS 9 C 
2VINNMs VOUTNS » VLONM s VVONMs HVAVNS sHVAVNM sSLOTNSs SLOTNMY 9 CSLINS 9 CSLINM 

BO 13 J=lsNOPINC 

FV=2 -0265FH+06*#(PIN(I IT )-POUT)#AVSVV(IPIN( 1) sPOUT) 

FL =2 eO265FH+O06# (PINOT )-~POUT)*AVSVL(PIN( 1) sPOUT) 

mee lZ2 M=153 

VIN=VINNSCI 9M) 

ARATIO=(1e-QUAL(M) }*#DNVOUT/ (QUAL (M)*DNLOUT ) 

FIN=le/(DENSIN(I)* C1e+ARATIO)) 

DA=NDNLOUT#ARATIO 
VVONM(T 9M) =SQRTF(FV4+VIN#RVIN) 

VLONM( IT 9M) =SQRTF(CFL4+VIN#VIN) 
VINNM(ToM)=FINX (VLONMCT 9M) *DA+VVONM*DNVOUT ) 

DELTA=VINNM( TI »M)-VIN 

POHECK=ABSFI(DELTA)/VINNMCT 9M) 

VIN=VINNMCT 9M) 

TF(e0901=-VCHFCK) 19232 

IF (VLIM=VVONM(IT9M)) 329494 

BLOINM( 1»M)=06 

moO 10 12 

BONST=e04/SLOTNS(I »M) 

G=VINNM(T »M)*#DENSTIN(T) 

WFV=0¢ 
WFL=G 

TSUM=0. 

PREFSS=PIN(T) 

VLI=VINNMCT oM) 

VVI=VINNM(T oM) 

mo 8 J=1s25 

TSUM=TSUM+TEMP (PRESS) 

B= | e/SVOLL(PRESS) 

By=1./5VOLV(PRESS) 

H=HTVP( PRESS) 

WFL=WFL-CONST/H 
WFV=WFV+CONST/H 

BJ=26-J 
VL=VL1+(VLONM(T»M)-VL1)/AJ 

PP=WFL/(DL*VL) 

VV2=WEV/((1e-AL)#DV) 

VL2=SORTFI(DV¥(VV2*¥VV2—-VV1L¥VV1)/DL+VLI*¥VL1) 

WELCHEK=VL2-VL 

ME(VLCHEK) 6197362 

VL=VL4+VLCHEK/5e0 

mo !'0 63 

VL=VL+VLCHEK/10-20 937 





Met eDO1=—ABSF(VLCHEK)/VL2) 69797 
NDP=DV*¥(VV2#¥VV2-VV1¥VV1)/26 
PRFSS=PPFSS=$D0P/1.01325E+06 
mel=VL2 

VV1=VV2 
me=)V*(1.—-AL)/(DV¥(1.2-AL)+DL*¥AL) 
QCHFCK=X/QUAL(M)=-1l.e 

Mee QCHECK) 1091199 

meno? ~CONSTX( 1 e-QCHECK/10-)/(1e—-QCHECK/S,}} 
me} =CONST*(1.-QCHECK/5.) 

met eVO1-ABSF(QCHECK)) 5311911 
TAVNM(T 9M) =TSUM/25, 

SLOTNM( T9»M)=.e04/CONST 

CONTINUF 

CONTINUE 

RF TURN 

FND 
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SURROUTINF INLET 


NIMFNSTON A(4)93(16)9C(6) 9016) 9E(6) sQUAL (3) 9PIN( 10) 9DPNS(1093) »DPNM 


14(1093)9TIN(10)sTAVNS (1093) sTAVNM(1093) sDTNS(1093) 3D TNM(1093) 9DENSI 
2N(10) sVINNS(1093) sVINNM( 1093) »CVINNS (1093) sCVINNM( 1093) sVOUTNS(109 
33) 9VLONM( 1093) »VVONM(10 93 ) sHVAVNS (10 93) sHVAVNM(1093) 9SLOTNS(1093) 9 
4SLOTNM(1093) sCSLINS(1093) sCSLTINM(1093) 


COMMON AsBesCeDesEsGUALsVLIMsPOUT sDELTAP sNOPINC®sPINsDPNS sDPNMsTINoTA 


TVNSs TAVNMSDTINSsDINYs TOUT sDENSINsDNLOUT sDNVOUT » VINNSs VINNM3CVINNS 9 C 
2VINNM es VOUTNS » VLONM »s VVONMs HVAVNS sHVAVNM 9 SLCTNS 9SLOTNM s9CSLTNS sCSLTNM 


m1 )=1e 

MST D=TEMP (PSTN) 

DO 1 T=1lsNOPINC 

DH=HLIOQ(PIN( I) )-HLIQCPOUT )--20265#(PIN(I)-POUT)/DNLOUT 
TAVIN=(TEMP(PIN(T))+4TEMP( POUT) )/2. 

mo | M=15%3 

DLNS=DH*DENSINO I) #VINNS( IT 9M)*¥SLOTNS(T 9M) 

Mem iINS(J~sM)=SLOITNS (J »M)/(1.+DLNS) 

CVINNS(T oM)=VINNS( I oM)*DENSIN(I)/DNLOUT 
DPNS(Jo“)=PIN(IT)=POUT4+DENSIN(ITI*XVINNS(T 9M) *#O(VINNS(I9M)—CVINNS( I 9M) 


M7 1 e901 3225F+06 


— 


DINS( I »M)=(TAVNS(T 9M)+DLNS*TAVIN)/(12+DLNS)-TSTD 
DLNM=DH*DENSIN(T)*VINNMCT 9M) *SLOTNMCT 9M) 
CSLINM(T»M)=SLOTNM( JT 9M)/(1.+DLNM) 

CVINNM( T oM)=VINNM CT oM)*#DENSIN(O IT) /DNLOUT 
DPNM(ITsM)=PIN( I )=POUT4+DENSIN(IT)¥VINNMC I 9M) *#(VINNM(CI 9M) —-CVINNM(I 9M) 
Me 01325F+06 

DTINM( To“) =(TAVNMCT 9M)+DLNM*#¥TAVIN)Z(010¢+DLNM)-TSTD 
AA 5 T=lsNOPINC 

NO 5 M=17,3 

Met COLINS(I9M)) 39293 

my I NNS( IT «M)=0. 

DENS(T»M)y=%, 

MNTNS(TsM)=0, 

IF(CSLINM(I9M)) 59495 

CVINNM(T »M)=9, 

DONM( TT »M)=0, 

PINM( T»sMjy=HaN, 

CONTINUF 

RETURN 

FNA 
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BONCTION SVOLL(P) 


DIMENSION Al4)9B(6) 29C(6)9D(6) 9E(6) sQUAL (3) sPIN( 10) »DPNS(1093) »DPNM 
1(1093)sTIN(10)sTAVNS(1093)9TAVNM(1093) sDTINS(1093) sDTNM(1093) »DENSI 
9N(10) »VINNS(1093) »VINNM(1093) sCVINNS( 1093) sCVINNM(1093) »9VOUTNS(109 
43)s3VLONM(1093) 9 VVONM(10 92) sHVAVNS (1093) 9HVAVNM(1093) 29SLOTNS(1093) 9 
LSLOTNM(1093)sCSLTNS(1093) sCSLTNM(1093) 

COMMON AsBsCeoNsFEsQUALsVLIMsPOUT sPDELTAPsNOPINC sPINsDPNS»sDPNMsTINo TA 
TVNSsTAVNYsDINSSDINMs TOUT sDENSINsDNLOUT sDNVOUT sVINNSsVINNM9CVINNS 9C 
2VINNMs VOUTNS » VLONM 9 VVONM sHVAVNS sHVAVNM 9SLOTNS 9 SLOTNM 9 CSLINS 9 CSLITNM 

D?7=D#P 

P32=P?%P 

SVOLL=AA(1)4A02)*P4A(3)*¥P24A(4) #23 

PETURN 

ENN 
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FUNCTION SVOLV(P) 


DIMENSION A(4)9B(6)9C(6) 9D(6) 9E(§) sQUAL (3) sPIN( 10) sDPNS(1093) sDPNM 
nO 53) #1 1N( 10)» TAVNS (1053) »TAVNM(1093) »DTNS( 1093) »DTNM( 1093) »DENSI 
2N(10) sVINNS(1093) sVINNM( 1093) sCVINNS (1093) »CVINNM( 1093) sVOUTNS(10>5 
23)sVLONM( 1093 ) »VVONM( 10 93) sHVAVNS( 1093) sHVAVNM( 1093) sSLOTNS(1093) 9 
4SLOTNM(10932)2CSLTNS(1093)sCSLINM(1093) 

COMMON AsRsCeoDeEsQUAL »9VLIMsPOUT sDFLTAPsNOPINC SPIN» DPNSsDPNMsTINg TA 
IVNSe TAVNMsDTINSSDINMs TOUT »DENSIN®DNLOUT sDNVOUT sVINNSoVINNM9CVINNS9C 
IVINNMs VOUTNS » VLONM 9 VVONMs HVAVNS sHVAVNM 9 SLOTNS » SLOTNM9CSLINSsCSLTNM 

D2=PxXP 

mo =P2xP 

SVOLV=B(1)/P4+B(2)4+B(3)*P+B(4)*P24+B(5)*P34+B(6)*P2*P2 

RETURN 

END 
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FUNCTION TFMP(P) 


DIMENSION A(4)9B(6)2C(6)3D(6)2E(6) sQUAL (3) »PIN(10) sDPNS(10923)sDPNM 
1(1093)s9TIN(10)»TAVNS(1053)9TAVNM(1093) »DTNS(1093)9DTNM(1093) sDENSI 
9N(10) »VINNS(19093) sVINNM(1093) sCVINNS( 10923) sCVINNM(1093) sVOUINS(109 
33)sVLONM(1093) »VVONM(10 93) sHVAVNS( 1093) SHVAVNM(1093)2SLOTNS(1093) 9 
BSLOTNM( 1093) 9CSLTINS(1093) »CSLINM(1093) 

COMMON AsBsCoNesFEsQUAL »VLIMsPOUT »DELTAP sNOPINC® PINs DPNS*SDPNMsTINe TA 
IVNSs TAVNMsDINS 9 DINMs TOUT sDENSINsDNLOUT sDNVOUT sVINNSs VINNMSCVINNS9C 


2VINNMs VOUTNS s VLONM 9 VVONM 9 HVAVNS s HVAVNM 9 SLOTNS »@SLOTNMs CSLINS eo CSLINM 
O07=PxdD 


OD3=P7xD 

meMP=C(1)+C(2)*P4+C(3) *¥P24+C(4)¥P34+C(5)*P2*¥P2+C(6)¥P3#P2 
RETURN 

ENT 
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BONCTION HLIQ(P) 


DIMENSTON A(4)9B(6)9C(6)9D(6) 9E (6) sQUAL (3) sPIN( 10) s0PNS(1093)sDPNM 
1(1093)sTIN(19) »TAVNS(10 53) »TAVNM(1093) sDTINS(1093) 9DTNM(1093) 9DENS] 
9N(10) sVINNS(1093) 9VINNM( 1093) 9CVINNS( 1093) sCVINNM( 1093) sVOUINS(109 
33)9VLONM(1093) »VVONM(1093 ) sHVAVNS( 1093) sHVAVNM( 1083) 29SLOTNS(1093) 
4LSLOTNM(1093) xCSLTNS(1093) 9CSLTNM(1053) 

COMMON As89C9DsF2QUAL sVLIMsPOUT sDELTAP sNOPINCs»PINsDPNSsDPNMsTINgTA 
1VNSs TAVNMaDTINS »DINMs TOUT sDENSINsDNLOUT sDNVOUT sVINNSs VINNMs CVINNS 9 C 
2VINNMs VOUTNS sVLONM s VVONMs HVAVNS 9HVAVNMs SLOTNS 9 SLOTNMs CSLINSsCSLTNM 

BD2=P#D 

D2=P?x*P 

HLTQ=D(1)4D(2)*P4+D(3)*P24+D(4)*P34D(5)*P2*P24+D(6) *P3%P2 

RPETIJON 

FAD 
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FUNCTION HIVP(P) 


DIMENSION Al(4)98(6) 9C(6) 9D0(6) 9E (6) sQUAL (3) sPIN( 10) »9DPNS(1093)9DPNM 
1(1093)9TIN(10)sTAVNS(1093) 9 TAVNM(1093)9DTNS(1093) »DTNM(1093) sDENSI 
IN(19) 9 VINNS(1093) sVINNM(1093) »CVINNS (1093) sCVINNM(1093) 29VOUINS(109 
43)9VLONM( 1033) »VVONM(10 93) sHVAVNS( 1093) sHVAVNM( 1093) »SLOTNS(1093) 9 
ASLOTNM(1093)9CSLTNS(1093) sCSLINM(1093) 

COMMON AsRsCaD9F sQUAL sVLIMsPOUT sDELTAPsNOPINCsPINs DPNSsDPNMoTINe TA 
IVNS s TAVNMsDTINS eDINMos TOUT sDENSIN es DNLOUT sDNVOUT sVINNSs VINNMSCVINNS 9C 
IVINNM sVOUTNS 9 VLONM 9s VVONM 9 HVAVNS sHVAVNM 9 SLOTNS 9 SLOTNM 9 CSLINSsCSLINM 

P2=PxD 

pa-Pp7xP 

HTVP=F(1)4+E(2)*P+E(3)*P2+F(4)*934+E(5)*P2*#P2+E(6) #P3*P2 

PETURN 

END 
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mencCTION AVSVL(PHsPL) 


DIMFNSION Al4)2B(6)9C(6) 9006) 9E (6) sQUAL (3) sPIN(10) s3DPNS(1093)sDONM 
1(1023)seTIN(10)sTAVNS(1053 ) sTAVNM(1093) sDTNS(1093) 3DTNM(1093) sDENSI 
2N(10)sVINNS(1093) sVINNM( 1093) *CVINNS(1093) sCVINNM( 1093) sVOUTNS(10>2 
33) sVLONM(1093) »VVONM(1093) sHVAVNS( 1093) sHVAVNM(1093) s9SLOTNS(1093) 9 
m1 O1NM( 1053) sCSLINS(1093) sCSLINM(1093) 

COMMON AoBsCeDsEsQUAL sVLIMsPOUT sDELTAPsNOPINC®SPINsDPNSsDPNMsTINo TA 
1VNSsTAVNMsDINSSDTNMs TOUT sDENSIN sONLOUT sDNVOUT sVINNSs VINNMSCVINNS 9 C 
IVINNMs VOUTNS 9 VLONMs VVONMs HVAVNS s HVAVNM 9 SLOTNS 9 SLOTNM9 CSLINS sCSLTNM 

PH?=PH*PH 

Bea =PH2*OH 
PL2=PL¥PL 
DIL 2=P] 2¥PL 

AVSVL=(A(1)%(PH-PL)+A02)%(PH2-PL2)/2 6+A(3)*#(PH3—-PL3)/3e¢t+A(4)*(PH2* 
TPH2-PL2*PL2)/4e)/(PH-PL) 

PETUPN 

FAID 
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FUNCTION AVSVV(PH»sPL) 


DIMENSION A(4)9B(6)9C(6) 9016) 9E (6) sQUAL (3) »9PIN( 10) 90PNS(1093) sDPNM 
1(41093)sTIN(10)2TAVNS(10 93) sTAVNM(1093) »DTNS(1093) sD TNM(1093) »DENST 
2N(10) 2»VINNS(1093) sVINNM(1093) 9CVINNS( 1093) sCVINNM(1093) sVOUTNS(10>5 
33) sVLONM(1093)9VVONM( 1093) sHVAVNS(10 93) sHVAVNM(1093) 29SLOTNS(1093) 9 
&4SLOTNM( 1093) sCSLTNS(1093) »CSLITNM(1093) 

COMMON AsBeCeDsEsQUAL »VLIMsPOUT »DELTAP »sNOPINCSPINsDPNS sDPNMsTINo TA 
1VNS» TAVNMeDTINSsDINMs TOUT sDENSIN »DNLOUT »DONVOUT »VINNSs VINNM9CVINNS9C 
2VINNM s VOUTNS » VLONM 9» VVONMs HVAVNS sHVAVNM 9SLOTNSs SLOTNM a CSLINSa9CSLTNM 

PH?=PH*PH 

PH3=PH2*PH 

PLA=PL¥PL 

PL3=PL2#*PL 

AVSVV=(B(1)*#LOGF(PH/PL)4+B(2)*(PH=PL)4+B(3)%(PH2-PL2)/2e+B(4) *(PH3—P 
3 )/3e+B(5)*(PH2*PH2-PL2*PL2)/4,4B(6)* ( PH3#PH2—-PL3#PL2)/5e)/(PH=PL 
E) 

RETURN 

ENN 
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SURROUTINF PRINTS 


DIMENSION A(4)9B(6)9C(6) 9D(6) 9E (6) »QUAL(3) sPIN(10) sDPNS(1093) »DPNM 
1(1093)9TIN(10)sTAVNS(10 93) »TAVNM( 1093) sDTNS(1093) 3DTNM(1093) sDENS] 
2N(10) sVINNS( 1093) »VINNM( 1093) »CVINNS( 1093) sCVINNM( 1093) »VOUTNS(10¢ 
33) 9VLONM(1023) »VVONM(1093 ) sHVAVNS( 1093) sHVAVNM(1093)*SLOTNS(1093) 5 
4SLOTNM(1093)9CSLTNS( 1093) »CSLTNM(1093) 

COMMON AsB2CsD2ExsQUAL sVLIMsPOUT sDELTAP »sNOPINC®SPIN»DPNS*DPNMsTINg TA 
1VNS» TAVNMs DITNSsDTNMs TOUT sDENSIN»sDNLOUT sDNVOUT sVINNSs VINNMsCVINNS 9 C 
2VINNMs VOUTNS s VLONM » VVONMs HVAVNS sHVAVNM 9 SLOTNS » SLOTNMsCSLINSsCSLTNM 
PRINT 1 

FORMAT(1H1957Xs17HOUTLET PARAMETERS) 

PRINT 29 POUT »sTOUT»DNVOUT »s DNLOUT 

FORMAT(1HOx16HPRESSURE (ATM) =9F60397X921]HTEMPERATURE (DEG K) =9F6 
1e297TX922HVAPOR DENSITY (G/CC) =sF8e697X923HLIQUID DENSITY (G/CC) = 
2?9F8e5) 

PRINT 23 

FORPMAT(1HO/1H 936X259HFLOW PARAMETERS BASED ON NO=SLIP MODEL ARE T 
1ARULATFED B3FLOW.) 

mo 15 M=1 4.2 

PRINT 4% QUAL(M) 

FORMAT(1HO/1H 952Xs16HOUTLET QUALITY =9F4e2) 

PRINT 5 

FORMAT (1HO+20X s5HINLET 99X sSHINLET912X 94HMEAN 9 9X9 SHINLET »8X95HINLET 
199X »s6HOUTLET »6Xs LOHNORMALIZED) 

PRINT 6 

FORMAT(1H »18Xs8HPRESSURE s5X*1l1LHTEMPERATURE s5Xs1l1HTEMPERATURE 95X97 
JHDENSTTY »6Xs8HVELOCITYs6Xs8HVELOCITYs8Xs4HSLOT ) 

PRINT 7 

FORPMAT( 1H 920Xs5H( ATM) 9 8Xe7HI(DEG K)99X97HIDEG K) 28X96HI( G/CC) 6X9 8H 
1(CM/ SEC) s6X98HI(CM/SEC) 7X s6HHEIGHT/I1H ) 
mo 12 IT=1sNOPINC 

IF(SLOTNS(1T9M)) 1091398 
PRINT 9o PIN(T)*TIN(T)*TAVNS(I 9M) »DENSIN(T) »VINNS(T 9M) »VOUTNS(T oM) 
TsSLOTNS( J oM) 

FORMAT(1H 919X9F60398X9F6 oe 291OXs FH e2SBX9F 7 oD 94X9ELO 4H 4X9E100495X» 
MP1N.4) 

BO 10 1? 

PRINT lis PIN(T) 

FORMAT(1H 219X9F6e3938X 349H--=PROGRAM ERROR=-- COMPUTED SLOT HEIGH 
1T NEGATIVE) 

CONTINUF 
BO 10 15 
PRINT 14 
FORMAT(1H 937Xs58HHIGHER INLET PRESSURE RESULTS IN EXCESSIVE OUTLE 
Mm VELOCITY) 

CONTINUF 
PRINT 1 
PRINT 29 POUT »sTOUT sDNVOUT sDNLOUT 
PRINT 16 
FORMAT(1HO/1H »35Xs61HFLOW PARAMETERS BASED ON NO@=MIXING MODEL ARE 
PE TABULATED BELOW.) 

DO 25 M=193 _ 24a 





PRINT 45 QUAL(M) 

PRINT 17 

FORMAT (1HO911X95HINLET 99X s5HINLET912X s4HMEAN 9 9X sSHINLET 9 8X 9 SHINLET 
199X s6HLIQUIDs11X»5HVAPOR»> 8X» 1LOHNORMALIZED) 

PRINT 18 

FOPMAT(1H » 9XsB8HPRESSURE 9 5X s11IHTEMPERATURE 95Xs LIHTEMPERA TURE 95X97 
LHDFNSTTY »6Xs8HVFLOCITYs5Xs11HOUTLFET VEL. sS5Xes11HOUTLFT VEL. 98X%94HSL 
BOT) 

PRINT 19 

FORMAT( 1H 911X95H( ATM) s98Xs7HIDEG K)99X97THIDEG K) 98X96H( G/CC) 96X98H 
1(CM/ SEC) 96X98H(CM/SEC) 9 8X »8H(CM/ SEC) 9X s6HHEIGHT/1H ) 

BO 23 T=lsNOPINC 

IF(SLOTNM(T9M)) 22924920 

PRINT 219 PIN(I)sTIN¢O I) sTAVNM(T 9M) sDENSIN(I) sVINNM( IT 9M) »9VLONM(I 9M) 
LsVVONM( I 9M) *sSLOTNM(T 9M) 

“FOPMAT(IH 910X9F66398X9F6 02910X 9 F602 98X96 745 94X 9F1 004 95X9F1004 6X9 
1F102496X9E10 04) 

GQ TO 23 

PRINT lls PIN(T) 

CONT INUF 

mee 1 25 

PRINT 14 

CONTINUF 

PRINT 26 

FORMAT(1H1942Xs48HFLOW PARAMETERS ADJUSTED FOR SUBCOOLING AT INLET 
i ) 

iy 31 M=152 

PRINT 27sPOUT sQUAL (M) sDNLOUT 

FOOMAT(IHO/1H 912Xs17HOUTLET PRESSURE =9F6.395H ATMes1l1Xs1l6HOUTLET 
BUALITY =9F4.2911X»24HLIQUID DENSITY AT PUMP =9F7.5:5H G/CC) 
PRINT 28 


FODPMAT(1H0510X 951 H------------ e+ NQ SLIP MODEL -+--------- 
fas 1 OX 9 51 H-—--+2-_-___--_ NO MIXING MODEL =--===—==——=—— ) 
PRINT 29 


FORMAT(C1LHOs1OXs9HDELTA (P)*5Xs9HDELTA (T)95X*91OHNORMALIZED 96X s6HV 

LOIN) sLIXsS9OHDELTA (P)25Xs9HDELTA (T)95X91OHNORMALIZED96X96HV (IN)) 

PRINT 30 

FOPMAT(1H 913X93HATM910X95HDEG Ks8Xs8HSLOT HT es 7X96HCM/SEC914X93HA 
1TMs1O0Xs5HDEG Ks8Xs8HSLOT HT es7Xs6HCM/SEC/I1H ) 

PRINT 3229 (DPNS( 79M) »DINS( 19M) sCSLINS(I 9M) sCVINNS(T 9M) »DPNM( I 9M) 9D 
LTNM( I 9M) sCSLTNM( I 9M) s9CVINNM(I 9M) sIT=l1 »NOPINC) 

FORMAT IH 912X9F5_399X9F5 2 96X9ELO e449 4XsFlLOe4s12X9F5 0 399X 9F 502 96X» 
1F10.494X9F10,4) 

RETURN 

FENN 
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